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Theory of Optimal Signal Filtration 
in the Presence of Internal Noise 


in Systems With Variable Parameters 


THER papers (1)! have been devoted to the synthesis 

of an optimal correction device for an automatic con- 
trol system (ACS) in the presence of stationary random 
perturbations (or internal noise) applied to various points of 
the system. 

This paper gives a complete solution of the problem of the 
synthesis of an ACS correction device for the class of systems 
with variable parameters in the presence of nonstationary 
random perturbations applied to various points of the svstem. 


1. Formulation of the Problem 


We shall assume that the ACS is subject to the interference 
n(t) and the internal random noise u(t), . . . , Un(t) applied 
to various points of the system, and that at a certain moment 
ty the input of the system is supplied with a control or useful 
signal s(t) consisting of a nonrandom component g(t) and a 
random component m(é), i.e., s(t) = g(t) + m(t). See Fig. 1. 

Thus, a signal y(t) = s(t) + n(t) is applied to the input of 
the system. In the presence of interference and internal 
random noise, let the system, for which we are seeking an 
optimal correction device, reproduce in the best way the 
signal 


s*(t) = H(p,t)s(t) = | x(t,a)s(a)do [1.1] 


(p = d/at) 


where H(p,t) is the reproduction operator, and x(t,7) is the 
weight function of the operator H(p,t), i.e., H(p,06(t — 7) = 


Translated from /zvestiia Akademii Nauk SSSR, Energetika 
i Avtomatika (Bulletin of the Academy of Sciences USSR, Power 
and Automation), 1960, no. 2, pp. 137-152. Translated by Re- 
search Information Service, New York. 

1 Numbers in parentheses indicate References at end of paper. 
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x(i,7). By “best” is understood the signal reproduction for 
which the following optimization conditions are fulfilled: 


Ele(t)] = 0 [for any form of nonrandom function g(t)] [1.2] 


E} gle(t)|} = min (at any moment of time ¢) [1.3] 
where 
e(t) = s*(@) — [1.4] 


w(t) is a process at the output of the system; E[e(é)] is the 
mathematical expectation of e(¢), and g(é) is some definite 
non-negative function of ¢. 

We shall assume further that: 

1 m(t), n(t), w(t), w(t), and u,(t) are un- 
correlated random (in general, nonstationary) functions with 
zero mean values and known co-variation functions. 

2 The nonrandom function g(t) has the form 


v=0 


where ¢, (v = 0,1,.... 2) are known constants, { ¢,(t)} is a 
certain arbitrary (known) complete system of functions. 

3 All the operators describing the behavior of part of the 
controlled member W;(p,t) (¢ = 1,2,..., m) are known. 

On the basis of the weight function of the whole system. 
found from the solution of this problem, it is possible to deter- 
mine the weight function of the optimal correction device. 


2. Derivation of the Integral Equation 


A. Case in Which the Optimal Correction Is Effected 
in Parallel and in Feedback, i.e., When W (p,t) Is Known 
and W (p,t) Is Unknown 

Before deriving the basic integral equation, we shall con- 
sider the question of reducing this problem to the problem of 
optimal filtration of a signal, when all the influences are ap- 
plied only to the input of the system. 


u,(t) v,/t) Up lt) 
«| | °°» ° | 
“t,r) K,(t,t) MK, (t,t) Min 


Fig. 1 4—Controlled member; B—known series correction; C—parallel correction 
in feedback 
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It is easy to see that the output 2(¢) of the real system, 


illustrated in Fig. 1, is determined by the two relationships 


c(t) = | I | W (p,te(t) + 
t=] 
| Il {2.1] 
v=1Li=p 
eo(t) = g(t) + m(t) + n(t) — W@(p,t)x(t) [2.2] 
where 


mM 
| = Walp) W [2.3] 
t=j 


(9 S41, 2 


Note here that, generally speaking, for operators with 
variable coefficients there is no commutation, and, accord- 
ingly, in Eqs. 2.1 and 2.3 the asterisk is used to emphasize the 
completely definite order of action of the operators on the 
function. 

From the operator equations 2.1 and 2.2 we get an expres- 
sion for x(t) 


= 


| 


I Wi(p,t) 


i=1 


| W (pt) W © (p,t) 


1 Mx 
rt) = | II | Wi(p.t) 
i=1 i=1 


or 


{1 — H(p,t)][g(t) + m(t)] + + Il 


Introducing the following new functions 


M(t) = s*(t) = H(p,t)[g(t) + m(t)] [2.6] 


= [1 H(p,t)|s(0) + + [2.7] 


v=1 


Ut) = {| wan ti W(p,t) 
t=1 
[2.8] 


(v 
we shall write Eq. 2.5 in the form 


my 
a(t) = +] IT Wilp,t) | WO (p,t) 


i=1 


| | WOp)IMO + NO] [2.9] 
t=1 


Je 


Note that 1/() and N(t) are completely definite and known 
functions, since the operators W;(p,t) and 
H(p,t) are given, and W“(p,t) is considered known for this 
case. 

Then, in accordance with Eqs. 1.4, 2.6 and 2.9, the repro- 
duction error is 


II W.(p,t) 


t=1 


e*(t) = M*(t) — x*() = M*(t) 1 + | 


e(t) = 


- {1 + | 
i=1 


t=] 


W (p,t) 4 


| [MO + [2.10] 
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t) 
x 
A*(t,T) 


wp t) 


Fig. 2 


| W (pt) [g(t) + m(t) + n(t)] + 


f 


m 


v=] 


Mx 


IT Wi(p,t) 


i=pv 


my 


IL Witp,t) 


| WO (pb) 
i=l 


| [2.4] 


ta + m(t)|] + 


—] 
| | | 


Let us consider the tracking system, shown in Fig. 2, which 
is subjected to the influence of a useful signal 7 *(t) and the 
interference N*(t). 

It follows from Fig. 2 that 


m 


v=1 


W(p,0) 


Mx 


II Witp,t) 


i=1 


[2.5] 


Mx 


r*(t) 1 + W(p.t) 
t=1 


Mx 


Let this tracking system be the best (in the sense of fulfilling 
the criteria of optimization, Eqs. 1.2 and 1.3) way of tracking 
the useful signal J/*(¢). 

Then, by definition, the reproduction error for this problem 
e*(t) is 


| + [2.11] 


Mz 


1 
| 
t=] 


(pt) + N*(O] (2.12) 


Comparing Eqs. 2.10 and 2.12 we see that if J/*(¢) = WW/(6) 
and N*(t) = N(t), then e*(t) = e(t) and 2*(t) = x(d). 

This means that the problem of optimal filtration of the 
signal s(t) in the presence of internal noise w(é), ... 
u,(t) and the interference n(é) has been reduced to the prob- 
lem of the optimal filtration of a signal ./(¢) in the presence 
of the interference N(t) (tio internal noise), i.e., constitutes 
the ordinary problem of optimal filtration of a signal. 

Henceforth, to simplify the calculations, we shall take 
g(t) =t%. It follows from Fig. 2 that 


[2.13] 
1307 


x(t) = [Md + NOI 


* ¥ * 
(2) 
A“ (t,t) 
= 
: 


where 


P(p,t) = + | Il | W (p,t) x 


t=1 


i=1 


| Tl | [2.14] 


MW) + NW = git) +m) + (2.15) 


v=1 
We shall derive the basic integral equation for the case in 
which the optimal system sought is a system with a finite 


memory 7’. 
Assuming that the correlation function of the processes 


mi) + nt) + 


v=l1 


has 2q¢ continuous derivatives, and, hence, that we differenti- 


Integrating Eq. 2.16 by parts, interpreting the substi- 
tution of limits during integration in the following way 


| ig(t,r) | 
T «—0 T 
r=(-—T r=t—T+e 
7) | 7) 
Or? 1=t Or! T=t—e 
(e>0) 
and taking into account the relationships 
(So) 
iol | 
lim = f glt,r)6O(t — r)dr 
(e>0) 


we obtain from Eq. 2.16 the following expression for x(t) 


x(t) = an + m(r) + n(r) + Ute) [2.19] 
v=1 
where 
q q 
+ — T — 7) + — 7) 
Kir) = i(t,7) (t)d™( T) (v)\(t — 7 re[ 12.20} 
0 t= 
t — T,t 
0 7 (t— 
The functions A,(t) and B,(t) (v = 0,1, ..., g) characterize 


ate the process 


m 
g(t) + m(t) + n(t) + Ub) 
v=l1 
q times, we shall make use of the following random process 
transformation formula (2,3) 


di 
= — 1 g(r) + m(r) + n(r) + 
v=1 


where W,(t,7) (j = 0,1, ..., q) are functions with limited 
variation characterizing the transformation operator ®(p,t). 
It is assumed that the functions W7,(t,7) (7 = 0,1, . . . , q) are 
continuous in the interval t — 7 < 7 <t, and, possibly, have 
jumps of the first class when 7 = t — T and 7 = ¢, so that the 
Stieltjes integrals in Eq. 2.16 exist. 


possible jumps of the weight function of the system when 
7 =t— Tand7 =:t. Replacing M*(t) and N*(t) with 
M(t) and N(t) in Eq. 2.12 and taking into account Eqs. 
1.1 and 2.19, we obtain from Eq. 2.12 the following expression 
for the reproduction error in the reduced problem 


e(t) x(t,7)[g(7) + m(r)]dr — 


K(t,7) Ke + m(r) + n(7) + Ute) [2.22] 


v=1 
Now by satisfying condition 1.2 or by requiring that the 
dynamic error be equal to zero, we get the following expression 
from Eq. 2.22 


»=0 


S(t} 2) w, (2,t) Wr X,(t) 
Fig. 3 
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sought 


T K(t,7) ¢(r)dr f x(t,7) ¢g(r)dr 


From Eq. 2.23 we get n + 1 instantaneous equations, which must be satisfied by the weight function of the optimal system 
= 0 
When Eq. 2.24 holds, formula 2.22 takes the form 


UU 


= f x(t,r)m(r)dr — K(t,7) | 


From Eq. 2.25 it is easy to get the following expression for the mean square reproduction error 


[2.24] 
m(r) + n(r) + vio) | dr [2.25] 
v=1 
v=1 
ditions (instantaneous equations) 2.24. 


Since the problem is a variational problem of the conditional extremum of the functional 2.26, then, following Lagrange’s 
n 

I[K(t7)] = — 2 Axe) 

where A,(t) (7 = 0,1, .. 


+o 4 
| —2 f [2.26] 
in t-—T 
method, we must solve the variational problem of the absolute extremum of the functional 
730 


The problem consists in finding the function A(¢,7), which leads to the minimum of the functional 2.26 and satisfies the con- 


, n) are indeterminate Lagrangian multipliers. 
criterion of the optimum 1.3, when g(t) = @, the basic integral equation 
t- 


K(t,o)[Rn(o,7) + R,(0,7) + R.(o,7) |do 
where 


i... 


n 


[2.2 
Solving the variational problem of the absolute ex- 


equation 


tremum for the functional 2.27 by the usual method, we get, as the necessary and sufficient condition for the fulfillment of the 
v=0 


(t(-T 
R.(t7r) = >> 
v=1 


and the n + 1 instantaneous equations 


where 


f K(t,o)[Rn(o,r) + R,(0,7) + R.(0,7) |do = f 


optimal system sought is a system with an infinite memory, then the optimal weight function A(t,7) must satisfy the integral 


[2.28] 
[2.29] 
By a completely analogous method it is possible to show that if we derive the integral equation for the case in which the 


K(t,7) ¢(r)dr — f )dr 


v=0 


(— © 


<7 [2.30] 

“0 n) [2.31] 

q 
K(t,r) = < pad 
Io 
q 

K,(t,7) = 
B. 


when 7 (—®, 


[2.32] 
when 7 € ¢) 
when 7 ¢ (—@, 0) 
Case in Which the Optimum Correction Is Effected in Series, i.e., When W (p,t) = land W'” (p,t) Is Unknown 
In this case we get an expression for x;(¢) from Eq. 2.4 
{ 
i=1 


My 


wars | Il 


i=1 


[2.33] 
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| WO(pt)[gQ) + + + 


i=1 


Il Wi(p,t) 


| 3 


W 
v=] Li=p 


fu [2.34] 
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In accordance with Eqs. 1.4 and 2.34 the reproduction error ¢,(¢) is 
mx 
a(t) = s*(t) — a(t) = {! + | II we | W M(t) - 
i=l 
{! ii oo | W(pQ)Ni(t) [2.35] 


i=1 i=1 
where 
MW) = + mp) — Vid (2.36] 
k=1 
Ni(t) = + [1 — + [2.37 } 
= Jue [2.38 | 
i=k 


Let us consider a tracking system, influenced by the useful signal ./;*(¢) and the interference V,*(t). From Fig. 4 it is easy 


to get the output of the tracking system 


4 Mx my 
a*(t) = ! | Il | Il | WO (p,t) + [2.39] 


i=1 Li=1 


Let the tracking system track the signal 1,*(¢) in the best way (in the sense of fulfilling conditions 1.2 and 2.3). 
Then, in accordance with Eqs. 1.4 and 2.39, the reproduction error is 


( my —1 


mx 
Jy | Ww | E | WO(pONi*(t) [2.40] 


\ i=1 i=1 
Since 
mx Mx 
{! + | II ti | WO (pt) *(t) = 
i=1 i=1 
+ | oo M,*(}) [2.41] 
i=1 
by comparing Eqs. 2.35 and 2.40, it is easy to see that e.*(¢) = e(é) and «,*(t) = 2,(0), if 


M*() = M0 = [2.42] 


This shows that in this case too the problem of optimal filtration of a signal s(t) in the presence of internal noise w(é), ue(t), 

. , Um(t) and interference n(t) can be reduced to the ordinary case of optimal filtration of a signal ./,(é) in the presence 
of the interference N;(t) only. 

In this case too, by applying the method of obtaining the integral equation outlined previously to the reduced problem of the 
optimal filtration of a signal, we get the following integral and instantancous equations: 

(a) In the class of systems with a finite memory 7 


+o 
K(t,o)[Rn(o,7) + R,(0,7) + R,(0,7) \do = + + p(t) 
t v=0 


(2.43 ] 
K(t,r) ¢r)dr — ¢Ar)dr = 0 (y = 0; 1,...,; 2) [2.44] 
where 
R(t,7) = [2.45] 
k=l 
(v = 0,1, ..., m)—indeterminate Lagrangian multipliers. 
(b) In the class of systems with an infinite memory 
© n 
+ R,,(¢,7) + R,(0,7) |do = x(t,0)R»(a,7)do + R,(t,7) + u(t) ¢,(7) (-— fo] < t) 
[2.46] 
t + 
fi K(t,r)¢(r)dr — f x(tr)¢(r)dr =0 =0,1,...,2) [2.47] 
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It is necessary to point out that the integral equations 2.25 
and 2.43 are equations of the same type, and therefore the 
method of solving these two integral equations is the same. 
Analogously, we may note that the integral equations 2.30 
and 2.46 are of the same type, and, accordingly, exactly the 
same method can be used to solve both of them. 


3. Solution of Integral Equations 2.28 and 2.30 


4. Method of Solving Equations 2.28 and 2.30 for a 
Random Process Reducible to a Stationary One (4) 


Let the random function 


m(t) + n(t) + 
k=1 


be a random function reducible to a stationary one (5), Le. 


m 


+ Ut) = [3.1] 
k=1 


m(t) + n(t) 


where f(t) is a nonrandom function, not having zeros in the 
interval [¢ — 7,t], and u(¢) is a stationary random function 
with a known co variation function R,,(7) [or with the spectral 
density S.(w)]. Then we have the equality 

m 


R,,(t7) + Ri(t,7) + R(t) 
k=1 


= f(Of(r)Ru(t — 7) [3.2] 


Bearing in mind Eq. 3.2, it is possible to rewrite Eq. 2.28 
1 


thus 


K(t,o)f(o)R.(7 — = 
t-— 


l 
+ v(t) 
f(r) x f(r) 


@-—-T<r<t) 


[3.3] 


Let the stationary process u(f) have a fractional-rational 
spectral density 


* am. 2 
Fig. 4 


Then it is possible to show (2) that the correlation of the 
stationary random process R(t — 7) with a spectral density 
of the form 3.4 is connected with Green’s function of the self- 
adjoint boundary value problem 


Lip) L¥(pyy(t) = 0 


lim y(t) = lim y(t) = 0 


t+ « 
by the relationship 
R(t — = M(p).M*(pyG(t 7) [3.6] 
where G(t — 7) is Green’s function of problem 35 and the 
solution of the equation 
L(p)L*(p)G(t — 7) = 6(t — 7) Br] 


Taking Eq. 3.6 into account, it is easy to get the following 
differential equation with constant coefficients from Eq. 3.3 


M,(p).M,*(p) K (ta) — ove | = 


+ gr) 


<r <0 
(jw) 
v=0 VE *( 7@) where 
aw” d d 
M(p)M,*(p) = M* 
dr dr 
Its general solution may be represented thus 
t e4 w R (r) | 
K(t,a)f(a)G(r — a)do = p) J .o)R,, PATIL 
(t,o), r a)do + p) + f(r) { 13.9] 
where Ay (uw = 1,2,..., 2k) are roots of the equation = 0, and (a = 1,2, ..., 24) are arbitrary constants 
relating to 7 and are subject to determination. 
Applying the operator L,(p)L,*(p) to both sides of Eq. 3.9 and taking into account Eqs. 2.20 and 3.7, we get 
Xu? + « 
f(r) f(r) 
= f(r) f v=0 v=0 
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In order to determine D,(t) (u = 1, 2,..., 2k), A,(t) 
(vy = 0,1,..., 9), (v = 0, 1, ..., g) and (» = 0, 
1,..., 7), we substitute Eq. 3.10 in Eq. 2.28 and require that 
it be fulfilled unconditionally with respect to 7. We then 
get 2l linear algebraic equations for determining the un- 
knowns. The remaining n + 1 equations are obtained from 
the instantaneous equations 2.24. 

By solving the system of 2/ + n + 1 algebraic equations, 
it is possible to determine the still unknown functions D,(¢), 
A,(t), B,(t) and A,(t). In a quite analogous way we get a 
solution of Eq. 2.30 in the form 


K 2k D l 
(t.7) = KO + 


flr) 


In order to determine D,(t) (u = 1,2,..., 2k), c(t) (v = 
0,1,...,q) and A,(t) (v = 0,1, ... , n), we substitute Eq. 
3.11 in Eq. 2.30 and require that it be fulfilled unconditionally 
with respect to 7. We then get / + k algebraic equations for 
determining the unknowns. 

The remaining n + 1 equations are obtained from Eq. 
2.31. By solving the system of 1 + n + k + 1 equations, it 
is possible to determine the unknown D,/(t), c,(t) and X,(¢). 


B. Generalized Shinbrot Method (8-11) of Solving 
the Integral Equations 2.28 and 2.30 


Since, according to Eq. 1.5, the functions ¢,(t), ge(t),..., 
¢,(t), .. . form a complete system of functions, it is possible 
to get a sufficiently good approximation of any continuous 
function of two variables with the aid of a linear combination 
of the functions | ¢,(t)}. Therefore we shall henceforth use 
the following representations of the co-variation functions 


pi 
= alOeAl7) = alter) 
v=0 v=0 
pe 
R.(t7) = = Yb (t 2 7) 
v=0 v=0 
ps Pp 
v=0 v=0 
+o 
f = > = 
v=0 
Pp 
(t 2 7) 
v=0 
where 


Pp = max {pi, po, ps, ps n} 
a,(t) =0 when p > v > p; 
=0 when p > v> pp 
c,(t) =0 when p > v > ps 
d,(t) = 0 when p > v > ps 


The quantities p; (7 = 1,2,3,4) depend on the form of the 
¢o-variation function and the accuracy of the approximation 
(in the case in which we replace the kernel with a close approxi- 
mate kernel). 

Then, on the basis of Eqs. 3.12, it is possible to rewrite 
Eq.'2.28 in the form 


| K (t,o) [E(r)®(o) — E(o)®(r)|do [3.13] 
((-T<7r<t) 
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L(p)Ler*(p) Vie) » 


q 
+ e(t)d(t — 7) 


where 

F(t) = Dit) + [3.14] 
E(t) = A(t) + BY) + C(t) [3.15] 
A(t) = fao(t), a(t), ..., ap(0)} 

B(t) = {bo(t), bi(t), ... , b,(t)} 

C(t) = {co(t), c(t), .. . , 

Dt) = }do(t), d,(t) d,(t)} 

®(t) = { got), gilt), ..., ¢p(t)} 


+ 


(-~<r¢tq=l—k-1) [3.11] 


We shall seek a solution of Eq. 3.13 in the form 
K(t,r) = g(tr)u(t — + — 7) (3.17 ] 
where 
when < 0 [3.18] 
1 when ¢ > 0 
g(t,7) = Owhens <t—T [3.19] 


Then from Egs. 3.13 and 3.17 we get the equation 


t—T 


[E(7)®(c) — E(o)@(r) <7 <t) [8.20] 
where 

G(t) = F(t) — A(HE(t) [3.21] 

Now we shall try to find g(¢,7) in the form of the scalar 


product of two vectors in a p + 1 dimensional space, i.e., in 
the form 


g(t.7) = [3.22 
where 
M(r) = Owhens <t—T [3.23] 
= .-., 7 p(t)} 


M(r) = {mo(r),mi(r), ...,m,(7)} 
Then from Eq. 3.20 we get the two integral equations 


— E(o)®(r) do = (3,25) 
t-— 
<r <b 


First we shall determine ./(7) from Eq. 3.25. 

It is not difficult to show that m,(7) (A = 0,..., p) is the 
essence of the solution of the following differential equations 
(8) 


= falr)u(r — + — &) [8.26] 
»=0 


(p= 
dr 


with the initial conditions 
= 0 [3.27] 
s—1;A\ =0)1,...,p) 


where L,(D,r) =a,,(7) D* + ag—1,,(7) +... + + 
a (7), aj(7) are known functions of 7 and c¢,, are known 
constants. 
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It is possible to show (12) that the differential equations 
3.26 with the initial conditions 3.27 give the same solutions, 
as are given by the following differential equations 


= fr(r) (r > &) [3.28] 
with the initial conditions 
= MralE*) [3.29] 


if in Eq. 3.29 mya(E*) is determined from the relationships 


(—1)*cka = 
n=0 dq 


1 


[3.30] 


where 


Ao*(q,€) = ain(&*) — + q’aza(é*).. 

— gazr(&*) + g?aan(E*) 

+ 


A,*(q,&) 


Although the differential equations 3.26 with the initial 
conditions 3.27 cannot be solved on analog or calculating 
machines, the differential equations 3.28 with the initial 
conditions 3.29 can. We shall now determine the vector 
I(t) from Kq. 3.24. From the vector equation 3.24 we get 


gt) = y(t) + > [3.31] 
= 


J= 


Substituting the functions m,(7) thus found in Eqs. 3.31 
and solving the system of equations 3.31, we completely de- 


l 


(t,7) dtor [A (t,7) K*(t,7)|] + 


7 ot 


Since G(t) depends on A(é) and h(t), the expression y;(0) 
contains the still unknown functions Ao(t), A(t), 
and A(t). The values of Ao(t), Ai(t), .. may be deter- 
mined from the instantaneous equations 2.24. 

We shall determine h(é) from the requirement that the ex- 
pression for A(é,7) contain derivatives of the Dirac function 
of the lowest possible order. 

This completes the process of determining A(¢,7). This 
method makes it possible to make an exhaustive computation 
of the weight function K(t,7). With this method it is also 
possible to solve the integral equation 2.30 (from the in- 
stantaneous equations 2.31). Note that this method holds 
when all the co-variance functions can be represented in the 
form of a sum of the separate functions. When the random 
process is white noise, the method cannot be applied in this 
form to solve equation 2.28. A method of solving 2.28 in this 
case is discussed in (11). 


4. Determination of the Weight Function 
K@(t,7) With K(t,7) Known 


We shall consider the problem of finding the weight func- 
tion of the optimal correction system AK (t,7), when the cor- 
rection is effected in parallel and in feedback (Fig. 1). 
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\ (7,7) 


f = K® — 


It is natural to assume that all the weight functions of every 
link of the controlled member K,(t,7), Ke(t,7), ..., Km(t7), 
the weight functions of the series correction A“(t,7), and 
the weight function of the optimal system A(t,7) are known. 
It is required to find the weight function AK (t,7). 

Since by definition the weight function of the system 
K(t,r) is the response of the system at the moment ¢ to a unit 
impulse, applied at the moment 7, we have the relationships 


K(tr) = fi 


e(t7) = 6(t — K® (t,r2)K (t2,7)dt2 [4.2] 


[4.1] 


where 


From Eqs. 4.1 and 4.2 we get an integral equation for de- 
termining the weight function A (t,7) 
K(t,7) = K*(t,7r) — 
It is necessary to note that in deriving Eqs. 4.2 to 4.4 we 
made use of the conditions of physical realizability 


K(t,r) = Ofort <r 


K,(t,7) = K,(t,7) 


We shall transform the integral equation 4.4 into an equa- 
tion for which it is possible to apply the method of successive 
approximations to find a solution of the integral equation. 

We shall first differentiate with respect to 7 and then with 
respect to ¢, finally obtaining the differential equation 


t OK 
t 0K *(t,c) rOA(K\,T) 


Iq. 4.5 permits the use of the method of successive approxi- 
mations, and, consequently, the problem of finding the weight 
function A“ (t,7) may be solved completely. 

We shall now consider the question of finding the weight 
function of the optimal correction system A“ (t,7), when the 
correction is effected in series (Fig. 3). We assume that 
all the weight functions of every link in the controlled member 
K,,(t,7) and the weight function of the 
optimal system A(¢,7) are known, and that it is required to 
find the weight function K“(t,7). 

As in the previous case, we have the relationship 


K(t,r) = [4.6] 
éo(t,7) = — — K(t7) [4.7] 

where 
= fi K“ (t,o) K**(0,7)do [4.8] 
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From Eqs. 4.6 and 4.7 we get an integral equation for 
determining A“ (t,7) 


K(t.r) = — 


T 


K f [4.10] 


Differentiating Eq. 4.10 with respect to 7, we get 


KOEr) = 


K(t,a)do 


hq. 4.11 is an integral equation for which it is permissible to 
use the method of successive approximations, and, conse- 
quently, the problem of finding A“(t,7) may be completely 

solved. 
—Submitted Dec. 7, 1959 
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Reviewer’s Comment 


This paper represents a comprehensive treatment of the 
synthesis of optimal compensation for automatic control 
systems operating in the presence of nonstationary random 
perturbations applied to various points of the system. The 
author treats a very general class of systems initially, and 
succeeds in demonstrating that this class can always be re- 
duced to a class of simpler form where all random perturba- 
tions are applied at the system input. The author then pro- 
ceeds to detail a method for determining the time-variable 
impulse response of the compensating element which will 
minimize the mean square error of the system response sub- 
ject to the constraints of physical realizability, fixed settling 
time and zero error in response to a specified nonrandom 
component of the input signal. Techniques are presented 
for determining the impulse response of the compensating 
element for the case where this clement is introduced in the 
system feedback path and also for the case where the com- 
pensation is performed in series with the fixed element in the 
feedforward path. 

The procedure indicated by the author involves three basic 
steps: The reduction of the system to a single input system, 
the determination of the optimum system impulse response, 
and the determination of the impulse response of the system 
compensator. As might be expected, the most complex step 
is that of determining the optimum system impulse response. 
Here the author develops an integral equation for the opti- 
mum time-variable system which represents an extension of 
Booton’s (1) work. The extension consists of including the 
constraints of finite system settling time and zero error in 
response to nonrandom input signal components. These 
constraints are included by the author in a manner similar to 
that used by Ragazzini and Zadeh (2) in their treatment of 
tracking system optimization. Two methods are described 
for solution of the integral equations defining the optimum 
system impulse response. The first method discussed em- 
ploys spectral factorization (3,4) and can be used only when 
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the random process involved can be reduced to a stationary 
process. The second method discussed represents an ex- 
tension of Shinbrot’s (5) method. 

The value of this paper would be increased if the author 
had included examples to illustrate the application of the 
somewhat complex mathematical processes discussed. This 
lack combined with a lack of complete rigor makes the dis- 
cussion of the generalized Shinbrot method of solution diff- 
cult to follow. Such difficulty arises in the extension of the 
Shinbrot method to the case where a finite settling time con- 
straint is involved. For example, the author in Eqs. 3.22 of 
the text writes g(t,7) as the vector product P'(t)M(r) where it 
is implied that M(r) is independent of t. However, Eqs. 
3.23 restrict M(r) to equal zero for 7 < t — T, thus requiring 
that M(r) be also dependent on t. This dependence of M(z) 
on ¢ raises a question concerning the general adequacy of the 
author’s indicated method of solution for T(t) and M(z): 
however, it appears that such a solution is obtainable from 
a direct extension of Shinbrot’s approach as explained in 
Shinbrot’s paper. Any dependence of M(zr) on t will be in 
the form of step functions, impulse functions, and derivatives 
of the impulse function all occurring at 7 = t — T. 


-N. W. TREMBATH 
Control Systems Department 
Space Technology Laboratories, Inc. 
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Stability of an Orthotropic Cylindrical 


Shell With Axial Load and 


External Pressure 


HIS paper deals with the stability of an orthotropic closed 

circular cylindrical shell acted on simultaneously by ex- 
ternal transverse pressure and axial compression in both its 
linear and nonlinear formulation. A method of computing 
the critical load and representing the results in graphical form 
is proposed. Graphs to facilitate computation are given for 
an orthotropic shell. 


Nomenclature 


L,Yi2 = coordinate axes along generatrix of middle surface 
of the shell, along are of centerline of shell cross 
section, and along normal to the middle sur- 
face in direction of center of curvature, respec- 


tively 

R = radius of curvature of the middle surface of the 
shell 

a = length of shell along the generatrix 

h = thickness of shell 

H = a’/Rh, geometrical parameter of the shell 

s = number of half waves along generatrix 

n = number of whole waves along are 

¢ = # = (rhs/an)*, square of ratio of length of a 
half wave along are to length of a half wave 
along generatrix 

n = n*(h/R), parameter of number of circumferential 
waves 

So, fifo = amplitudes of shell deflections 

x1 = fi/h, x2 = fo/h, relative amplitudes of deflections 

U,v,w = displacements of points on the middle surface 
along coordinate axes 2,7,z, respectively 

€2",€y",Yzy° == components of relative deformations of middle 
surface of the shell 

Ney = components of variation in curvature and 
torsion of middle surface of the shell 

== stress components in the middle surface 

. = function of stresses in the middle surface: 
oF oF oF 
— = = 0 
ox? ¥? dy? *? Ox Oy 

q = transverse load (pressure), normal to surface of 
the shell and directed toward the center of 
curvature 

Pp = axial normal stress, applied to end faces of the 
shell (compression corresponds to p > 0) 

k = mean tangential normal stress (compression 
corresponds to k > 0) 

p = p/k, parameter characterizing loading of shell, 
ratio of external axial stress to the tangential 
stress 

=— dimensionless parameter of axial 

VEE. h 


loading 


Translated from Design of Three-Dimensional Structures, no. 5, 
State Publishing House of Literature in Architecture, Con- 
struction and Construction Materials (Gosstroyizdat), Moscow, 
1959, pp. 523-536. Translated by Research Information Service, 
New York. 
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Riga 


k R 
V h V h 
parameter of transverse loading 


E,°,E.° = reduced (averaged with respect to thickness) 
moduli of elasticity [see (2 and 5)'] 

@ = modulus of rigidity 

M1", = coefficients of transverse deformation: = 

dD, = E,h3/12(1 — / cylindrical stiffnesses of the 

= Eyh3/12(1 — shell 

Dr = (h3/12, shear stiffness in torsion of the shell 

FE, E>,G,ui,u2 = reduced elastic characteristics of the material of 
the shell in bending and torsion: wk, = 

= 

Q = V 

w = VE\/E, dimensionless elastic 

4 = (4G/V E;E:)(1 — +) parameters of the 
2V material 

B = VEPE2/G — 2V 

6 = 1/12(1 — wipe) 


Thus, the elastic properties of a material inhomogencous 
with respect to the thickness of the shell may be characterized, 
for example, by the independent parameters (5) 


and for an orthotropic shell, homogeneous with respect to 
thickness, by the parameters 


Meaning of the indexes 


B = an index corresponding to the upper value of the 
critical load 

H = an index corresponding to the lower value of the 
critical load 

1,2 = indexes denoting that the elastic characteristics of 


material relate to the x and y axes, respectively 


* * 


We shall consider a closed circular cylindrical orthotropic 
shell of medium length, the material of which has at each 
point three mutually perpendicular planes of elastie sym- 
metry. At any point, one of the planes of elastic symmetry 
is tangent to the middle surface and another perpendicular 
to the axis of the cylinder. 

We shall find critical combinations of loads (axial com- 
pression and external transverse pressure) for loss of stability 
“in the small” and “in the large,’”’ confining ourselves to the. 
case of elastic deformations. 


' Numbers in parentheses indicate References at end of paper: 
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In investigating loss of stability “in the large,” we shall 


consider deflections w as comparable with the thickness of the 
shell, and displacements u and v as small compared with this 
dimension. 

We shall assume that the shell thickness is small compared 
with the radius of curvature of its middle surface. 

As expressions connecting the relative deformations of the 
shell with the displacements u, v and w, we shall take the 
relationships 


Ox 2 Ox? 

Ov ow\? w 

ou , Ov ow Ow 


We shall also use the equation of compatibility of deforma- 
tions (5) 


1 1 1 OF 


ox? “Oy? R ox? 
F We shall consider the shell freely supported at the ends so 
that axially symmetric deformation can exist (sections at the 
supports remain circular). The supports also permit the ends 
of the shell to slide freely in the direction of the z-axis; i.e., 
they do not create thrusts when the shell is bending. The 
axial stress, applied to the end faces of the shell, is then equal 
to the mean axial stress. 
We shall solve the problem by means of Ritz’s energy 
method. We shall therefore put the expression for deflection 
in the form of a series 


[3] 
i=0 


Here ¢(z,y) are functions satisfying the geometrical boundary 
conditions and f; are parameters to be determined. 

We shall express the energy of the systems 9 as a function 
of the parameters f;. Then the equilibrium states of the 
system are determined from the conditions of minimum 
energy of the system 


09/of; = 0 255, [4] 


The solution of the problem reduces to the solution of a 
system of nonlinear algebraic equations, Eqs. 4. 

As in many other treatments (1,3), we shall reckon that for 
loss of stability the deflection of the shell can be put in the 
form 


ny 


= sin — sin — ¢ sin? — sin? — 5 
fot h a R a R [5] 


In this form the deflection satisfies the geometrical bound- 
ary conditions: When xz = 0, aw = fo. Expression 5 does 
not satisfy the condition of equality of bending moments at 
the edges of the shell; i.e., the statical boundary conditions 
are not satisfied; however, this is not required when solving 
by Ritz’s method. 

We shall determine the energy of the system 


a = Ui. + — (Wp + Wy) [6] 
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The potential energy of bending and torsion of the shell is 


O*w 
[> + (Se) + 
(2 y (7] 


Substituting the expression for deflection, Eq. 5, and integrat- 

ing, we get 

= [Dat + + 2Dx)a?B? + DB) fi? + 
+ + + [8] 


where 
a = 
B =n/R 


The potential energy of elongations and displacements of 
the middle surface of the shell is 


dz? dy? GP (2) (9] 


We shall find an expression for the stress function F as the 
integral of Eq. 2, substituting in its right-hand side the ex- 
pression for deflection, Eq. 5 

x y? 

F=F,-—k--— 10 
where F, = F,(z,y) is a particular solution of Eq. 2; the 
coefficients of this function depend on fi, fo, a, B, R and the 
elastic parameters E,°, E2°, and The expression for 
F, is omitted because of its clumsiness. 

Substituting Eq. 10 in Eq. 9 and integrating, we get 


U, =} + + 


R (Agfi2fe + Asfo®) + 


k? 
aRah Bo pk + [11] 


R (Acfi? + Arxfe?)] + 


where A;, As,..., A; are coefficients depending on a, 8 and 
the elastic parameters. 
The work done by the external axial forces is 


where 


oF 
= 


2 
Then 

p -f, E (2 B,° dy? phdxdy 


[13] 


Aaphdy [12] 


Substituting values of w and F from Eqs. 5 and 10 and in- 
tegrating, we get 


w, = athe (72+ p + (2 — pe 
4 41 1 47° Pp Eo” 
[14] 
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4 
= 
= 
rT 


The work done by the normal pressure is 


W,= wadzdy {15] 
0 0 
Substituting w from Eq. 5 and integrating, we get 
W, = 2nRaqfo + rRagqfe [16] 
Making use of the condition that the shell is a closed shell 
=0 [17] 


where 

ov 1 w 

1 OF 1 (2 
\ox? oy? 2 \oy R 

and integrating, we get 


k = — E,°B? + fe + fo) [18] 


We shall introduce expressions for certain derivatives, which 
are necessary for writing out equations of the form of Eqs. 4 


ok 
[19] 
ok 1 
ig fe + 4R [21] 


Taking Eq. 19 into account, the first of Eqs. 4 
09/Ofo = 0 


gives the relationship between the tangential stress and the 
transverse pressure 


k = q(R/h) [22] 


After transforming to dimensionless parameters, the other 
two equations of Eqs. 4 


assume the form 
+ — + ag + — mk* = 0 [23] 


+ — — + (a7 + — 
= 0 [24] 


where 
fo + 1/fw + + + [27] 
[28] 


B+ i/fw 16fw 
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= dyn (s2 +A+t =a) [32] 
l, = 36yn + + [33] 
m= p+1/¢ [34] 


We shall transform the system of equations 23 and 24 in 
the following way: We shall eliminate &* from Eq. 24, then 
solve with respect to x2, and, substituting x,7 in Eq. 23, we 
shall obtain a single equation 


mk* =b + a4 + — [35] 
ay — 
where 
=at+h 
be = a2 — 
bs = 3a, — as 


bs = atl — Fa6+h) 
bs = a3 — fa, 


For an isotropic shell y = Q = w = 1,A = B = 2and6 = 
1/12(1 — The coefficients a,, a3, as, b1, be, bs, bs and 
take on the form obtained by Agamirov [(1), page 328]. 

When axial compression alone is applied (i.e., p = ©), we 
shall have p* on the left side of Eq. 35. This follows from 


* 
mk* = (> + + =) p* — p* when p> 
p pt 

The results of computations using formula 35 may con- 
veniently be presented in the form of graphs of the dimen- 
sionless load parameter k* for different p as a function of the 
geometrical shell parameter H. 

In order to plot the graphs we shall proceed as follows: 
We shall assign certain values to ¢ and for each value of ¢ 
shall compute the coefficients a1, a3, a4, b:, b2, bs, 6s and 6; for 
certain values of 7. Thus, for each pair of values (¢, 7) 
characterizing wave formation at the time of stability loss, 
we shall have definite values of these coefficients. Moreover, 
to each pair (¢,y) for a given s there corresponds a definite 
value of the geometrical shell parameter, since 

[37] 
fn 8? 


Then, assigning certain values to x, and employing formula 
35, we plot the graph of mk* = f(x2) and find the minimum 
value of this function (mk*)min. The values of 2. must be 
assigned so that the value of mk* would be close to the mini- 
mum for any fixed pair (¢,7). 

By way of example we shall present graphs of mk* = f(z2) 
for an isotropic shell for & = V¢ = 0.6 and for values of 7 
from 0.12 to 0.55 (see Fig. 1). 

In computing (mk*) nin We must bear in mind that the re- 
sult is obtained as a small difference between large numbers, 
and to insure an accuracy of up to three to four significant 
figures all the computations must be made with an accuracy 
of not less than six significant figures, and for small values of 
¢(¢ < 0.01) not less than seven to eight significant figures. 

The minima of mk* might be found analytically, but follow- 
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09 0 0 
of Of. 


ing such a course would be very complicated, since it necessi- 
tates solution of algebraic equations of high degree. 
Graphs of k* = f(H) for different p are constructed by first 


(¢,n). On the basis of these computations, we plot a family 
of curves k* = f(H) corresponding to different values of ¢. 
The envelope of these curves also gives the required depend- 


putting s = 1. Assigning the necessary values of p, we find ence of the lower critical stress parameter ky,* on the geo- 
m for the ¢ previously given. Then from the minimum metrical shell parameter ky* = f(H) for a given p. The 
values of mk* we find k* = (mk*)min/m for all the given computations are repeated for other p. 
When only axial compression is acting p* = mk*; ie., 
the value (mk*) min is equal to the lower critical stress param- 
eter p,* for given ({,n). Assigning values s = 1, 2, 3,..., 
and making use of Eq. 37, it is possible to compute depend- 
si encies py* = f(H) for different s. From these we must 
take only those parts of the curves py* = f(A) giving mini- 
a6 mum p,*. This will be proved in the following with reference 
to the design of an isotropic shell. 
In order to determine the upper critical load parameter 
kx*, corresponding to an equilibrium state infinitely close to 
€-06 the initial cylindrical equilibrium state of the shell, it is 
necessary to put 2; = 22 = 0 in Eqs. 23 and 24. 
a5 k* = (l, + ag)/m [38 | 
or 
79085 
0.5 l l 
= —— 106 5 
as 4 + + 
oe aw 1 
y ] [39] 
nto + B+ 1/fa) 
- Hence, it is possible to get dependencies kg* = f(H) for 
0.530 different values of p by a method analogous to that described 
for k,*, which, at the same time, is simpler. 
025 ‘ ; 
as ri ae > For an isotropic shell formula 39 may be written as 
V5 ot 1/? n v6 (& + 1/6)? 
[40] 
We find the minimum value of k* from the equation 
92 
2 3 4 5 z k* 
Fig. 1 on \V6 
id 2 
® | 
1) 
oY] N 
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NS b R | 
| SSS B xpyxckax npocmasnens snavenua 
01 
N 
~ 
4 
| 
= ~ 
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1 t ~ 
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| 
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Fig. 2 *Kg/cm?; ‘values of p are circled 
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We get 


1/8)? Vie + 2p 1 


= 


Substituting Eq. 42 in Hq. 40, we obtain 


V5 ot 1/8] + 2p 1 
[43] 


Assigning different values of &? = ¢ for a certain p. we deter- 
mine kn*/V6 from Eq. 43, 7 V ‘6 from Eq. 42 and the values 
of corresponding to pairs (¢,nV6). We take 
s = 1, and plot in logarithmic coordinates graphs of 


ke* 
V5 


or, What amounts to the same thing 


ketV1 — = f(HV1 — 


This form of presenting the results was chosen in order that 
the graphs might be convenient for computing shells of iso- 
tropic material with any values of F and yu. 

Such graphs are shown in Fig. 2 for p = —10, —4, —2. -1. 
—}, —},0, 3, 1, 2,4, 10 and 100. The value p = 0 corre- 
sponds to transverse external pressure acting alone, p = 3} to 
external pressure acting from all sides on a cylinder closed 
with end faces, p = © to axial compression acting alone, and 
p < 0 to the action of loads of unlike sign. 

When p = | it is possible to get a simple formula for /,* 
(the graph of kg* = f(H) in logarithmic coordinates is a 
straight line). 

From Eq. 42 we find 


whence 


and whens = | 


2 
(1 + 4] 
from Kq. 43 we find 
-(V3+ 
14+ Vv3/ 
or 
9 
= [45] 


V1 — + 1/8) 
and taking Eq. 44 into account 


= 


p?) 3/47 1/2 


when = 0.3, ka* = 0.916/WH. 


Passing from dimensionless parameters to dimensional 


quantities, we get 
h (h 
= 0.916E - | — 46 
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When F = 2 X 10®kg/cm? we have 
ge = 


Papkovich’s formula for critical pressure on all sides. 
When p = @, Eq. 40 yields the known value of the upper 
critical stress in compression 


ps* = 0.6 0r paR Eh = 0.6 


Results of computation of the lower critical stress parameter 
for an isotropic shell p,* are presented in Fig. 3. 

Experiments show that when shells which are not too short 
are subjected to axial compression, a number of waves is 
ordinarily formed along the shell or else the waves are con- 
centrated in a rather narrow band: i.e., for such shells we 
cannot assume that s = 1. 

The minimum value of pyR/Eh = py*, equal to 0.25, is ob- 
tained for € = Vt = 0.58, and for » = 0.15; which, when 
s = 1, corresponds to a geometrical parameter H = 195. 
When s = 2, the minimum of p,* occurs for H = 780; when 
s = 3, for H = 1750, and so on. H = s*(?/€n); see Eq. 
37. 

In Fig. 3 the coordinate H has a logarithmic scale; therefore 
for any s we can get the dependence p* = f(H) by shifting 
the whole curve for p* = f(H1) when s = 1 to the right by 
log s*.. The graph given is analogous to the known graph of 
the coefficient of critical stress of a plate compressed along 
one of its sides. 

For shells with a large parameter H the coefficient py* is 
practically constant with a minimum value of 0.25. 

The value pyR/Eh = 0.25 corresponds to the form of de- 
flection adopted here for loss of stability equation, Eq. 5. In 
other theoretical investigations this quantity ranged from 0.18 
to 0.334. Experimental values of this parameter fluctuate 
within the limits 0.23-0.42. 

Graphs of the parameter of lower critical pressure 


are presented in Fig. 4 for p = 0, 3, 1, 2, 4, 6, 10, 20, 40, 70 
and 100. Note one peculiarity of the variation in the lower 
critical stress parameter as compared with that of the 
upper critical stress. The latter is a monotonic function of 
parameter H for any p (see Fig. 2), and in the computations 


we may confine ourselves to s = 1. The parameter ky* 

varies monotonically only when p < 10, and when p > 10 

it has a minimum for values of H V1 — yu? from 300 to 


170 (for an isotropic shell). 
In this case, as was shown previously in investigating py*, 
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It would be pos- 
sible to reconstruct the ky,* curves for p > 10, putting s = 2, 


it is impossible to take s = 1 for large H. 


Fig. 4. 

In this case, when H > Ho, where Hp corresponds to the 
minimum of the parameter ky*, we recommend assuming that 
ky* is not greater than the minimum value. 

For the sake of comparison, Fig. 4 also shows (dashed line) 
values of kg* for p = 100. As is evident, for p 100 the 
greatest difference between the upper and lower critical 
stresses is obtained when HV1 — bw? = 170. 

We shall establish approximate limits of applicability of 
the graphs computed for an isotropic shell. The upper limit 
of the geometrical parameter H;,, is determined from the 
condition n > 4; i.e., we shall consider that the theory of 
shells of medium length is sufficiently accurate for more 
than four circumferential waves. The lower limit of H is 
defined by the condition of absence of edge effect. 

In the case of loading with an external transverse pressure 
and a comparatively small axial compression (p 1), when 
H is large the upper and lower critical stresses are close to 
each other, and for an isotropic shell we may use the simple 
relationship, Eq. 44, in which, when H is large, we may neglect 
unity compared with 1/&. Since 1/£ = Hn/z?, it follows 
from Kq. 44 that 


, and so on, as was done for pxz*, but this would obscure 


2 24/1] — R? 
When = 0.3 and n = 4 
R? ath 
= (0.22 — 
0.22 ie or < 0.22 


As Timoshenko (4) has shown, close to an absolutely rigid 
supporting contour (this being somewhat different from the 
boundary conditions assumed here), the shell wave length 
along the generatrix for loading with an external transverse 
pressure is equal to 


length of the shell must be at least greater than 4(A/2), i.e., 
a>98VRh. 

Hence Hiim, = 96. 

Thus, in any case, the graph in Fig. 4 may be used in the 
range 96 < H < 0.22 R?/h?. Evidently, it may also be used 
for shorter shells. Thus, the difference between calculated 
and experimental values, obtained by the author (for p = 
0-3), for fairly short shells (smallest value of H was taken 
equal to 35) was 5-9%. 

For loading with predominantly axial compression (p > 
10) at the time of stability loss, a sufficiently long shell has 
several half waves along the generatrix, and, thus, over part 
of each half wave behaves like a shell of medium length. 
However, for a sufficiently large length, the shell may lose 
overall stability, as Euler predicts for a long thin-walled pipe. 

An approximate value of the limiting H\;, may be deter- 
mined from the condition py = pe; the Euler critical stress is 


=C 
From the equation 
h 
R (a 
we get 
iim = 
Pu h? 


For a circular thin-walled pipe ~ QR/V8 = O.7R. 
Since the end sections of the shell do not rotate, c = 4. 
Substituting values of c, imin and py* = 0.25, we get 


\ = 4.9 VRh 
: 2 2 
If we assume that the effect of local bending at the support- Him = 77 - oor h <77 
ing contour is propagated through a distance \/2, then the eR 
a 
* gayent ~ L 
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R 
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Reviewer’s Comment 


In this paper the author develops, by a Ritz energy method, 
‘quations for stability of an orthotropic cylindrical shell 
under combined loading. Eq. 35 of the paper is the stability 
ondition. Starting with this equation, the author confines 
‘is discussion to the reduction of Eq. 35 to a form usable for 
sotropic cylinders. Although these results were interesting, 
the question of whether the stability equation can be re- 
luced any further for easier application to orthotropic shells 
-till remains unanswered. 


The results of this paper could be compared to the work 
of Becker (NACA TN 4237) and the writer (ARS Journat, 
Feb. 1961). In the former, equations are given for critical 
stresses in stiffened (orthotropically) cylinders, and in the 
latter, results are reduced to design curves which allow quick 
determination of critical loads on cylinders composed of an 
orthotropic material. 

—Tuomas E. Hess 
Missile and Space Vehicle Dept. 
General Electric Co. 


Method of Determining the Electrical 


Potential of Bodies in Plasma 


HE USE of rockets and artificial Earth satellites for 

geophysical studies requires measuring their electrical 
potential because the magnitude of the latter influences the 
readings of certain devices located in rockets and artificial 
satellites (1,2).!. In this connection a series of methods of 
determining the potential of bodies sounding the upper at- 
mosphere (1,3,4) has been proposed and employed. 

In (4) the potential of a rocket was determined on the basis 
of the influence of the potential magnitude on the readings 
of a radio frequency mass spectrometer, which determines 
the ion composition of the atmosphere and is located onboard 
the rocket. The effect of the potential shifts the real values of 
the mass numbers of ions vs. the values obtained by means 
of the mass spectrometer. It must be borne in mind, how- 
ever, that this shift can be caused by other factors as well; 
for example, the influence of movement of the sounding body, 
etc. A detailed analysis of the operation of the high fre- 
quency mass spectrometer will not be given [such an analysis 
is given in (2)]; note, however, that the relative error of 
measurement of the electrical potential of the sounding body 
is directly proportional to the magnitude of the error of meas- 
uring of the mass numbers of ions (without taking the clectri- 
cal potential into account), and inversely proportional to the 
magnitude of shift of the mass numbers of ions, dependent on 
the influence of the electrical potential of the sounding body. 


Translated from Artificial Earth Satellites, 1960, no. 4, pp. 
‘61-164. Translated by Astrid Werner, New York. 
' Numbers in parentheses indicate References at end of paper. 
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YA. M. SCHWARTZ 


Inasmuch as the error of measuring can be considerable and 
the shift small, then the specified method permits only an 
evaluation of the potential magnitude. 

It is proposed in (1 and 3) that the potential be determined 
by simultaneous measurements of the electrical field at the sur- 
face of the satellite and the concentration of positive ions 
in the environment. 

This method, as shown in (3), is based on the fact that the 
magnitudes measured by two devices (electrostatic fluxmeter 
and a device for measuring the concentration of positive ions) 
are closely related. In fact, the concentration NV. computed 
from the reading of the device for measuring the concentra- 
tion of positive ions is a function of the satellite potential, 
i.e., Ni =f(V). 

At the same time, the strength of the electrostatic field 
at the satellite surface (in the case where external electrostatic 
fields are absent) is a function of potential and the concen- 
tration of positive ions, Le., = f(N+,V). 

It is clear that this method of determining the potential 
requires the knowledge of the functional relations V. = 
f(V) and E = f(N.,V). The nature of the relation EF = 
f(N+,V) is strongly dependent first on the law of velocity 
distribution of charged plasma particles in the space under 
investigation, and, second, on the angle between the normal 
to that section of the sounding body surface, where the elec- 
trostatic field strength is measured, and the direction of mo- 
tion of the body (3). The complexity, and in a series of cases 
the impossibility, of finding the functional relation F = 
f(N+,V) is the principal disadvantage of this method and 
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greatly limits its applicability. 

The method proposed in (1) of determining the potential 
is based on clear-cut ideas and is simpler to carry out. With- 
out dwelling on a description of the method [it is given in 
detail in (1)], it is pointed out that the accuracy of deter- 
mining the potential is strongly dependent on the extent to 
which the sounding characteristics show the inflection, and 
this, in turn, is dependent on the relation between the ion 
velocities and the velocity of the sounding body and also on 
the nature and quantitative ratios in the mass spectrum of all 
ions. For the purpose of determining the potential by this 
method, it is necessary to know the mass number of heavier 
ions in the space under investigation. 

Another approach to the solution of the problem of deter- 
mining the potential of bodies in plasma, including rockets 
and satellites, is possible. The space charge theory (5) indi- 
cates that the electrostatic field strength at the surface of 
isolated bodies is determined in terms of the equation 


w= (1] 


where 


E., = electrostatic field strength at the surface 


V = potential of the body 
6 = layer thickness of bulk space charge at the surface of 
the body 


k,m = coefficients, which as a function of the gas density, 
may assume different values 


At the same time, the thickness of the layer of positive 
charges at the surface of the body? is determined by the fact 
that the ion current j, passes through this layer for a potential 
difference V, and the initial ion energy is low compared with 
V. The relationship between j,, V and 6, in agreement with 
(5), can be expressed in terms of the equation 


Jn = A(V"/6*) [2] 


It is clear that if the numerical values of coefficients A, m 
and k& are known, then, measuring the electrostatic field at 
the surface of isolated bodies (including satellites or rockets) 
in plasma and the density of the ion current at the point where 
the device for measuring the electrostatic field strength is 
located, it is possible to determine the electrical potential V 
of the body and thickness 6 of the space charge layer at the 
point where the device for measuring the electrostatic field 
strength is located. For example, the equation for the 
potential of the body can be expressed in the terms 


Experimental procedures for the proposed method are 
simple; for example, the method of determining the magni- 
tude of the ion current at the surface has been developed 
sufficiently already. The method of measuring electrostatic 
fields at the surface of bodies in plasma has been developed to 
a lesser extent; however, as shown in (3), measurements of 
this kind are simple in principle. It is emphasized that the 
present method, based on utilizing the boundary conditions 
at the surfaces of bodies in plasma, is largely free of errors 
which are brought about by the requirement of exact know!- 


* Let us note that the body placed in plasma, as a rule, is 
charged negatively as verified by experiments described in 
(4, 6 and 7) for rockets and satellites. 


edge of the parameters of the plasma (concentration, tem- 
perature) and the law of distribution of particles in accordance 
with velocity. 

Concerning the determination of the numerical values of 
coefficients k, m and A, it is quite clear that the calculation 
of coefficients k and m is simple. Under the conditions of 
ionospheric measurements by rockets and satellites, coeffi- 
cients k and m equal 2 and 1.5, respectively. The value of 
coefficient A is dependent on the ion composition of the plasma 
under investigation, on the ratio of layer thickness 6 to the 
radius of the electrode probe for measuring the ion current 
(5). However, in the case of a plane electrode probe, th 
value of coefficient A is not dependent on the geometry of 
the layer. For a plasma produced in laboratory conditions, 
determination of A is simple because the composition of such 
a plasma is known beforehand. In the case of the ionosphere, 
however, matters are more complicated, but the circumstanc« 
that atomic ions of oxygen and nitrogen, which are the prin- 
cipal components of the ionosphere (6, 7) at altitudes above 
230 km, have close mass numbers (16 and 14, respectively 
considerably simplifies the problem. 

If it is taken into account that A ~/1/M (M is the mass 
number of the ion) then the error which results when onc 
of the given molecular weights is selected as the base for the 
calculations is small. At altitudes less than 250-230 km 
the composition of the ionosphere is infinitely richer in var- 
ious kinds of ions. The accuracy of determining the numer- 
ical value of coefficient A at these altitudes is dependent on 
our knowledge of the composition of the atmosphere at 
these altitudes. Much of the data available at the present 
time is at variance. The presence of negative ions at these 
altitudes may also complicate the interpretation of the  re- 
sults of measurements of the electrical potential by the spe- 
cified method. It is clear that every individual case requires 
a detailed analysis of the conditions of measuring. The 
condition which, in principle, limits the application of the 
method of measuring the potential is that Eqs. 1 and 2 be 
satisfied at that portion of the body surface where the devices 
for measuring electrostatic field strength and current of 
positive ions are located. For a fixed sounding body, these 
conditions are observed over the complete surface of the body. 
In the case of bodies moving at a high velocity, satellites for 
example, the applicability of Eqs. 1 and 2 is limited. They 
are not satisfied on the rear portion of the surface of the 
moving body where an area with electrons primarily is formed 
(1). Therefore, for determining the potential, it is neces- 
sary to use only those data obtained by the device for measur- 
ing the electrostatic field strength and the device for measur- 
ing the current density of positive ions which refer to instances 
when the normal to their surface is not at angles greater than 
90 deg in relation to the course of the moving body. 

The specified method of determining the electrical potentia! 
of satellites and rockets can be used for determining the po- 
tential of any body placed in plasma. 
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Reviewer’s Comment 


In this article, the point is well made that the absolute 
determination of vehicle potential from radio frequency mass 
spectrometer data is a rather difficult process. If corrections 
for the effects of vehicle velocity and aspect upon the posi- 
tion and amplitude of the mass peaks are made, however, 
upproximate values for the vehicle potential may be readily 
determined. With care, values of the potential accurate to 


1 


within + 3 v may be obtained, from measurements of mass 


position. 


Although the measurement of the vehicle potential may 
‘¢ approached thus using the RF spectrometer, this instru- 
ment would, of course, never be chosen as the prime sensor for 
such data. Obviously a less complex device which is funda- 
uentally sensitive to vehicle potential would be the logical 
hoice for a potential sensor. The proposed method which 


requires the measurement of both electrostatic field and ion 
liffusion current also appears to be rather complex and 
.dmittedly of limited application. 

The author suggests the measurements of field and current 


should be made in close proximity along the vehicle skin. 
This in part infers the actual physical case, namely, that in 
many instances the vehicle skin is not an equipotential sur- 
face. This is due not only to the fact that the electron sheath 
is distorted because of vehicle velocity, but also to the skin 
potentials induced from traversing the ambient magnetic 
fields. In the satellite case, this situation has been investi- 
gated and good vehicle potential data have been obtained 
from a comparatively simple, modified Langmuir probe ex- 
periment.!_ From such studies it is quite apparent that the 
complex nature of vehicle electrification requires an integrated 
study of the several related parts that go to make up the 
whole, rather than a single measurement at single point, if 
the composite model of the vehicle plasma sheath is to be 
understood. 
—H. TayLor 
Goddard Space Flight Center 
NASA 


1R. E. Boureau et al., “Measurements of Sheath Currents 
and Equilibrium Potential on the Explorer VIII Satellite,’ pre- 
sented at American Astronautical Symposium, March 17, 1961. 


Correlation Method for Analysis of 
Large Heterogeneities 


of the lonosphere 


Lhe FACULTY of propagation, generation and transmis- 
sion of radio waves of the Physics Department of 
Moscow State University is doing a great deal of work on 
analysis of large (dimensions of the order of tens and hundreds 
of kilometers) heterogeneities of the ionosphere. This 
analysis is being done by the phase shift method. It has 
been found that correlation analysis is the most suitable 
method because of large variability of the ionosphere and, 
correspondingly, large variability of the phase of received 
waves. Although the basic ideas of this analysis, founded on 
veometrical interpretation, have been published (1,2),! 
exceptional complexity and insufficient evidence in separate 
cases (especially during statistical anisotropy of heterogene- 
ities) has made it necessary to work out a detailed analytical 
ipparatus for the correlation method of treatment of experi- 
mental data. 


Translated from Bulletin of Moscow Universiiy, 1959, no. 6, pp. 
‘7-98. Translated by Ernst Volgenau, Los Angeles, Calif. 
1 Numbers in parentheses indicate References at end of paper. 
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In analyzing large heterogeneities of the ionosphere the 
measured parameter is the phase of the wave ¢(z, y, t) for 
which the correlation function p(& 7, 7) is defined by the 
relationship 


nea yb — OF 
[y(z, y, OP — y, OF 


7, 7) = 


where (x, y) and (2 + & y + n) are the geometrical coordi- 
nates of placement of the measuring apparatus, and tand ¢ + 7 


are temporal moments of comparison of measured phases. 
The bar indicates the statistical mean. 


From the definition of the correlation function, it follows 


that the function is limited to |p! < 1 and is even, that is, 
p(é, n, T) = (—& — n, — 7). Considering £, 7, 7 as rec- 
tangular coordinates, it is possible to obtain all of the prop- 
erties of the correlation function from the analysis of the form 
of the surface p(&, 7, 7) = constant. It is known that if two 
symmetrical points with coordinates (£ + & mo + 7, To + 7) 
and ( — & mo — 7, To — 7) lie on one surface, such a surface 


valled a central surface, and the point (£0, 70, 70) is its center 
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(3). From the even nature of the correlation function, it 
follows that the surface of the correlation function is a central 
surface, and the point (0, 0, 0) is the center of this surface. 
Hence it is possible to represent the correlation function in the 
form p = f(u), where f is an arbitrary function satisfying the 
condition |f(u)| < 1 and u is polynomial of even degree in 
£, 7, r. In order to obtain the fundamental characteristics 
of the geometry and the time changes of the correlation func- 
tion, a sufficient approximation of its argument is the quad- 
ratic form in £, n, 7. During this calculation, higher degrees 
obviously give more refined properties of the correlated 
function. 
Thus, let us assume that 


flan? + a22n? + + 

+ 2aisét + [1] 
where ax, ~ 0 (i, k = 1, 2, 3). Since the coordinate axes 
being utilized are oriented arbitrarily relative to the principal 
axes of the given quadratic form, it is possible to estimate the 
shape of this central surface of the second order by the nature 
of the intersections with three mutually perpendicular planes 
— = constant, 7 = constant and 7 = constant. Obviously 
the form of these sections is defined by the sign of the cor- 
responding invariants 6,, 6,, 5¢ 


6, = — [2] 


etc., and the following relationships occur 


= FOF On? 3] 


etc. But since (9) 


Of 
1 
On? t=n=r=0 Or 
| o* \Ox Oy 


etc., for the definition of the signs of the invariants it is pos- 
sible to utilize the inequality of Boonyakovsky (4) 


< 
| 
Then on the basis of relationships 3 and 4 we obtain 
6, 2 0 6,20 &20 


But because of the even nature of the correlation function it is 
necessary to exclude noncentral surfaces from consideration; 
finally we obtain 


5 >0 6,>0 &>0 


That is, all sections are ellipses. If three arbitrary mutually 
perpendicular sections of the central surface of the second 
order are ellipses, then that surface is an ellipsoid. Let us 
pause here to discuss further the coefficients of the positively 
definite quadratic form obtained previously. During re- 
flection of radio waves from the ionosphere, the time varia- 
tion of the wave field and, in particular, of its phase, are deter- 
mined in the first place by the statistical variability of the 
ionosphere, and, secondly, by drift in the ionosphere, that is, 
by the translational motion of heterogeneity. 

During this drift the time variation of the phase is greater 
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the greater the drift velocity of the heterogeneity. It follows 
that the smallest variation of the phase with time will be in 
the system of coordinates connected with the moving hetero- 
geneities of the ionosphere. Utilizing the transformation of 
Galileo (5) 

t=2' + oat y=y' + [5] 
where x, y are the coordinates of the fixed system, and 2’, y’ 
are the coordinates of the system moving with the speed 
Vo(Ura, Yya). Taking into account the fact that the time 
variation of the phase of the waves is characterized by the 
quantity (dg/dt)?, we will obtain for the moving system 


(te + Uraty y! + Yyat 
dt 


) 2 


+ Zayv vy + Zai3vz + + as3)f’(0) 


According to the preceding, v,(v;, vy) must make (d¢/dt): 
a minimum. Thus for the determination of vz, and vy, we 
obtain two equations 


Hence 


— 


— 


Ura 

From Eqs. 5 and 6 it is easy to obtain the correlation 
function, Eq. 1, of the phase of the field in the system of 
coordinates fixed with respect to the ionosphere 


po = flan’ + + a'33T?) 


In addition, if the coordinate system is rotated in the &, 7 
plane so as to eliminate the coefficient associated with the 
£’»’ term in the new coordinates, the correlation function 
acquires its simplest form 


fla'né” = a'son’? + a’33r?) 


Where a’, a'22, a’33 are new coefficients. We shall designate 
= 1/a?, = 1/b?, a’33 = 1/7,2, where a and 6 are geo- 
metrical characteristics, and 7, is the time characteristic of 
the correlation function. These characteristics define ‘“di- 
mensions,” or radii of the correlation function in space and in 
time. It may be pointed out, however, that a and 6 charac- 
terize the variation of the phase of the wave in space, and 7, 
characterizes the variation of the phase of the wave in time. 
Thus 7, may also be called time of relaxation or time of spread- 
ing. In this manner, the correlation function of the phase 
of the field in the system of coordinates fixed with respect 
to the ionosphere is 
12 12 2 
T 
f b2 
Expression 7 easily enables one to use corresponding ex- 
pressions in the system of coordinates fixed with respect to 
Earth. For that, it is necessary to utilize relationships 5 and 
6, expressed in polar coordinates 
=rcosa n = 7 sina 
= Ura COS @ — Veg SIN 
= Urd SIN @ + COS @ 


where r is the magnitude of the base separation between 
the two measuring points in the direction a from the major 
axis of the ellipse a; v,d0, Uyao are the projections of the velocity 
vector of drift on the major axis of the ellipse in the plane 
£, 13 Urd, Yaa are the radial and transversal components of the 
velocity of drift relative to the base separation. After simple 
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transformations, Kq. 7 is given by 
p(r, a, 7) = f(Ar? — 2Drr + cr?) [8] 


COS @ sin @ 
Vrdo + 2 Vydo = Abra Fvea 
a b 


a b Te Te? 


1 cos? a sin? a 
a? + b2 


1 l 
F = sin @ cos @ 


COs? Qa sin?aq 
a? 6? 


Ay 


vhere aq is the direction of the drift. 

Thus one must determine six unknown characteristics of 
1c correlation function: a, b, ao, Ua, Qa, Te, Where ap is the 
lirection of the major axis of the ellipse in the &, 7 plane 
relative to the north-south direction. The indicated charac- 
‘eristics are defined from the set of functions p(r, a, 7), in 

which r and @ are parameters, and 7 is a variable quantity. 

When the measuring equipment is located in several fixed 

points on Earth, the above formulation is characteristic of 
-tudies of statistical properties of phase of scattered radio 


waves, 


NI 


The geometric characteristics a, 6, ao are determined from 
the comparison of the space and time correlation functions 


p(r, = p(T.) 
or 


Hence 
= tVC/A (9] 


According to Eq. 9, the end of the vector v’, describes an 
ellipse, called the characteristic ellipse, and v’, is the charac- 
teristic velocity. Since the characteristic ellipse is a central 
one, then for its entire definition it is necessary and sufficient 
to define v’, in three non co-linear base directions. Since the 
orientation of the base directions with respect to north-south 
is known, @ can be determined from three values of v’,. 
However, the formula for its definition in the general case 
turns out to be quite complicated and of little use for practical 
calculation. Ordinarily, therefore, ao is defined by geo- 
metrical construction according to the characteristic ellipse 
obtained. It also follows from Eq. 9 that 


VC VC 
the coefficient C may be obtained from the definition of the 
time radius of correlation 
Cro? = 1 [11] 
which satisfies the equation 
p(ro) = f(l) = ¥ 


ind 0 < y < 1 is defined by the form of the function f(u). 
Having substituted Eq. 11 in Eq. 10 it is possible to obtain 
the final expression for the dimensions of the semi-axes of 
the space characteristic ellipse 

max T0 b= v’. min 70 


where v’, max and v’, min, for the same reason as Qo, are ob- 
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tained graphically from the ellipse of characteristic velocities 

During the study of the space characteristics of the cor- 
relation function, the important parameter is the parameter 
of anisotropy e which may also be defined from the graphic 
constructions of the ellipse of characteristic velocities 

@ _ We max 


e= = 


b min 


§2 


The determination of the characteristics of drift vg and a 
is carried out with the help of the introduction of the so- 
‘called apparent velocity v’. Since r and @ are parameters 
and 7 is a variable, let us consider the intersection of the cor- 
relation function p(r, a, 7) with the straight line p = con- 
stant parallel to axis 7. By the use of equation p(r, a, 7) = 
constant, we obtain the equation 


Cr? — 2Drr + Ar? = const = B [12] 


from which it is possible to determine the coordinates of the 
point of intersection 71, 72 of the straight line with the cor- 
relation function 


Dr + V D*r? — ACr? + cB [13] 
1,2 = C 


From Eq. 13 it follows that the points 7; and 72 are located 
symmetrically with respect to the straight line 


r=r' = (D/C)r [14] 


passing through the maximum of the correlation function. 
Indeed, the apparent velocity v’, defined as v’ = r/r’, is 
connected with the maximum of the correlation function. 
As a result of the substitution of v’ in Eq. 14, and by the use of 
Eq. 8, the equation for the determination of v’ is obtained 


COS @ sin @ 
( a + be =C [15] 
This may be expressed as 
v’ = p/cos (a — yp) 
where 
By 
Uydo 
tan y = — 
Urdo 


Eq. 16 is the equation of a straight line in polar coordinates. 
This straight line is called a front of apparent velocities. It is 
evident that on the front of apparent velocities lie all the 
ends of the radii vectores of apparent velocity v’. From 
Eq. 15, with the utilization of Eq. 9, it is easy to obtain 
(v'/v’.)? = AC/D*. But Eq. 12, according to the foregoing, 
is the equation of an ellipse; hence it is necessary that 6 = 
AC — D?>0Oandv’? > v’.? for all values of a. It follows that 
the front of apparent velocities does not intersect the ellipse 
of characteristic velocities. From Eq. 16 it follows that the 
front of characteristic velocities forms with the polar axis 
(the major axis of the characteristic ellipse) the angle a; = 
y + 7/2; but then, Eq. 17 transforms to 


tan ag tan ay = —b?/a? [18] 
because 
Vydo/Urdo = tan ag 
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Consequently, on the basis of Eq. 18 the direction of drift 
and the direction of the apparent velocities are conjugate 
diameters of the ellipse of characteristic velocities. In 
addition to this, the end of the drift vector points in the direc- 
tion of the front of the apparent velocities. Eq. 18 indicates 
also the obvious method of geometrical construction for the 
determination of the direction of the drift of large hetero- 
geneities. For the determination of the magnitude of the 
drift velocity, we will write Eq. 15 utilizing Eq. 8 in the form 


v'(Avra + Fvaa) = C 
It is clear that in the direction of drift veg = 0; hence 
= 


or 


C cd 
— = {19] 
A 
where v’.q and v’¢ are characteristic and apparent velocities in 
the direction of drift. Since F = 0 along the principal 
axes of the characteristic ellipse 


U min tan @ lv a 


where v’, and v’, are apparent velocities along the corre- 
sponding principal axes of the characteristic ellipse. 


$3 
After the determination of the magnitude and direction 
of drift, by utilizing Eqs. 8 and 11 it is easy to find the time 
of relaxation 7, of the heterogeneities 


l 


or 


To [20] 


V1 V1 va 


Together with the time of relaxation 7,, it is convenient to 
introduce the so-called velocity of relaxation v, defined by 
the relationship 


1/7.VA 


which, based on Eqs. 9, 11 and 20 is equal to 


; 

Ud Vd 
Consequently, the characteristic ellipse of relaxation velocity 
is similar to the ellipse of characteristic velocities v’, and is 
distinguished from the latter only by dimensions. Besides 
the absolute characteristic of the time change of the hetero- 
geneities 7,, as seen from Eq. 21, it is convenient to utilize 
the relative characteristic of that variation vy = v'.q@/v'q 
changing within the limits 0 < <1. 


dg(z, y, t) 


= lim 


dt dt to—>0 
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The method of determining all the basic parameters of the 

heterogeneities of the ionosphere was given in the foregoing 
on the basis of the analysis of the various properties of the 
correlation functions. However, one obvious relationship 
still remains. Here we are concerned with the investigation 
of the equation 

er) | P(Ter) 
corresponding to the equation for the argument of the cor- 
relation function 

Ar? — 2Drr., = CTer [22] 
Designating by v,, the so-called crossing velocity v., = 1/27-, 
and utilizing Eq. 8, relationship 22 may be rewritten in the 
form 
SIN — Vydo)® _ 


(Yc, COS @ — Urdo)? 
b? 


y2 2 

a’ Te 
Eqs. 22 and 23 indicate that the end of the radius vector 
v,, describes an ellipse with semi-axes smaller than those of the 
ellipse of the characteristic velocities, and displaced relative 
to the latter in the direction of drift by the amount 3vg. 
That is, the ellipse of crossing velocities is not a central ellipse. 
Since the center of the ellipse of crossing velocities is in this 
case an unknown quantity, its complete definition requires 
four points (aside from the origin of the coordinates). In the 
case of the ellipse of the characteristic velocities, three points 
are sufficient. Accordingly, the analysis of the ellipse of 
crossing velocities is not carried out in practice, since it 
requires at least four base directions, even though the utiliza- 
tion of this analysis would be convenient for the definition 
of the magnitude and direction of drift. However, one of 
the formulas following from Eq. 22 is of definite inter- 
est, from the viewpoint of the verification of the entire theory. 
Specifically, on the basis of formula 8, the following relation- 
ship holds 


This may be rewritten in the form 


= [24] 
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Other methods are also used for the investigation of large 
heterogeneities in the ionosphere: 

a Measurement of the time variations of the derivative of 
phase with respect to time dg/dt [see (6) ]. 

b Measurement of the time variations of angles of arrival 
0¢/dx, 0¢/0, of the radio waves scattered by large hetero- 
gencities (7). 

By virtue of the extreme complexity of such measurements, 
ordinarily only one of the methods of recording large hetero- 
geneities is utilized, i.e., ¢ or dg/dt, or and 0¢/Oy. 
A comparative analysis of the measured data is possible only 
when the mutual dependence of the given data is correctly 
understood. 

The following section is devoted to this mutual dependence 
of various parameters obtained during the correlation proc- 
essing of time variations in the space-scattered points ¢, 
dg/dt or 0¢/Oy. 

Let us start with the expression 


t+ — o@+ hy + at t+ Di le@yt+b) — _ 


to? 


g(a + &y +7, t+ r)e(z,y,t) [25] 
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Therefore the correlation function of the velocity of the 
phase with respect to time p,(&,7,7) will be equal to 


0*p(E, 7, 7) ==7T=0) 
= — - = ~~ Yt 
p-(&, 7, 7) [26] 
Similarly 
7, 7) = = 7 = 0) 
PHS, 0, = Oe? of 
Py\9, 7, on? on? 


\Vhen using polar coordinates, it is possible in place of 
o:(&, 9, 7) and p,(&, 9, 7) to consider p,(r, a, 7) and p,(r, 
« + 2/2, 7). In connection with this, only the expressions 
o-(r, a, 7) and p,(r, a, 7) will be analyzed further. First of 
oll, it is necessary to point out a very interesting result con- 
erning the possibility of determining the dispersion of angles 
of arrival during the recording of the time variations of the 
hase of the reflected wave in the space-scattered points. 
rom Eqs. 25 and 8 it follows that 


(% =o = 0) _ 
Or? 


2 2 ( 
(5°) = — Af"(0) 
or? 


—o?Cf'(0) 


From this 


or \ dt 


According to the general considerations stated, it is possible 
to represent the function p, in the form 


fi( Ayr? 2Dirr Cyr)? [29 


where the parameters A,, Di, C; may be distinguished from 
the analogous parameters A, D, C for the function p(r, a, 7). 
Also, the apparent velocity, which is equal to v’, = C,D,, 
may be introduced for p,(r, a, 7) as well as for p(r, a, 7) 
Thus 7’ = r/v’,, as in Eq. 14, corresponds to the time maxi- 
mum of p-(r, a, 7), that is 


p-(7, a, 7’) = 0 [30] 
But then from Eqs. 8, 26 and 30 it follows that 
+ f’"(u) 0 [31] 


where u = Ar? — 2Drr + Cr*. Since Eq. 31 must hold for 
the arbitrary form of the function f(w), the only solution that 
remains is 


du 
= Cr’ — Dr) = 0 

or 
a 


Consequently, the apparent velocities of the arbitrarily 
oriented bases for the functions p(r, a, 7) and p,(r, a, 7) 
coincide. Having utilized the functional representation of 
iq. 8, the functions p(r, a, 7) and p,(r, a, 7), Eq. 832 may be 
represented in the form 


C C 
D, = Aira + Fitaa = =A Va + F Vad 
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Since the magnitude and direction of the drift may have any 
arbitrary values, then 


Cie. 
A, = A F,= ae 
But from here it also follows that 
= = pt? 
cr A; A c 


That is, the ellipses of characteristic velocities p,(r, a, 7) 
and p(r, a, 7) completely coincide, and v., = v, and €, = e€. 
The parameters will be different provided they are determined 
by the time radii of correlation zo, and 7». Since to; < 7», 
then a, < aand 6b, <b. The basic properties of the function 
p,(r, a, tT) may be obtained from the investigation of the 
crossing velocity t,,,,. If one considers p(r, a, 7) as the 
function of rand @ with 7 as a parameter, then the maximum 
of this function is attained for r’, which is defined by the equa- 
tion 


(0/Or’) p(r’, a, 7) = 2(Ar’ — 2Dr)f'(u) = 0 
From here 
= D/A 


The analogous position must also exist for p,(r, a, 7); that is, 
= But, according to Eq. 27, v,,,,can be defined 
from the equation analogous to Eq. 31 


0% Ou \* OU Ou 


or correspondingly from 
Ou/oOr = 2(Ar’ — Dr) = 0 
From here 


dD, D 

for an arbitrary direction of the bases; that is, the ellipses of 
crossing velocities for p,(r, a, 7) and p(r, a, 7) completely 
coincide. By virtue of Section 4 (in which ellipses of crossing 
and characteristic velocities are similar and_ identically 
oriented), the magnitude and direction (1, = J) of the ani- 
sotropy coincide. Moreover, from Eq. 27 it follows directly 
that a, > a,b, > b and 7,, < 7. Information concerning 
quantities 7,,,, v’,,, and v’, and their connections with 7,, 
v’. and v’, cannot be obtained within the limits of the theory 
being described. However, it is possible to get some informa- 
tion from general considerations of the pertinent processes 
of spreading or relaxation. First of all, on the basis of Eq. 8 


v2 + + a\? 
= — = —=h+k 
Te 


similarly 
=h+ k(a,/te, +)? 


If one is to make the natural assumption that the spreading 
of heterogeneities of the ionosphere is governed by the dif- 
fusion equation (8), then the initial dimensions of a and 7, 
are related by a/WV7, = constant. [In (8) the case of the 
spherical heterogeneities is considered. In the case of ellip- 
soidal heterogeneities, it is possible to show that the solution 
has a similar form.] From this relationship when a, < a, 
Tey < Te, it follows that (a,/7.,,)? > (a/r-)?. It also follows 
that v’,? > v’?, and on the basis of Eq. 24, v’., > v’c. 

In Sections 1 to 3 methods are given for determining various 
parameters of the heterogeneities based essentially on the 
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assumption of concentricity of the ellipses obtained in the 
cross section of the correlation function. It is evident that 
this assumption may not hold in all cases. Consequently 
the space characteristic ellipses in the various sections may 
have a different orientation and different parameters. In 
particular, an asymmetrical run of the function p(r, a, 7) 
relative to 7 may develop. In these cases 7, v’., v’ for the 
different bases (determined according to methods stated in 
Sections 1 and 2) will correspond, in the general case, to 
various cross sections of the correlation function, and there- 
fore are not suitable for the comparison and determination of 
all the parameters sought. In such cases v’, = C/A and 
v’ = C/D should be determined from a single standard cor- 
responding to the standard used in determining time radius 
of correlation 7>. Comparison of the cross correlation func- 
tions p(r, a, 7) with the autocorrelation p(7) on the indicated 
standard gives a system of two equations for the definition 
of v’ and v’, 


p(ri, a1, 71) = p(r)o + Ar? — 2Drr, + Cr? = Cro? 


p(n, Qi, 71) = p(n, a, + A, 2Drr, + Cr,? = 
Ar? — 2Drr2? 


where 7; and 7, are the abscissas of the points of intersection 
of the horizontal straight line with the correlation functions 
by the standard of definition of 7>. From Eqs. 33 it follows 
that 


A To" + TiT2 D 


=——_ [34 
T1 + T2 


Utilizing Eqs. 34 it is possible to carry out the determination 
of a, b, va, Qo, Tc, @a by the methods described in Sections 1 to 3. 
In connection with this, the only, but very sensitive, criterion 
of legitimacy of utilization of the correlation function in the 
form p = f(Ar? — 2Drr + Cr?) is the rectilinearity of the 
front of apparent velocity. That is, in the presence of a 
three-base measuring system the ends of all three apparent 


Reviewer’s Comment 


The basic concepts developed in this paper appear to be 
the standard ones for the analysis of spaced-receiver fading 
data. The mathematical approach is different, however, 
is less geometrical and makes less appeal to the physical 
considerations involved than is typical of earlier develop- 
ments. The relaxation, or diffusion, time of the irregulari- 
ties is discussed. It is perhaps unfortunate that no observa- 
tional details were given to illustrate the method of treat- 
ment. 

Other workers also have commented on the large amount 
of work involved in the full correlation analysis, but the avail- 
ability of fast computational facilities can largely neutralize 
this feature. 

The more recent papers in this field have tended to be 
mainly experimental in nature (1-5). Kent (6) and Frihagen 
and Tréim (7) have taken advantage of the special properties 
of signals from artificial satellites in obtaining height data on 
irregularities. Mitra et al. (8) have presented data on drifts 
during a magnetic disturbance. 

There still appears to be a need for more long-term sta- 
tistics on drifts; the data on the geometry of the irregulari- 
ties appear to be conflicting, and there are few data on the 


1328 


velocities must lie on one straight line. 


Conclusions 


A large amount of work has been performed on the proc- 
essing of the material from observations with respect to the 
time variations of phase of the reflected wave in three space- 
scattered points. The preliminary results of determination 
of various parameters of large heterogeneities are discussed in 
(6). Analysis of the experimental material has shown that, 
within the limits of accuracy of determination of various 
parameters, theoretical considerations used as the basis of 
this article hold quite well in the majority of cases. Com- 
parison of large heterogeneity parameters obtained by means 
of the phase and the angle method (7) was carried out in (6) 
on the basis of the analysis in Section 5. Inthe comparison, 
the good coincidence of these parameters was pointed out. 

In conclusion, it may be mentioned that all the results ot 
the present article may be applied for the correlation analysis 
in the investigation of small heterogeneities by the amplitude 
method, and in the study of average heterogeneities by the 
radio-astronomical method when angles of ariival of the radio 
waves are near the vertical. 
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distribution of heights in a particular experiment. Then 

there are always the final questions of the origin of the 
irregularities and their lifetimes. 

—A. W. 

Defence Research Telecommunications 

Establishment 

Ontario, Canada 


References 


1 Reid, G. C., **The Variation with Siderial Time of Radio Star Scintilla- 
tion Rates,”’ Canadian J. Phys., 1957, vol. 35, p. 1004. 

2 Chivers, H. J. A., “The Simultaneous Observation of Radio Star 
Seintillations on Different Radio Frequencies,"’ J. Atm. Terr. Phys., 1960, 
vol. 17, p. 181. 

3 Jones, I. L., ‘Further Observations of Radio Stellar Scintillation,” 
J. Atm. Terr. Phys., 1960, vol. 19, p. 26. 

4 Chivers, H. J. A., ‘Observed Variations in the Amplitude Scintilla- 
tions of the Cassiopeia Radio Source,” J. Atm. Terr. Phys., 1960, vol. 19, p. 
54. 

5 Gruber, §., ‘Statistical Analysis of Radio Star Scintillation,"”’ J. Atm. 
Terr. Phys., 1961, vol. 20, p. 59. 

6 Kent, G. 8., “High Frequency Fading Observed on the 40 Mc. Wave 
Radiated from Satellite 1957a, J. Atm. Terr. Phys., 1959, vol. 16, p. 10. 

7 Frihagen, J. and Tréim, J., ‘‘Scintillations of Signals Transmitted from 
Earth Satellites,” Norwegian Defence Research Establishment, Oct. 1960, 
Tech. Note no. 1, Contract AF61(052)-344. 

8 Mitra, S. N. et al., ‘Horizontal Drift in the Ionosphere over Delhi,” 
J. Atm. Terr. Phys., 1960, vol. 19, p. 172. 


ARS JouRNAL SUPPLEMENT 


7 
~ 


» the 
yace- 
ition 
ad in 
that, 
rious 
is of 
‘om- 
eans 
(6) 
Ison, 
ts ot 
lysis 
tude 


the 
‘adio 


1953 


» 481 
izdat 


Mos- 
SCOW 
\iyar- 

Nauk 


Soviet 


‘hen 
the 


DEY 
ions 
nent 
radia 


tilla- 


Star 
1960, 


ion,”” 


tilla- 
19, p. 


Atm. 
Vave 
0. 
from 


1960, 


ENT 


Magnetometric Equipment of the 
Third Soviet Artificial Earth Satellite 


EASUREMENTS of Earth’s magnetic field in the upper 

atmosphere, as envisaged in the International Geophysi- 

il Year program, were accomplished in May—June 1958 on 
‘he third Soviet artificial Earth satellite (1). 

For the first time, detailed information was obtained on the 
-patial distribution of Earth’s magnetic field in the altitude 
ange of 230-1880 km. Although the experimental data cov- 
red mainly the Eurasian region, south of the 65th parallel, 
ihey nevertheless permit the solution of some general ques- 
‘ions of geomagnetism, since the magnetic field within the ter- 
ritory of the Soviet Union contains all the peculiar features of 
-trueture of Earth’s magnetic field—global anomaly, regional 
ind local anomalies. 

The period during which the magnetic measurements were 
uade (May 15-June 5) was marked both by disturbed and 
calm states of the magnetic field and hence was suitable for 
the investigation of questions relating to the variable part of 
-arth’s magnetic field. 

The experimental results obtained were processed and ana- 
lyzed. The results will be published as the work is completed. 

The present paper describes the construction of the mag- 
netometer, questions relating to its mounting in the satellite, 
information characterizing the stability and precision of the 
apparatus, the nature of the experimental data obtained, and 
the method of allowing for the effect of magnetic deviation. 


|. General Apparatus Requirements 


The principles of measurement, the arrangement, and met- 
rological characteristics of the magnetometer were dictated 
by the geophysical aims of the experiment and by the condi- 
tions of measurement on the satellite. 

Magnetic field measurements in the upper atmosphere are 
being carried out in order to obtain experimental data on the 
spatial distribution of Earth’s magnetic field intensity, to de- 
termine its homogeneity, and the nature and absolute values 
of gradients. 

Another problem is to discover the nature of the changes in 
the terrestrial magnetic field at high altitudes during intense 
magneto-ionospherie disturbances. 

The parameters of the magnetometer (sensitivity, meas- 
urement range, speed of response) were fixed with reference to 
the main aims of the experiment (2). For evaluation of the 
parameters, the values of the normal field gradients, known 
from calculations and terrestrial maps, could be used. As we 
know, the vertical gradient of the total vector of the dipole 
part of the terrestrial field is given by the formula A7T/AR 
= 37/R, where R is the distance from Earth’s center. The 
latitudinal gradient of the normal field at Earth’s surface is 
~3.5 y per km (1 y = 107 oc.) For the fixed orbit we 
should expect that the minimum field values (over the equator 
region) would reach 15,000 y. From geophysical and practi- 
cal considerations, the maximum field value measurable by 
the magnetometer was chosen as 60,000 . 


Translated from Artificial Earth Satellites, 1960, no. 4, pp. 
135-160. Translated by Research Information Service, New 
York. 

' Numbers in parentheses indicate References at end of paper. 
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Some uncertainty arose in the evaluation of the possible ef- 
fects of fields of ionospheric current systems, since an analysis 
of the surface variations of the magnetic field does not give a 
unique value for the intensity of equivalent current systems. 
On Earth’s surface the amplitude of daily variations on quiet 
days (S,) attains an average value of 40-50 y, i.e., approxi- 
mately 0.1 % of the value of Earth’s constant field. Current 
systems, the thickness of which is small in comparison with 
their transverse dimensions (low latitude current systems), 
resemble a thin homogeneous current layer. 

For this reason the field intensity close to such a layer should 
not differ appreciably from that observed at Earth’s surface. 
However, magnetometer measurements at altitudes of 90-105 
km revealed current system effects around 400 y (3). 

The magnitude of the disturbance daily variation S, varies 
in the range 40 to 300 y depending on the magnetic activity 
and the latitude of the point. 

In polar zones, according to the data of spherical analysis, 
currents of high intensity and extremely inhomogeneous struc- 
ture flow in comparatively narrow bands. The intensity of 
the field of such current systems will vary in inverse propor- 
tion to the distance from the current system. Hence, in high 
latitudes we can expect current system effects reaching sev- 
eral thousand y. 

The amplitudes of the fluctuations of magnetic elements 
during magnetic storms have a very wide range. The maxi- 
mum effect of the current system of a magnetic storm, accord- 
to Chapman’s estimates (4), may reach 27,000 y. 

Thus, it becomes obvious that a satellite magnetometer 
must possess a very wide range and high sensitivity. If the 
velocity of the satellite is 8 km per sec, the instrument must be 
able to measure magnetic field changes occurring at a rate of 
more than 300 y per see. 

Another important characteristic of the apparatus is the sta- 
bility of its zero. A comparison of the magnetic data ob- 
tained at the beginning and end of the period of operation of 
the magnetometer can only be made if the zero has high sta- 
bility, or if its variation, studied in the laboratory, is regular 
and lawful. A zero shift of ~ 50 y per day can be regarded as 
quite admissible, since with this instability of the zero the in- 
determinacy in measuring the field values (after introduction 
of the appropriate corrections) will not exceed the errors of 
terrestrial magnetic maps. The average error of modern 
magnetic maps, as we know, is 200-300 y. 


2. Selection of Measuring Apparatus 


As we know, Earth’s magnetic field at any point in space is 
characterized by a magnitude and a direction. The best op- 
portunities for further analysis and comparison with the 
ground values of the magnetic field will naturally be provided 
by an instrument which measures the magnitude and direc- 
tion, or particular components, of the magnetic field. The 
majority of magnetically sensitive detectors measure the field 
components directly. 

Since satellites in motion in orbit are continuously changing 
their orientation in space (and relative to the magnetic field 
vector), the measurement of individual components, the 
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vertical component, for instance, raises considerable technical 
difficulties. 

For this reason magnetometers measuring the scalar value 
of the field are most suitable. Such instruments are nuclear- 
resonance magnetometers and total-vector magnetometers 
with saturable-core detectors. 


Nuclear Resonance Magnetometer 


The nuclear resonance method of measuring magnetic 
field intensity is based on the use of the phenomenon of free 
precession of protons in an external magnetic field. The 
free-precession frequency of protons, which possess a mag- 
netic and mechanical moment, is determined, as we know, by 
the Larmor relationship w = y,H, where y, is the gyro- 
magnetic ratio and H is the magnetic field intensity. 

By measuring the free-precession frequency of protons in 
Karth’s magnetic field we can determine the field intensity. 
The method of observing free precession of protons in Earth’s 
magnetic field [as suggested by Packard and Varian (5)] is 
extremely simple. A sample (liquid with high proton con- 
tent in the field Ho) is subjected for a short time to a strong 
field H (about 100 oe), which is then sharply cut off. The field 
H is created in an excitation coil surrounding the sample, 
and is directed approximately perpendicular to the field Ho. 
Under the action of field H the sample acquires a macroscopic 
magnetism, the intensity of which is equal to J, = x,//, 
where x, is the nuclear susceptibility. When the field H 
is cut off, the macroscopic magnetic moment arising in the 
sample begins to precess freely around field Ho with a fre- 
quency » = y,f>. The magnitude of the macroscopic 
magnetic moment of the sample gradually decays, but the 
relaxation time is about 3 see and is adequate to allow meas- 
urement of the frequency of the alternating voltage induced 
by the precessing moment of the sample in the signal coil, 
which consists of an excitation coil connected in this case to 
an amplifier. The alternating emf induced in the coil will 


be 
E = sin? 


where 


ll 


constant depending on the coil parameters 
6 = angle between field H and Ho 

t = time elapsing from cutoff of the field 

T. = relaxation time 


Field measurement based on nuclear induction possesses 
outstanding advantages over other methods: 

1 Field measurement reduces to frequency measurement. 

2 The precision of measurement depends on the design 
parameters of the detector and the channels shaping the 
signal. 

3 The magnetometer readings are given in absolute 
quantities. 

4 The detector and signal shaping channels are free in 
principle from zero creep; the precision of measurement is 
determined only by the accuracy of the frequency standard. 

5 The results of measurement for a constant position of 
the detector do not depend on the orientation of the detector 
in the field being measured. 

The detector signal, amplified by an audiofrequency ampli- 
fier, can be transmitted by telemetry and accurately meas- 
ured at a distance (the signal can be used directly for modula- 
tion of the transmitter). The precision of measurements 
with proton magnetometers is an observatory is estimated as 
1 y (10> oe) and is limited only by the accuracy of deter 
mining Y>. 

When this method is used on moving objects certain diffi- 
culties arise, and these complicate the measurements and 
increase their error: 

1 Although the precession frequency does not depend on 
the orientation of the excitation coil in relation to the field 
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being measured, the signal amplitude is proportional to sin? 6 
and at small 6 will be close to zero. 

2 When the coil rotates together with the moving object 
at angular velocity ¢ about an axis perpendicular to the coil 
axis, the field is measured with an error AH = +3.7 ¢ 
(AH is measured in gammas, and ¢ in radians per sec) (6). 

3 Inhomogeneity of the magnetic field within the volume 
of the magnetometer detector greatly reduces the relaxation 
time 7.. Even with field gradients 3-4 y per em the signal 
decays so rapidly that measurement of the field is practically 
impossible. 

Another limitation of this method arises due to the ex- 
tremely high sensitivity to interference from extraneous 
alternating voltages. The best properties in this respect are 
possessed by toroidal detectors, which in the ideal case 
should be free from the effect of extraneous interference if the 
latter is uniform within the detector volume. In fact, even 
toroidal detectors are not free from interference from alter- 
nating voltage of either low or high frequency. 

In ground magnetometers designed for field measurements 
in a very narrow range (in this case the detectors are located 
in sites free from interference), an adequate signal to noise 
ratio can easily be attained with comparatively small polar- 
izing fields. For instance, with the amplifier band 100 cps, a 
signal to noise ratio of 20 can easily be obtained with a 
polarizing power of about 20 w. 

Ina rocket magnetometer (7) which had a relatively narrow 
measurement range (hence, a fairly narrow amplifier band 
also), the polarizing power had to be increased to 70 w in 
order to obtain a signal to noise ratio of 12. 

To attain the same ratio in a magnetometer intended for 
magnetic field measurements from a_ satellite in the field 
range 15-55 y & 10° and, hence, with an amplifier band of 
around 1600 eps, would require an increase in the polarizing 
power to 200 w (7). 

In the case of measurements from a satellite, the proton 
magnetometer becomes more complex by the necessity, 
when the number of ground stations is limited, of a memory 
device and a generator of reference marks, which must be 
remembered together with the magnetometer signals. Lastly, 
in the immediate proximity of the detector, there must be no 
sources of extraneous magnetic fields with gradients 3-4. 
perem, 

In practice all these requirements can be met only in an 
autonomous experiment and with a container of special de- 
sign. 


Total-Vector Magnetometer With Saturable-Core Sensors 

The first measurements of the Earth’s magnetic field on a 
rocket (8) were made by a magnetometer with saturable-core 
detectors. In the arrangement used, signals from three 
mutually perpendicular detectors were fed to the input, and 
at the output there was obtained a value proportional to the 
square of the total vector. The dimensions and power 
supply of such an instrument are adequately small. How- 
ever, this arrangement is applicable only for measurements 
in a limited range of field variation and a particular, though 
not critical, orientation in space, as was the case in the cited 
rocket flight. In measurements from a satellite, such a 
simple arrangement would give only extremely coarse results. 

The absolute accuracy of squaring elements in a wide 
range of signals is not better than 5% (9). 

Magnetie field measurements in a wide range with high 
sensitivity and the transmission of these data by telemetry 
with an accuracy much better than the accuracy of tele- 
metry are possible if a compensation method of measurement 
is used. 

The compensation method in a three-component mag- 
netometer in a rotating satellite can only be used if one of the 
magnetometer detectors (the measuring one) is automatically 
kept oriented in the direction of the total vector of the mag- 
netic field. 
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An analysis of the actual operating conditions of the mag- 
netometric equipment on the third artificial Karth satellite 
Jed to the conclusion that experimental measurement of the 
magnetic field was perfectly realizable in practice by the use 
of a magnetometer with saturable-core detectors, orienting 
itself along the total vector. The following circumstances 
were taken into consideration: 

1 With the satellite containing a large amount of scien- 
tific and telemetry equipment, which inevitably creates 
magnetie and electrical interference, observation of the signal 
from the nuclear-resonance magnetometer detectors would be 
extremely difficult. 

2 To attain a measurement accuracy of 0.05-0.1 % in a 
wide range of measurable fields a compensation method of 
measurement would have to be employed. 

3. Since the spatial orientation of an orbiting satellite is 
continuously changing, the compensation method cannot be 
used unless the magnetometer detector is automatically 
oriented in the direction of the field. 

4 By the employment of semiconductors and _ efficient 
circuitry, it is possible to construct a light magnetometer 
which has a low power requirement, which can automatically 
orient itself along the total vector, and which can function at 
any orientation of the satellite and at any magnetic latitude. 

5 With such a magnetometer it is possible not only to 
measure the value of Earth’s magnetic field, but also to 
determine the orientation of the satellite in space and the 
nature of its motion relative to its own center of mass; such 
sata are necessary for interpretation of the data of other ex- 
periments. 


3. Magnetometer of the Third Artificial Earth 
Satellite 


General Description of Measurement Principles 


The magnetometer of the third satellite was an automatic 
device, which continuously performed two operations: 

1 Orientation of a special measuring detector in the direc- 
tion of the magnetic field vector at any orientation of the 
satellite. 

2 Measurement of the magnetic field intensity by means 
of this detector. 

Orientation in the direction of the field and field measure- 
ment were effected by means of sensitive elements, the opera- 
tion of which was based on the same physica! principle. The 
detector signals were merely used in different ways for the 
orienting and measuring channels of the magnetometer. 

A schematic diagram of the magnetometer is shown in 
Fig. 1. It consisted of four main parts: Measuring channel 
circuit, circuit of orienting channels, mechanical orientation 
unit, and automatic control elements and power supplies. 

The measuring and orienting channels contained several 
similar functional elements: Detectors, selective amplifiers, 
phase-sensitive rectifiers. These clements served for con- 
verting the permanent magnetic field signals into d-c electric 
signals of adequate power. 

The elements sensitive to Earth’s magnetic field consisted 
of saturable-core detectors of the second-harmonic type. 
In the simplest case, a saturable-core detector consists of a 
plate or wire of material of high magnetic permeability (per- 
malloy) and with primary and secondary coils wound on it. 

When the wire is magnetized by a field H = Hy + Hy, sin 
wt (here Ho is Earth’s magnetic field, and H/,, sin wt the auxil- 
iary sinusoidal excitation field of frequency w, magnetizing 
the wire to saturation twice during one period of field varia- 
tion) a voltage Us» of doubled frequency arises in the secon- 
dary winding. The amplitude of the voltage U2» of doubled 
frequency is proportional within certain limits to the intensity 
of the field being measured. When Hy) = 0 the voltage Usa 
is also zero. When Hy changes sign, the phase of the second 
harmonic voltage also changes. 
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A thin and adequately long permalloy wire possesses 
another remarkable property: It is easily magnetized by a 
field directed along its axis, and it is practically impossible 
to magnetize it with a field perpendicular to its axis. Thus, 
saturable-core detectors are unique devices for converting 
a permanent magnetic field signal into an electrical alternat- 
ing voltage signal of frequency 2 w. Such detectors possess 
two important properties: Direct sensitivity to magnitude 
and sign of the magnetic field and selectivity in relation to 
direction of field. 

The total-vector magnetometer of the self-orienting type 
employs three mutually perpendicular detectors mounted on 
a stage in the special orientation unit. One of these detectors 
is the measuring detector, and the other two, fixed in the plane 
of the stage, serve to orient the stage perpendicular to the 
total vector. 

The excitation voltage of frequency 2 ke per sec, providing 
periodic magnetization of the cores, is taken from a special 
generator of frequency 2 ke per sec. The detector signals 
(second-harmonic voltage) are amplified by selective ampli- 
fiers tuned to frequency 4 ke per sec, after which they pass to 
the input of the phase-sensitive rectifiers. The constant 
component of the output voltage of the phase-sensitive recti- 
fier of each channel is linearly dependent within certain limits 
on the magnitude of the permanent magnetic field being 
measured, and its sign depends on the sign of the field. 

Subsequently the systems of the measuring and orienting 
channels differ substantially. The d-c signals from the phase- 
sensitive rectifiers of the orienting channels are fed to the 
inputs of the magnetic amplifiers of the servo system of the 
orienting channels, where they are converted into an alter- 
nating voltage of frequency 400 cps. After power ampli- 
fication, this voltage is imposed on the control windings of the 
low-inertia motors of the mechanical orientation unit. 

The motors rotate the stage, in whose plane the orienting 
detectors are mounted, until it comes into a position perpen- 
dicular to the total vector of the magnetic field. When this 
occurs, the second-harmonic signal from the orienting de- 
tectors will be equal to zero, and the motors will stop. 

The measuring detector will be oriented along the total 
vector, and the signal from it is used for measurement of the 
magnetic field intensity. The direct current passing to the 
compensation winding on the detector of the measuring 
channel gives a measure of the magnetic field and its varia- 
tions. 

The main part of the magnetic field is compensated from a 
special current source of high stability. The compensation 
current flowing through the detector winding can be varied 
discretely (in definite steps) by means of the automatic 
range switch. Continuous field variation in limits somewhat 
exceeding the individual compensation steps of the auto- 
matic range switch are automatically compensated by the 
introduction of large negative feedback; compensation is 
effected by the current obtained after the phase-sensitive 
rectifier in the measuring channel from the signal produced 
by the uncompensated part of the magnetic field. The feed- 
back current compensates almost the whole field, with the 
exception of a very small fraction which serves for maintain- 
ing the compensation current. The use of a compensation 
method based on negative feedback greatly increases the 
stability of the magnetometer. 

After power amplification by the d-c amplifier the self- 
compensation current passes to the windings of the polarized 
relays connected in series with the load resistance of the meas- 
uring channel. When the self-compensation current ex- 
ceeds the limit for a particular range of continuous compensa- 
tion of the field, the polarized relays close the appropriate 
contacts of the automatic range switch. In this way the 
total measurement range is widened, and field measurements 
can be made in a broad interval of values with constant high 
sensitivity. 

The position of the movable platform relative to the mag- 
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Fig. 1 Block diagram of SG-45 magnetometer 


a. Orienting and measuring channels: A,—amplifiers for fre- 
quency 2w; PR—phase-sensitive rectifiers; G.—generator of fre- 
quencyw; SV—synchronization voltage of frequencyw; DC A—d-c 
amplifier; \/ A—magnetic amplifiers for frequency w; Gw,—gener- 
ator of frequency w:; R:—RP-5 relays; M:—motor for fork 
orientation; /.—motor for stage orientation; RS—electromagnetic 
range switch; RD—reducing gears; P—potentiometers; /SE— 
stage with magnetically sensitive elements; 1 W—control windings 
of motors; (P—constant-phase windings of motors; R7S—telem- 
etry. Circuits: 1—signal from measuring detector (frequency 2w); 
2—excitation of detectors by frequency w; 3—signal from orienting 
detectors (frequency 2w); 4—damping of self-oscillations of stage; 
5—damping of self-oscillations of fork; 6—synchronization voltage ; 
7—feedback; 8—field compensation (—H); 9—field compensation 
(+H); 10—stepwise compensation; frequencies w = 2000 cps, w: 
= 400 cps 

b. Automatic control components and power supplies: R»—relay 
for cutting in circuit of measuring channel; /A)—relay for cutting in 
circuit of orienting channels; P.\/—circuits from programming mech- 
anism; PS—power supply circuits 


netic field vector is determined by the position of the moving 
contacts of the two potentiometers, the frames of which 
are fastened to the platform body, and the wipers of which are 
connected by a gear train to the spindles which cause the 
rotation of the stage relative to an axle mounted in the fork, 
and the fork itself relative to a perpendicular axis. One 
turn of the potentiometer corresponds exactly to one turn 
of the stage or fork. 


Description of Operation of Individual Units and 
Construction of SG-45 Magnetometer 
Schematics of the magnetometer are shown in Figs. 2 and 3. 


Magnetometer transducers 

A harmonic analysis of the output voltage on the trans- 
ducer winding (when the core is magnetized by an excitation 
field H,, sin wt and a constant field Ho) shows that the output 
voltage contains a discrete spectrum of odd and even har- 
monies, the amplitudes of which are functions of the fields 
H, and H,,. The even harmonies are an odd function of the 
field Hy and an even function of the auxiliary excitation 
field H,,; a change of phase of H,, does not alter the phase of 
the even harmonics. On the other hand, the odd harmonics 
are an even function of the field //) and an odd function of 
field H,,. 

The detector for measuring Earth’s magnetic field consists of 
two cores (permalloy strips) contained in coil forms. On each 
of the forms are wound the excitation coils W,, which are 
connected in series and are wound in opposite senses (Fig. 4). 
On top of the coils W, in the detectors of the orienting chan- 
nels are wound signal coils W, connected in series and wound 
in one sense. With this arrangement the odd harmonics 
will be compensated in the coil W, (if the half elements of 
the detector are identical). In the measuring channel 
detector one coil W, is wound round both strips. 

Besides these coils the feedback coils Wy are wound on the 
frames of the orienting channel detectors, and a compensation 
coil W,. is wound on the measuring channel detector. 

Since the cores are constructed in the form of thin plates, 
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Fig. 2. Schematic of measuring channel of magnetometer 
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Fig. 3 Circuitry of orienting channels 


the length of which is much greater than the cross-sectional 
perimeter, or in the form of a fine wire, the demagnetization 
coefficient of the cores along the major axis is thousands 
of times smaller than the demagnetization coefficient in the 
transverse direction. Such a core is not magnetized by the 
whole field, but only by the projection of the field on its major 
axis. For this reason the voltage of the second harmonic in 
the coil of a detector set perpendicular to the field is zero. 
When the detector deviates to one side or another a voltage 
of corresponding phase arises. 


Viagnetometer generator 


The auxiliary voltage for excitation of the transducer is 
provided by a generator of sinusoidal oscillations of fre- 
queney 2000 eps. 

For stable operation of the magnetometer it is extremely 
important that the generator voltage contains no even har- 
monies, particularly the second harmonic. For this reason 
the generator is connected in a push-pull arrangement, 
and the detector is supplied after passage of the output voltage 
through low-pass filters C3, Cos, Cos). 

The generator consists of two stages: A master generator, 
incorporating P2B germanium triodes (K7’s and AT’,) and a 
buffer power amplifier, consisting of germanium triodes of 
the type P3V (AKT, and AT;) in a common-base circuit. 
Che tank circuit of the master generator (7’r. 5) utilizes an 
isifer core. The voltage from the buffer stage transformer 
Tr. 4) passes to the input of a low-pass filter and to the drive 
f{ two power amplifiers incorporating the triodes P2B 
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Fig. 4 Connection of magnetometer coils. Thin lines—W, 
windings, thick lines—-W, windings 
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(KT; and KT. on Fig. 3), the voltage of which is used for 
synchronization of the phase-sensitive rectifiers of the 
orienting channels. The synchronization voltage of the 
phase-sensitive rectifier of the measuring channel is supplied 
from the secondary winding of transformer Tr. 4. The 
magnetometer detectors require a current of 75 ma at 4 v. 
The 2000-cps generator consumes from the battery a current 
of 60 ma at 25 v. 


Magnetometer amplifiers 


With cores of length 40 mm and a practical and convenient 
number of turns of coil W,, the sensitivity of the detectors is 
approximately 10 wv per y. Due to the difference in 
thickness of the cores, somewhat different shapes and sizes 
of frames, the distributed capacitances, and leakage fluxes, 
the signal coil contains odd harmonics whose potential may 
reach 50-100 mv. 

For suppression of the odd harmonics and for power and 
voltage amplification of the useful signals, three identical 
selective amplifiers, tuned to 4000-cps frequency, are used 
in the circuit. The input element of each amplifier is a 
passive band filter with maximum attenuation at frequencies 
2000 and 6000 eps (first and third harmonics of excitation 
voltage). The purpose of the filters is to provide a favorable 
signal to noise ratio and thus to prevent the possible ap- 
pearance of spurious even harmonics due to nonlinear effects 
in the amplifier. 

The amplifiers each comprise three stages consisting of type 
P6D triodes in a common-emitter circuit. The tank circuits 
are made with “Oksifer 200” cores. 

The total voltage amplification factor is 20,000. The 
bandwidth at level 0.7 is 350 eps. The filter attenuation for 
the first harmonic is 54 db. The relatively broad band is 
necessary for stable operation of the magnetometer. The 
supply voltage of the amplifiers is 11 vy. All three amplifiers 
are tuned to attain precisely identical phase characteristics. 


t 
uow 


Fig. 5 Addition of even and odd harmonics 
displaced in phase 


Phase-sensitive rectifiers 


The action of the phase-sensitive rectifier used in the mag- 
netometer is based on the fact that when odd and even 
harmonics displaced in phase are combined, there appears a 
voltage wave of distorted shape, the negative and positive 
peaks of which have different heights (Fig. 5). The direction 
of the excess peak depends on the phase of the even harmonic 
and the difference in peaks depends on the signal magnitud 
and the phase shift between the harmonics. By changing th: 
phase of the excitation voltage of the detector, the phas 
shift between the synchronization voltage (2000 eps) and th: 
useful signal voltage (4000 eps) can be made equal to 90 
deg. 

The difference in amplitudes of the half-waves, which in 
certain limits is proportional to the field being measured, is 
measured by a detector consisting of a differential peak 
voltmeter. 

The circuit (see Fig. 2) consists of an input transforme: 
(Tr. 2), the primary winding of which constitutes the res 
onance circuit of the last amplifier stage, and the secondary 
winding (center-tapped) is connected through two germaniun 
diodes to a system of capacitors and resistors (Ci4, Cis, 
Ri, Ris, Rig, Ro). To the central arm of this balancing net- 
work is connected the transformer winding, through which 
the synchronization voltage of frequency 2000 eps is fed. 

In this case the mean values of the voltage on the capacitors 
will be equal if the amplitudes of the transformer voltagi 
half-waves are equal, and will differ if the half-wave am- 
plitudes are different. 

In the latter case a current will pass through the load 
resistance of the phase-sensitive rectifier. The sign and 
magnitude of this current will depend on the value of the 
difference voltage which, in turn, will depend on the sign of 
the magnetic field being measured. 


Compensation circuits of magnetometer measuring 
channel and output to telemetry 

As we already mentioned, the measuring detector operates 
by the null method, the main part of the field compensated 
by means of a stable field source, and the small continuously 
changing part of the field compensated by the introduction 
of large negative feedback. 

The voltage of the phase-sensitive rectifier is fed through 
filters (L7, Ly and Cis) to the feedback winding. By changing 
the number of turns of the feedback winding, the sensitivity 
of the instrument can be varied within wide limits. The 
self-compensation current changes sign when the magnetic 
field changes sign. 

A voltage up to 6v of one sign can be applied to the 
telemetry system. To satisfy this requirement the output 
to the telemetry is effected in the following way. The self- 
compensation current is amplified by a d-c amplifier, consist- 
ing of triodes P6V (KT, A7),) in a balancing network. 
The d-e amplifier load consists of the control windings of two 
polarized relays of the type RP-5, the purpose of which will 
be clear from the following. The windings are connected in 
series, each of the windings being shunted by resistors in 
such a way so that the relays make contacts when the current 
is | ma. Connected in series with the relay windings are 
resistors Rey and Res, chosen so that at a current of 1 ma the 
voltage drop on the whole circuit will be 12 v. This voltage 
also changes sign when the field sign changes. 

The voltage from the ends of the circuit of the relay 
windings and resistors Ros, Reg is fed to a circuit consisting of 
two diodes D3 and D, and load resistors Ryo and Rs; connected 
as shown in the diagram. The voltage from resistors Ry 
and 3; is fed through emitter followers K7; and KT’; to the 
alternately operating telemetry channels MF, and MF», 
one measuring the positive changes of field, and the other the 
negative. In this case a voltage of one sign, varying from 
0 to 6 v, passes to channels MF; and MF». 
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The automatic compensation of the main part of the field 
(40,000 + 24,000 y) is provided by a circuit fed from an 
extraneous source of stable current. The compensation 
coil of the measuring detector is connected through a special 
switch with a system of parallel resistors which are chosen 
so that the compensation field in the coil can be changed by 
steps of 3000 y for each successive switch of the ranges 
(rig. 6). To prevent the resistance of the whole system 
changing when the compensation coil is shunted, the com- 
pensation coil is connected simultaneously by means of the 
~econd plate of the special switch to a system of series- 
mnected resistors chosen so that the total resistance of the 
ompensation circuit remains constant. 
The current is fed into the compensation coil through a 
ad resistor Rian from a battery of zinc-mercury cells 
f 60 v and capacity 3 amp-hr. Connected in_ series 
ith the compensation coil are two filter chokes and an 
dditional resistor, which is chosen so that with a nominal 
irrent of 5 ma through the coil the difference in potentials of 
oints a and 6 is equal to the voltage of a standard battery 
hich, as shown in the diagram, is connected as a buffer to 
he load. The standard battery consists of a battery of 
inc-mercury cells of voltage 2.7 v and capacity 30 amp-hr. 
his battery supplies power to the compensation coil only 
vhen the voltage of the operating battery begins to fall. 
\n analysis of the operation of the circuit (10) showed that 
he stability of the current in the compensation circuit was 
ligher, the greater the ratio Ro,/Rs: (where R,; is the internal 
resistance of the standard battery) and the more stable the 
voltage of this battery. 
This circuit allows the compensation current to be held 
constant without any regulation and with a known degree of 
accuracy. 


Fig. 6 Circuit of range switch 


The whole system of automatic compensation works in the 
following way. When the voltage drop on the relay windings 
reaches 12 v (6 v on the resistor from which the signal for the 
telemetry channel is taken), the right or left contact of 
relay RP-5, depending on the sign of the current, is made. 
The sign of the current depends on whether at this value of 
the magnetic field intensity there is overcompensation or 
undercompensation of the field in the coil of the measuring 
channel detector. 

The making of the right or left pairs of contacts of relay 
RP-5 leads to closure of one of the range switch electro- 


4 


Fig. 7 Diagram of orientation unit 


M,—DID-0.5 motor (for rotation of fork); .4.—DID-0.5 motor (for rotation of stage); Pi,,.—inclination potentiometer; 
Piec—declination potentiometer; 1—bevel gears; 2—-gears; 3—pinions (small driving gears) 
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magnets, which then shifts the movable contacts of both 
plates of the range switch through the necessary number of 
ranges until the relay RP-5 opens. This corresponds to field 
compensation within the range of continuous variations. 
To eliminate hunting of the range switch, feedback propor- 
tional to the rate of change of the signal is introduced by 
means of the differentiating circuits Cig. and Rgo, Cio. 


Mechanical orientation unit 


Fig. 7 shows a diagram of the orientation unit. The stage 
carries the three mutually perpendicular detectors, one of 
which is the measuring detector and the two others, lying in 
the plane of the stage, are the orienting detectors. 

The stage can rotate round an axle mounted in the fork. 
In its turn the fork can rotate round a perpendicular axis. 
The rotation is effected by means of two gear trains. The 
driving gears are fastened on the driving spindles, which 
transmit the rotation through reducing gears from two low- 
inertia motors of the type DID-0.5. One of the motors 
(Mz) effects rotation of the stage round an axis perpendicular 
to the plane of the magnetic meridian; the other (.V/,), 
the azimuthal, brings the system into a position where the 
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Fig.8 Measurement error as function of stage orientation 


Fig. 9 Photograph of magnetometer 


1—Orientation pack; 2—electronics pack 
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stage with the orienting detectors is perpendicular to the 
total vector of the field. 

The voltage is transmitted to the detector coils by means of 
sliding conductors. In view of this construction, the ori- 
entation unit allows measurements when the movable plat- 
form is repeatedly rotating. The driving spindles of the 
orientation unit carry pinions which engage with the gears 
of two ring potentiometers, which receive a voltage from a 
constant voltage source of 6v. The voltage taken from the 
movable contacts of the potentiometers depends on the 
position of the contacts and characterizes the mutual position 
of the movable platform and the stage on which the detectors 
are attached. 

Fig. 8 shows the values of the signal on the orienting 
detectors and the error in the measuring detectors as functions 
of the angle of misalignment. The misalignment signal, 
which arises in the orienting detectors when the normal to the 
stage deviates from the direction of the total vector of the 
field by an angle a, is proportional to sin a and the corre- 
sponding change of the signal in the measuring detector is 
proportional to (1 — cos a). With a 1-deg error in orienta- 
tion of the stage, the error of the measuring detector does 
not exceed 4 ¥ (in field measurements at Earth’s surface). 


Magnetometer servo system 


The electronic units of the orientation pack are designed for 
power and voltage amplification of the misalignment signals 
generated by the coils of the orienting detectors, and the shap- 
ing of the signals for controlling the motors of the mechanical 
orientation unit. 

As we already mentioned, the analogous elements in the 
orienting and measuring channels are the detectors, selective 
amplifiers phase-sensitive rectifiers. Otherwise the 
systems differ substantially. 

The signals from the phase-sensitive rectifiers of the 
orienting channels arrive at the input of magnetic amplifiers 
MA, and M Az connected in a push-pull arrangement. The 
magnetic amplifiers operate at a frequency of 400 eps. The 
400-cps voltage is supplied from a special generator consisting 
of a master stage incorporating germanium triodes P2B 
(KT\g, and KT\s), and two power amplifiers connected in a 
push-pull arrangement with germanium triodes P3V in a 
common-emitter circuit (A 7, AT and The 
output windings of the transformers of these amplifiers are 
connected to the constant phases of the motors of the ori- 
entation unit. From these amplifiers a voltage is fed to the 
excitation coils of the magnetic amplifiers. The output 
winding of each magnetic amplifier is connected through a 
matching transformer (7'r. 15 and Tr. 18) to the inputs of the 
power amplifiers supplying the control windings of the 
motors. The power amplifiers of the controlling phases 
consist of triodes P3V in a push-pull arrangement. The 
voltages arriving at the windings of the constant and con- 
trolling phases of the motors are shifted through 90 deg. 
Thus, the low-inertia motors of the orientation unit operate 
at a different frequency from the excitation frequency of the 
detectors and the frequency of the useful signal, and this 
prevents the measuring system from being affected by the 
motors, which are located in the immediate vicinity of the 
detectors. 

Since the servo system contains no narrow-band selective 
elements, its operation is affected very little by the instability 
of the supply voltage. The servo system rate of handling is 
about 45 deg per sec. 

The magnetometer is provided with separate switching-in 
arrangements for the orientation system and the measuring 
channel; this is achieved by the use of two independent 
relays DP-10. The instrument is switched on by a special 
programming mechanism; the orientation system switched 
in first, and then the measuring channel, thus eliminating 
the possibility of overloading the latter. 
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Construction of magnetometer 


The magnetometer consists of three packs connected by 
means of plug connectors (Fig. 9): Orientation pack, elec- 
tronics pack, and power pack. 

The orientation pack is of light construction and consists of 
three main parts: A front part with the magnetically sensitive 
elements and current collectors; a central part with the 
reducer, potentiometers and motors; and a rear flange 
carrying the plug connector (Fig. 10). 

The low-inertia motors of the type DID-0.5, which effect 
the setting of the measuring detector along the total vector, 
can develop up to 13,000 rpm under no load, whereas under a 
load of 3 gm per cm the motors give 7000 rpm, corresponding 
to a rotation of the stage or fork at a rate of 14 rev per sec. 
Repeated clockwise and counterclockwise rotation in two 
mutually perpendicular directions is effected by the use of 
bevel gears. 

The motors (the only magnetic component in the pack) 
are placed at a considerable distance from the detectors, so 
that the latter will not be affected. The transmission of 
rotation to the head end is effected by long driving spindles 
connected to the corresponding axles by flexible couplings. 

The orientation pack is reliably shockproofed, thus insuring 
operation of the instrument under vibration. 

The electronics pack is assembled in the form of a stack. 
The bottom plate carries the 400-cps generator, the servo 
system and power amplifiers. The second bottom plate 
bears the 2000-cps generator, the low-pass filter and power 
amplifier; the next plate, the three selective amplifiers and 
the phase-sensitive rectifiers. The top plate carries the 
range switch, the two cut-in relays DP-10, the relay RP-5, 
voltage distributing plate, chokes and other components. 

The range switch is provided with differential gearing, 
which allows the kinematic uncoupling of the two rotations. 
For greater electrical reliability, reduction of contact re- 
sistances and their variations, the bars and contacts of the 
switch are twinned and palladium coated. 

All the connecting cables of the magnetometer are attached 
to the top cover of the electronics pack. The plates are all 
held together by four rods. When the rods are removed, the 
electronics pack can be spread out without unsoldering the 
plates. The electronics pack is completely enclosed in a 
casing. Four bolts are provided in the base of the elec- 
tronics pack for attachment. 

The weight of the electronics unit is 12.5 kg. 

The power pack contains 12 batteries, the leads from which 
are connected to the plug connector. By the provision of 
separate supplies for the generators, amplifiers and servo 
system, the effect of the channels on one another and inter- 
ference between the individual electronic units are removed 
with no significant increase in weight or dimensions. It was 
of particular importance to separate the supplies of com- 
ponents consuming great power, but not requiring high 
stability, from the supply of components requiring little 
power but a supply voltage of adequate stability. 


Metrological Characteristics of Magnetometer, Nature of 
Record and Its Decoding 


All the magnetometer specimens were carefully tested for 
their basic metrological characteristics: Sensitivity, stability 
of zero in time, temperature coefficient, rate of handling of 
signals, and dynamic error. 

The readings of experimental specimens of the magnetom- 
eter were given absolute values by comparing them with 
observatory magnetometers. The zero of a satellite magne- 
tometer specimen was determined absolutely by means of a 
proton magnetometer. Measurements by both instruments 
were made at the same point on a special stage. 

Information on the operation of the magnetometer on 
board the satellite was transmitted through five telemetry 
channels, two of which (KY, and KY.) transmitted informa- 
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Fig. 10 Photograph of orientation pack with 
casing removed 


tion about the position of the orienting detectors; the other 
three transmitted information on the positive and negative 
changes of the magnetic field (WF, and MF.) and on the 
position of the range switch (RS). The magnetometer 
sensitivity was such that the total range of field measure- 
ments in each of the two telemetry channels was +2400 y. 
Channel MF, transmitted positive variations of field and 
channel MF», negative variations, in relation to the initial 
level. When there were readings in channel ./F, the readings 
in channel MF, were zero, and vice versa. The total range 
of the two channels (4800 y) exceeded the range of one step 
of the discrete compensation of the range switch (3000 y). 

The nature of the voltage traces from the telemetry 
channels in the case of linear variation of the magnetic field 
is shown in Fig. 11. As the graphs reveal (Fig. 11a), in- 
crease of field is recorded by an increase of voltage in channel 
MF,. When the limiting value is reached, range switching 
occurs; it is marked by a simultaneous voltage kick in 
channel RS. At the same time there appears a voltage 
in channel MF», this voltage corresponding to the difference 
between the value of the range of channel .7F, and the value 
of the range of one step of discrete compensation (this 
difference was 600 y). With further increase in field the 


signal on MF, falls to zero, then appears in channel WF;, 
and so on. 

When the field decreases, the nature of the change in the 
form of the trace is similar to that just described. The 
sequence of changes in channel RS is reversed; the voltage 
increases in channel WF:, and decreases in channel JF; 
(Fig. 11b). 
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Thus, the range switching is recorded in the indications of 
three telemetry channels. 

The intensity of Earth’s magnetic field from the readings 
of the three telemetry channels is given by the formula 


T = Ty) — (n — m) K + & (@ — ao) — & (8B — Bo) 


where 
To = magnetometer zero determined by comparison 
with the proton magnetometer on the reference 
stage 


Y 
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Fig. 11 Voltage traces from telemetry channels 


a—Field increase; b—field decrease; c—voltage traces of detectors 
controlling angle of rotat on 
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Fig. 12 Variation of compensation current with time 
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No,Qo,80 = readings from channels RS, MF, and MF, 
on the reference stage 

n,a,8 = readings from channels RS, WF\, WF: in the 
measurements 

€1,€2 = sensitivities of channels .1/F, and 

K = range value 


Fig. 1le shows the voltage traces of the angular position 
control detectors for a constant rate of rotation. The 
orientation unit was designed so that rotation along the 
longitudinal axis (rotation of fork) caused a change in the 
position of the stage on which the transducers were mounted. 
Hence, the indications of the angular position control potenti- 
ometers (K.V.and KY.) are connected by the relationships 


AV xy, => kAg 
AV xy, ko Ay 


The stability of the magnetometer zero depends on the 
stability of operation of the electronics units and the stability 
in time of the current compensating the main part of the 
field. The nature of the compensation current. variation 
over a long period of time is shown by the graph given in 
Fig. 12. This graph pertains to an arrangement in which the 
standard battery (B,,) had a capacity of 12 amp-hr and a 
voltage 4v. The mean current variation in this arrangement 
was equivalent to a zero shift of 135 y in 21 days. As we 
have already pointed out, the stabilizing properties of the 
cireuit depend on the internal resistance of the standard 
battery. The capacity of the satellite standard battery was 
increased to 30 amp-hr, and the voltage reduced to 2.7 v. 

Fig. 13 shows the compensation current variation for a 
circuit with parameters identical to those of the compensation 
circuit of satellite specimens. These observations were 
part of a series conducted to determine the temperature 
coefficient. The temperature varied fairly frequently. 
The readings were brought to one temperature. Despite the 
fact that the internal resistance of the standard battery was 
in this case reduced by a factor of nearly 4, the mean zero 
drift was 10 y per day. This can be attributed to the high 
temperature coefficient of the batteries and the thermal 
inertia of the system, and this resulted in a somewhat greater 
instability of the system. 

Fig. 14 shows the compensation current as a function of 
temperature. The measurements were conducted over a 
period of 80 hr in order to reduce the effect of thermal inertia. 
The temperature coefficient was variable. In a first approx- 
imation two values of the temperature coefficient can be 
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Fig. 14 Compensation current as function of temperature 
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Fig. 15 Graph of temperature variation from May 14-31 for 
three temperature detectors 


taken for evaluation of the temperature errors: 6 y per deg 
in the range 0-20 deg and 2 y per deg in the range 20-50 
leg. 

Owing to the reliable system of heat regulation, the tem- 
perature on the satellite varied within narrow limits, as Fig. 
15 reveals. 

The stability of the zero depends also on the stability of 
operation of the electronics units. The general variation of 
the zero of the entire magnetometer in continuous operation 
is apparent from a consideration of Fig. 16, which shows the 
magnetic field changes from the data of observatory magneto- 
graphs and from the magnetometer of the third satellite. 
The continuous curve on Fig. 16 shows the field variation 
AZ and AH from the readings of observatory magnetographs. 
The reading accuracy from the visual instrument in the latter 
case was approximately 10 y. The lower curve shows the 
difference of the first two curves. It illustrates the stability 
of the magnetometer during the first 12 hr of continuous 
operation. 

The curves reveal that an appreciable instability is mani- 
fested in the first few hours of operation of the instrument. 
Subsequently the scatter of the points of the difference curve 
are of a rather random nature, and in absolute value lie 
within the range of errors of reading from a visual instrument 
of class 1.5 %. 

A very important characteristic of a magnetometer mounted 
on a rotating satellite is the rate of handling of the signal. 
This applies both to the servo system and to the elements of 
the measuring channel. The servo system rate of handling 
was determined by measuring the time taken to set the stage 
in a position perpendicular to the meridian for a pre-set 
angle of deviation. The rate of handling at small mis- 
alignment angles was determined by similar measurements 
with a homogeneous sinusoidal field of known frequency 
applied perpendicular to the vector of the terrestrial field. 
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Fig. 16 Variation of magnetic field with time 


The current source was a generator of very low frequency 
oscillations. These measurements showed that the servo 
system was capable of handling a misalignment signal at a 
rate of 45 deg per sec with angular accelerations up to 130 
deg per sec. The range switch of the measuring channel 
allowed handling of a misalignment signal of 48 y xX 10° 
in 1-1.5see. 

The time constant of the output circuit to the telemetry 
Was approximately 0.08 sec. The occasional hunting of the 
range switch associated with such a time constant of the out- 
put circuit and with the rapid action of the range switch 
was eliminated by the introduction of negative feedback 
proportional to the rate of change of the signal. 

At nominal voltage and in the quiescent state the instru- 
ment consumed about 12 w. At maximum rate of operation 
the power consumed rose to 20 w. 


Investigation of Deviation and Its Elimination 


The orientation pack with the sensitive detectors was 
placed in a special compartment, as far away from the other 
apparatus and power supplies as possible. However, this 
still did not completely remove the effects of the magnetic 
components of the apparatus for other experiments. The 
elements creating most interference were the permanent 
magnets of the magnetic manometers, electrostatic fluxmeters 
and magnetic components of other instruments. The relative 
positions of the sensitive detectors of the magnetometer 
and the elements creating most interference are shown in 
Fig. 17. 
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Fig. 17 Relative positions of instruments on third artificial Earth satellite 


1,2—Magnetometer; 3—magnetic and ionization manometers; 4—electrostatic fluxmeter; 5—solar batteries 
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Special measurements .were made to determrine the mag- 
nitude of the total magnetic interference and its nature. The 
satellite with all its equipment was placed on a rig permitting 
rotation of the satellite around a longitudinal horizontal 
axis and around a vertical axis. An example of these meas- 
urements is shown in Fig. 18. The maximum value of the 
magnetic deviation was 3500 y. 

Exclusion of the magnetic deviation from the results of 
measurements was made possible by the following cir- 
cumstances: With the given value of magnetic deviation, the 
magnetic interference alters the field intensity but has 
practically no effect on the field direction; the satellite 
completes a regular precessional motion with a large angle of 
precession, the speed of the precessional motion being so 
high that the magnetic interference creates changes of field, 
the gradients of which differ substantially from the gradients 
of the terrestrial magnetic field at high altitudes. 

If H, is the vector of Earth’s magnetic field, H; the vector of 
the magnetic interference, and ¢ the angle between them, 
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Fig. 18 Deviation curve and readings of orientation channels 

with simultaneous rotation around longitudinal and vertical axes. 
Circles—direct motion; crosses—reverse motion 


then the scalar value of the total vector may be denoted by 
the formula 


IT. + T; = + T — 2T,T cos¢ [1] 


In view of the smallness of value 7’; this expression can be 
expanded in a series of powers of 7’ 


9 


T, = T, — Tose + 2 T, sin? ¢ [2] 


Taking only the first term, we obtain 
T, = T. — Tose 


Thus, the satellite magnetometer in practice measures Karth’s 
magnetic field and the projection of the vector of the magnetic 
interference on the direction of the magnetic field. Owing 
to the precession of the satellite, the angle ¢ is continuously 
changing in magnitude and sign. 

Fig. 19 shows a specimen of the magnetometer trace. The 
slow smooth variations of the terrestrial magnetic field are 
modulated by the magnetic interference, which periodically 
varies in magnitude and sign. By graphical or analytical 
averaging of the curves the effect of the magnetic interference 
can be excluded to a certain degree. As follows from formula 
2, the third term of the expansion is not eliminated by this 
procedure. Disregarding the third term can produce a 
maximum error of about 100 y. 
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1—Altitude variation of satellite with time; 2—magnetic field intensity from data of ground measurements; 3—magnetic 
field intensity from measurements on satellite ; 4—magnetic field intensity measured on satellite and corrected for deviation 
caused by satellite itself 
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Reviewer’s Comment 


The main body of the paper, Section 3, is a detailed, well- 
written record of the magnetometer experiment on the third 
Soviet satellite. The record leaves no doubt that the experi- 
ment was carefully planned and executed by competent 
people. There is a frank appraisal of the limitations of the 
instrument and of the expected accuracy of measurement. 

The magnetometer design is not a new one. This type of 
instrument, a magnetically stabilized (or oriented) saturable- 
core magnetometer, was used extensively for submarine 
detection during World War II and for geophysical pros- 
pecting (1). The authors reject the use of a proton precession 
magnetometer—inherently a more precise and_ reliable 
instrument—on the grounds that their satellite produced 
excessive electrical and magnetic interference. This problem 
is always a difficult one in using proton precession magne- 
tometers. One feels that by proper use of shielding and by 
mounting the magnetometer sensor on a long boom a proton 
instrument could have been used. 

An interesting comparison can be made between these 
first Russian magnetic measurements at relatively low 
satellite altitudes and the first U.S. experiment for the same 
purpose. The Vanguard III satellite was chiefly intended for 
magnetic field) measurements (2). A proton precession 
magnetometer was used in conjunction with a worldwide 
network of ground magnetic observatories. A measurement 
of the magnetic field at the satellite was obtained by radio 
command as the satellite passed near a magnetic observatory. 
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Measurements accurate to 1 y (107° gauss) were expected 
under ideal telemetry conditions. The weight of the mag- 
netometer, including the battery supply for magnetometer 
and telemetry, was about 9 Ib and the lifetime at the rate of 
50 measurements per day was expected to be some 3 months. 

The Russian experiment described in this paper was 
designed to achieve measurements accurate to about 0.05 % 
of the total field (15,000 to 60,000 y) or from 7 to 30 y at 
best. The weight of the electronics pack of the magnetom- 
eter was stated to be 12.5 kg. The orientation pack was 
said to be of light construction, and the battery pack consisted 
of 12 batteries of unspecified capacity. The power con- 
sumption of continuous operation was between 12 and 20 w. 
A total weight of 50 lb is probably a modest estimate for 
comparison. This one experiment on the Soviet. satellite 
weighed more than the complete Vanguard III satellite 
(approximately 20 Ib). 

A final comparison of the two experiments awaits complete 
publication of the measurement results. In summary, the 
Russian experiment appears carefully done, although the 
instrument was somewhat outmoded and cumbersome. 
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Investigation of Meteoric Particles 


by the Third Soviet Satellite 


| hoe ADDITION to the planets and their satellites, asteroids 
and comets, the solar system contains a multitude of 
smaller bodies known collectively as meteoric material. 

Up to now the bulk of our information about meteoric 
bodies entering Earth’s atmosphere from interplanetary 
space has been obtained by astronomical methods. 


Translated from Artificial Earth Satellites, 1960, no. 4, pp. 165 
170. Translated by Research Information Service, New York. 
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T. N. NAZAROVA 


Visual, photographic and radar tracking of meteors, ob- 
servations of the zodiacal light and of the Fraunhofer com- 
ponent of the solar corona, together with various methods 
of studying meteoric material reaching Earth’s surface 
have made it possible to formulate a physical theory of 
meteors and have been used to determine the space density 
of meteoric bodies in the solar system (1-7).' The first three 


! Numbers in parentheses indicate References at end of paper. 
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methods can be used to determine the number of particles 
entering Earth’s atmosphere, as well as their mass and 
distribution with respect to size. At present, however, these 
methods are only useful for registering particles with masses 
of 10-4 gm and more. 

Information on finer particles can be obtained from the 
integral characteristics of meteoric material, derived from 
observations of the zodiacal light, and also by resorting to 
extrapolation. 

As Table 1 shows, photometric investigations give highly 
divergent values for the density of dust in the vicinity of 
Earth’s orbit. 

Besides being of great interest to astronomers and geo- 
physicists, the study of meteoric material—its cosmic nature 
(numbers of particles, their space density and mass, etc.), 
distribution and the part it plays in the solar system—is 
also important in solving problems connected with the 
evolution and origin of planetary systems, since it helps to 
clarify a series of questions basic to the hypotheses of modern 
cosmogony. 

A work recently completed (8) deals with some new 
aspects of the influence of micrometeoritic particles on geo- 
physical processes in the atmosphere. 

An investigation of meteoric material is also necessary for 
the solution of certain problems of an applied nature, in 
particular, those connected with the study of the conditions 
under which rockets and artificial satellites move in cosmic 
space (9-11). 

As is known, the risk of a spaceship encountering a meteoric 
body capable of destroying it or damaging the cabin pres- 
surization system is small. The space density of meteoric 
bodies capable of doing serious harm to a rocket has been 
thoroughly investigated. On entering Earth’s atmosphere 
these bodies produce a glow, ionization, or penetrate to 
Earth’s surface in the form of meteorites. By analyzing 
the data from observations of the encounter between meteors 
and Earth’s atmosphere, it is possible to determine the space 
density, the mass of the meteoric bodies and a whole series 
of other important facts. 

The prolonged action of very fine meteoric particles on the 
skin of a rocket or satellite leads to the gradual erosion of the 
surface and may also cause damage to solar cells and optical 
instruments. Thus, the study of micrometeoric particles is 
a matter of great scientific and purely practical interest; 
moreover, it is of immediate concern, since, as previously 
mentioned, our only source of information at present about 
particles with masses less than 10~4 gm is in integral charac- 
teristics and extrapolation data. 

Evidently, data characterizing individual particles can 
only be obtained by direct methods, using rockets and, in 
particular, artificial satellites designed to spend long periods 
in the upper atmosphere. 

The impact of a meteoric body on an obstacle is attended 
by certain physical phenomena which can be utilized to design 
a recording instrument for mounting in rockets and satellites. 
Among the devices used in such instruments are special 


Table 1 


Space 
density p, 


Author gm per 
V. G. Fesenkov (1947) 
Allen (1947) 4.10728 
Van de Hulst (1947) 3-1072! 


Behr and Siedentopf (1953) 10-28 


Elsiisser (1954) 2-10-25 
Minaart (1955) 6-10~22 
Siedentopf (1955) 2-4-10-22 


crystal microphones, photomultipliers for registering the 
glow on impact, detectors consisting of a series of thin wires 
which break under the impact of a particle, detectors working 
on the principle of the accelerometer and piezo-electric de- 
tectors of various types. The erosion of the surface of a 
satellite as a result of the impact of micrometeoric particles 
is determined from the variation in the resistance of a thin 
metallic layer and also from changes in the intensity of beta- 
radiation (12-18). 

The third Soviet artificial satellite contained apparatus 
for recording the number of particle impacts and the momen- 
tum of the material of the detector ejected in the collision 
between the particle and its surface (19). There are various 
theories on the relationship between the “ejection impulse” 
and the mass, velocity, diameter, ete., of the particle. 

Theoretical calculations, made by Staniukovich (20), 
have shown that for high velocities the impulse recorded is 
proportional to the energy of the incident particle 


I ~ Eo 


where 
6 = coefficient depending on unsteadiness of the process 
and particularly on angular distribution of ejected 
mass (for a typical collision 6 ~ 4) 
€, = energy density of the crystal lattice of the vaporable 
body 
a = dimensionless coefficient whose value depends on 


the properties of the vaporable medium 


Thus, ~F (£)). 
According to the work of Lavrent’ev (21) 


I = A av! —0.3 


where 
A = constant 
a* = mass of the body 
ry = velocity of the body 
T) = minimum heat density necessary to convert substance 


to the gaseous state 


If we proceed from the simplest theoretical law relating 
the impulse, obtained with the detector, and the energy of the 
particle, while assuming that the mean velocity of the particle 
is equal to 40 km per sec, then the masses of the particles 
recorded can be determined from the relationship J = Al. 

In the period following the fifth IGY assembly, the pro- 
portionality factor A between the ejection impulse and the 
kinetic energy of the particle £y was established more ac- 
curately. Naturally, a change in the factor A involved a 
change in the computation of the mass of the recorded 
particle. The factor A has now been determined correct 
to a half order of magnitude. 

In the third Soviet satellite the “ejection impulses” were 
measured with the aid of a ballistic piezo-electric transducer 
(22). Electric signals, generated in the transducer as a 
result of collision with a meteoric particle, were separated 
with respect to amplitude into four ranges (channels) by 
means of an amplifier-converter, which counted the number 
of impulses in each channel. 

The apparatus was calibrated to record particles with 
masses within the ranges 8-10~° to 2.65-10-8, 2.65-10-8 to 
1.5-1077, 1.5-1077 to 5.6-10~*, and greater than 5.6-10~* 
gm. A signal was passed to the telemetering system after 
the accumulation of 32 impacts in the range of particles of 
smallest mass, after 16 and four impacts in the next two 
ranges, and after each impact in the last range. 

Vibration and noise from within the satellite was not picked 
up by the apparatus, since its thresholds were established 
beyond the level of this interference. 

It should be noted that the piezo-electric detector in the 
satellite, besides recording the particles impacting on its 
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-urface, also had a certain sensitivity to impacts on the body 
of the satellite. This was taken into account in processing 
the experimental data. 

Four piezo-electric detectors with a total area of 3410 ¢m? 
(including the housing) were mounted on the bottom of the 
satellite. 

If, as pointed out, we proceed from the theoretical law of 
dependence of the impulse, obtained from the detector, on 
the energy of the meteorie particle, and assume that the 
mean velocity of the particle is 40 km per sec, then the de- 
tectors recorded impacts from particles with masses of from 
S-10~° to 2.65-10-> gm, possessing energies of the order of 
104-10° ergs. 

During the telemetering periods of May 15 the impact 
frequency was from four to 11 per m? per sec; on May 16 
and 17 it fell to 4-107% impacts per m? per see and then be- 
came 5.3-1074 impacts per m? per sec. Then for eight days 
not a single count was registered; that is, the amplifier- 
converter did not accumulate enough impacts to transmit 
them to the telemetering svstem—the intensity of the flux 
of meteoric particles was less than 10~4 gm per m? per sec. 

Some of the features of the phenomenon observed on May 
15 should as noted: 

1 In spite of the considerable scatter (especially at 
apogee), it is possible to detect the existence of minima with 
a period of the order of 150 see (Fig. 1). The minimum 


Table 2 
Number of 
No. of impacts, 
loop Altitude, km 
2 1300-1500 5 
2 1700-1880 11 
3 400-600 7 
5 400-700 
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number of impacts is recorded at times when the satellite 
is turned with its bottom toward Earth. On these occasions 
it is mainly the cone acting, the base being only partially 
involved (because of its convexity). 

2 The impact frequency varies with time as the satellite 
continues in its orbit (Table 2). The location of the sections 
of the orbit in question is shown in Fig. 2. It is clear from 
the figure that the number of impacts varies with changes in 
the position of the satellite in its orbit but not with changes 
in altitude. 


Fig. 2 Location of sections of the satellite’s orbit over which 
the impacts of meteoric particles were recorded. a: Second 


loop; b: third loop; c: fifth loop 
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3 A large increase in the impact frequency was observed 
only on May 15. On the following days the impact fre- 
quency was three to four orders less. 

At the moment it is hard to say just what caused this in- 
crease in impact frequency. The phenomenon cannot be 
ascribed with certainty to the fact that, at the time, the 
satellite was cutting across a meteoric swarm [by reason of 
the excessive impact frequency (3)], although it is known 
from observational data that the space density of meteoric 
bodies in the Draconid shower is thousands of times greater 
than the density in the showers of active streams, where, in 
its turn, it fluctuates within the limits of two orders. Further 
investigation will help to clarify this point. 

In conclusion the author wishes to thank O. D. Komisarov, 
L. N. Neugodov, A. A. Trukhachev, L. Z. Rusakov, A. K. 
Bektabegov and G. M. Kurtev for their active participation 
in the work. 
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The crucial points involved in evaluating the information 
from this paper revolve around the calibration of the de- 
tectors. Granting that one may assign to the dust particles 
a mean velocity, and granting even that it is as high as 40 
km per sec, we cannot evaluate and compare the Russian 
data unless they explain exactly how they “established more 
accurately”’ the proportionality factor A between the ejec- 
tion impulse and the kinetic energy of the particle. This 
cannot be done in the laboratory as yet; therefore they should 
simply state the sensitivities of their detector instead. 

Equally obscure is the statement that “impacts on the 
body had a certain sensitivity, and this was taken in account”’ 
(italics mine). What was the “effective’’ detecting area? 

Nor are the results completely reported. One presumes 
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that the observed impacts were all recorded in the lowest 
amplifier channel; but what about the higher ones, what did 
they show? Did the observed impacts follow a Poisson 
distribution, as might be expected on the basis of random 
arrival of the particles? What is the statistical error in the 
reported rates? And what was the intensity before the re- 
ported increase of May 15? 

It is generally conceded that measurements of meteoric 
particles are among the most difficult and uncertain of 
satellite experiments. But this makes it all the more in- 
cumbent on the investigators to report their findings in the 
greatest detail possible. 


FS. SINGER 
Department of Physics 
University of Maryland 
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. Rotational Motion of Artificial 
58. 
SS, 
Bodies in Space x. ncnosm 
ter 
nd 
ite In this article the rotational motion of an artificial satellite whose center of mass follows a circular 
— orbit about the central planet is examined. For certain precisely stipulated assumptions, this 
problem has simple particular analytical solutions that may be of practical interest. In these 
op particular solutions, which we term “‘regular motions,”’ the center of mass of a body moves along a 
2. circular two-body orbit, and the body retains a constant orientation relative to that circular orbit. 
Furthermore, the stability of the regular motions relative to the Euler angles of the body’s axis of 
m symmetry and their first derivatives with respect to time are examined. ‘Two types of artificial bodies 
are specifically considered: A ‘*prolate body’’—a limiting case of which is a homogeneous rod, and 
R, an “oblate” body—a limiting case of which is a homogeneous flat disk. It is found that in the 
Pe, majority of cases the prolate and oblate bodies behave in an opposite manner with respect to the 
stability of the regular motions; that is, when a rod is stable, a disk is unstable, and vice versa. 
Finally, the effect of foreshortening of a central body on the regular motions of an artificial celestial 
body rotating around the central body is considered. 
st WE SHALL consider two absolutely rigid bodies VW and Coo — (Ao — Bo)p09 = OU/O¢o 
a Mo, the elementary particles of which mutually attract ; ; 
1 according to Newton’s law. Let m and A, B, C designate p = sing sindy + cos od 
mn mass and the principal moments of inertia of the body M; Po = SiN go SiN Doo + cos goto 
1e and mo and Ag, Bo, Co be similar quantities for the body Mo. ; ; 
oo Let us take a system of rectangular coordinates with an q = cos¢g Sindy — sin gd 
| origin at the center of inertia Go of the body and with 
ic fixed directional axes. If & 7, ¢ are coordinates of the center 
yf of inertia G of the body M, and y, 3, ¢, Yo, % and go are r=cosd~+e 
ie Euler angles of the bodies M and Mo, then the equations of : 
ce motion of the general problem of two bodies are written in the Tr = cos Doo + ¢o 
form (2)! ' f f f dmodm 
‘ _ The expression for the force function U is derived given a 
mo + modu sufficiently general specification of the shape and the con- 
7° an ae stitution of the bodies 1 and My. U may be presented in 
the form of a series, the first terms of which are written 
om m ou 
“or 
Ap — (B — C)qr = E cos 9 + cos A+ 3] 
Bq -(C — = -, E — cos Here r = VE + 7? + is the distance between Gy and G, 
and /, and / are the moments of inertia of the bodies My and M 
Cr — (A — B)pg = 0U/d¢ relative to the straight line GG. Hence (See Eq. 2) each of 
sin gy [OU these moments are dependent, aside from the coordinates 
AoPo — (Bo — Co)qoro = = | — cos Do ; ] a n, € of the point G, only upon the Euler angles of their own 
sin I LOW O¢o 
body. 
nits ou Let us assume now that the body M is very small in mass 
OW, and linear dimensions compared with the body Mo. A case 
: : such as this will be interesting to us, for example, when the 
|- body M is an artificial satellite and the body Mo is Farth or 
sin Jo Ovo the moon. 
; i Since each moment of inertia is a quantity having the form 
SIN fo 54 of the product of the mass of the body by the square of an 
= extreme linear dimension, the second term of the quantity 
Translated from Bulletin of the Institute of Theoretical Astron- [(mp + m)/mom]U has, for its first part, an expression of the 
omy, 1960, vol. 7, no. 7(90), pp. 511-520. Translated by 
Ernst Volgenau, Los Angeles. type f(mo + m)d?/r®, and:the.second part of this term is an 
1 Numbers in parentheses indicate References at end of paper. expression of the form f(mp) + ‘m)d?/r3. Consequently, this 
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second item is smal! enough to be disregarded in the case of 
an artificial satellite orbiting in the vicinity of the central 
planet. It follows that this term (and of course all following 
terms) may be neglected, and the equation of the motion of 
the point G will contain only the unknowns &, 9, ¢. We shall 
integrate the remaining unknowns from the equation of 
rotational motion. 

We shall consider the body Mo as a figure of revolution 
such that 4g = By. We shall assume that this body rotates 
uniformly around its axis of symmetry and that it will be 
possible to show this motion to be compatible with the 
equations of motion. As a consequence, the equations of the 
postulated motion will be written in the form 


ow 
g + rs 
. wn OW 
On 
ow 


where 


f(my + m) 


3%) 
r 


a = quantity having the order of flattening of the body 1/) 


Taking up next the equations of rotational motions, we will 
observe that in the derivatives of U with respect to the Euler 
angles, the first and second terms of U vanish, since they are 
not dependent on these angles; but the third term does de- 
pend upon the Euler angles, and appears in the expressions 
for the derivatives. Here also all the remaining terms of U’ 
are negligible. Likewise it is customary to consider M as a 
body of revolution with the moments of inertia A and C. 
At this time we shall introduce a distinction between two 
types of artificial bodies; that is, we shall call the body 
“prolate” if C< A, and “oblate”? whenC> A. Since B = A 
and 0U/d¢ = 0, then the third of the equations of rotational 
motion gives the integral cos dy + ¢ = w = constant, and 
the first two, after the elimination of p and q, can be exhibited 
in the form 


= —2cotan dW + cosec + 3fmoo 
r4 sin 


= sind cosdy* — csin + 3fmyo x 
[p cos sin (Y — v) — sind] 
where 
p= VE + 2! 


v = arctan n/t 
vy = (p/r) sind sin (¥y — v) + cosd 


o=1-C/A 
e = (C/A)w 


We observe that » is the cosine of the angle between the 
radius vector r to the point G and the positive direction 
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(which can be set up arbitrarily) of the axis of revolution of 
the body. Furthermore, the parameter o is positive for a 
prolate body, zero for the sphere, and negative for an oblate 
body. We also note that for the limit of the prolate body, 


for example, for the homogeneous rod, C = 0 and o = 1. 
For the limit of the oblate body, for example, for the homo- 
geneous circular disk, C = 2A and o = —1. 


The solution for the problem of rotational motion of an 
artificial body is carried out in the following manner: First, 
we integrate the equations of orbital motion, Nqs. 1. After 
that we determine the angle of precession y, and the angle 
of nutation d, from Eqs. 2. Finally, we obtain by integration 
the angle of spin of the body 


— go = wt cos ddt 


2. Integration of Eqs. | constitutes a classical problem, 
which we may consider solved and thoroughly investigated. 
Hence, our attention will be turned only to Eqs. 2 which, 
however, may not be integrated in terms of known functions 
and must, therefore, be solved by iterative methods. Con- 
sequently, partial solutions of the problems which correspond 
to motions possessing simple mechanical properties play a 
prominent role in the present investigation. Among such 
solutions we will mention here those having a great deal of 
practical significance, i.e., ‘‘regular’’ motions and also motions 
nearly regular with respect to the initial data. 

Inasmuch as the body Mo possesses symmetry relative to a 
plane passing through its center of mass and perpendicular 
to the axis of rotation, we shall always be able to select this 
plane as the fundamental coordinate plane &y. In this case, 
as is well known, Eqs. 1 admit of a particular solution corre- 
sponding to circular two-body motion of the point G in this 
plane, that is, in the equatorial plane of the central body 


r=p=a v= nt 


n? = (u/a*)(1 + 3a/a) 


In order to define the rotational motion of an artificial 
body, the center of mass of which moves on a circular orbit, 
we shall turn to Eqs. 2 which assume the form 


vy" = —2 cotan + (€/n) cosec + 
30 sin (W — 7) cos (W — 7) 


= sind cos 0p” — (e/n) sin dy’ + 
30 sin 3 cos sin*(y — 7) 


where in place of ¢t the quantity 7 = nt is utilized as the in- 
dependent variable so that the primes designate differentia- 
tion with respect to 7, and it is assumed that fmy = n?a'. 

Eqs. 3 define the rotational motion of the body M for 
any initial conditions except for the singular case for which 
sin & = 0. (Henceforth, we shall designate with Wo and 3 
the initial values of the angles ¥ and #, corresponding to the 
instant for which t = 0.) 

In order not to exclude from consideration the indicated 
special case, we shall transform Eqs. 3 with the help of the 
substitution 


x = —sind cos (vy — 7) y = sind sin (vy — 7) 


so that 
sin? = 2? + y? y —7 = — arctan y/r 
sind cos 0 = ax’ + yy’ 
— 1 = — ay’)/(@* + 
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‘The equations defining the new variables xz and y can accordingly be written in the form 


2- 
( nV 1 


€ 


— xy’)? — y” 

(ue! — xy’)? — — 


For these equations the point x = y = 0 is obviously an 
ordinary point; but, on the other hand, these equations 
present another peculiarity, i.e., at points for which x2? + 
/? = 1, which in the old svstem of variables corresponds to the 
case of cos 3 = 0. 

In such a manner Eqs. 3 and 4 encompass all the possible 
cases of rotational motion of the body, whose center of mass 
describes a circular two-body orbit lying in the equatorial 
plane of the central planet. 

As is well known, Eqs. 1 do not have other solutions for 
which the circular orbit lies in planes other than the equatorial. 
However, if a is very small, Eqs. 1 admit of solutions slightly 
different from two-body circular solutions; such solutions 
are also of considerable interest. 

In order to obtain such solutions we shall seek integrals to 
Iqs. 2 in the form of the series 
a4... 
expanded in powers of the parameter a. It stands to reason 
that the functions satisfying Eqs. 1 are also presented in a 
form similar to the foregoing series, as is, at any rate, well 
known from classical theory. To obtain the zeroth approx- 


imation, i.e., the functions ¥ and @, it is necessary to assume 
that a = 0 in Eqs. 1 and 2. After this assumption, Eqs. 
1 will become those of a two-body problem and will, accord- 
ingly, admit of a solution representing a circular orbit in an 
arbitrary plane through the origin of the coordinates, i.e., 
through the center of inertia Go of the body Mo. It is now 
permissible to take the plane of this circular two-body orbit 
as the plane &; then, as easily seen, the equations defining 
the functions ¥ and @ will coincide with Eqs. 3. In a number 
of cases, it will be possible to satisfy the zeroth approximation 
at least in regard to the representation of the initial character 
of the motion. But it also will not be difficult to consider 
the influence of the terms of first degree in a. The discussion 
concerning this influence will be presented in the following. 

Let us now consider regular and nearly regular motions, 
ic., examine the simplest particular solutions of Eqs. 3 
or 4. 


3. It is not difficult to prove that Eqs. 3 have the par- 
ticular solutions 


=H 0 cos = €/n 

In order to establish the nature of the motion corresponding 
to this solution, let us observe that the quantity v and the 
cosine of the “angle of attack,’ ie., the angle between 
the direction of velocity of the point G and the direction of 
the axis of rotation of the body M [see (2)] for the cir- 
cular orbit of the point G are defined by the formulas 


vy = sind sin (y — 7) cos y = — sind cos (W — 7) 
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Hence in the particular solutions being considered, which 
we have previously called “the pointer’ solution (3), cos y 
= + sin & and vy = 0. In this regular motion the axis 
symmetry of the body M remains perpendicular to the 
radius of the circular orbit of the point @ (i.e., the body 
maintains a horizontal attitude) and will form with the ve- 
locity direction of point G the angle y equal to either 7/2 — 
d or 34/2 — %. Also, it is not difficult to see that for this 
motion the angular velocity of the naturally rotating body 
remains constant. 

Let us now consider the question of the stability, for the 
regular motion of the body with horizontal attitude relative 
to the quantities y, y’, J and 3’. For this purpose let us 
assume 


and form the equations of perturbed motion, i.e., the motion 
slightly differing, according to the initial conditions, from the 
regular motion of the body in horizontal attitude. These 
equations are presented in the form 


x", = 30x, — y’; cotan J) + X, 


= sin? + 2x’; sin cos + Yi 


where X, and Y, designate the terms above the first order in 
the expansion of the equations of the perturbed motion with 
respect to the powers of 2, 21, y1, yi. As is well known, the 
solutions of the problem of the stability of the horizontal 
attitude problem depend, first of all, on the nature of the 
roots of the characteristic equation, which may be written in 
the form 


x? — 30 


x cotan Jo 
. => 0 
— x sin J cos Vo 


x? + sin? I 


or in expanded form 


(1 — — sin? = ( 


From this 


It is easy to see that if the body W is elongated, i.e., if o> 0, 
then for any value of sin J # 0 both values of x? are real, 
where one is positive and the other is negative. Hence the 
characteristic equation has obviously a positive, a negative 
and two complex conjugate roots. From here it follows that 
the regular motion of a body having a horizontal attitude is 
unstable for ¢ > 0, but may possess conditional stability 
depending on one arbitrary parameter. If, however, the 
body M is foreshortened, i.e., if ¢< 0, then both the values of 
x” are real and are negative, for any value as before of sin 3 
~ 0. Hence the characteristic equation has in this case a 
pair of complex conjugate roots; moreover, it has two imagi- 
nary roots. It follows that the regular motion of a body with 
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a horizontal attitude for ¢ < 0 is stable at least in the first 
approximation, so that librations will result. It is impossible 
to decide whether this motion will actually be stable or un- 
stable on the basis of considerations of only the first approxi- 
mation. 

The angular velocity of proper rotation ¢ in the perturbed 
motion will be defined by the formula 


¢ = w — n(1 + cos (Io + m1) = 
Go — nx’; cos + ny SiN Do + ... 


From here it follows that for o > 0 the angular velocity 
increases without restriction secularly with t, and for « < 0 
this angular velocity, at least in the first approximation, 
remains bounded. 

It will now be possible to define the motion near to regular 
motion at a horizontal attitude. According to the mathe- 
matical theory of Liapunov (4), it is possible to find the 
solution of the system Eqs. 5 in the form of a series expanded 
according to the degree of the perturbation obtained, i.c., 
according to the power of the quantities 210, x'10, yio and y's. 
The series is absolutely convergent for significant values of 
T not exceeding known limits defined in terms of the maxima 
of the quantities |210|, |x’10|, |y1o| and {y’10|. The first terms of 
these series are obtained by integration of the equation of 
first approximation, i.e., Eqs. 5, in which the terms of the 
higher order in X, and Yj are discarded. For illustration, let 
us use the initial terms of such series in two special cases—for 
o = +1. Moreover for simplication of the results, we shall 
restrict ourselves to the case in which cos 3 = 0 so that 
¢@ = 0. Then we shall have: For ¢ = 1 


= cosh V3r sinh V3 7 
yocost + sint +... 
foro = —1 
= cos V3.7 + 3) sin V3 7+... 


Yi = yocost + y'osint +... 


4. Let us turn to the following particular solution which 
we shall refer to as the vertical attitude solution. In this 
solution 


3m 
4 
y=r't+ y’ =1 =H 
2 
cos Bo = €/n(3o0 + 1) 


Hence cos y = Oand »y = ¥ sin dp; i.e., in this regular motion 
the axis of symmetry of the body M remains perpendicular 
to the direction of the velocity of point G and forms with the 
radius vector of this point an angle equal to either 37/2 — 
or — vd. The angular velocity of the naturally 
rotating ¢ in this motion also remains constant. 

For the purpose of examining the questions regarding the 
stability relative to the quantities y, y’, 3, 3’, for the case 
of the motion with vertical attitude let 


y=rt +2. 


2 
us 
2 


The equations of the perturbed motion can be presented in 
the form 

x", = — + (30 — 1) cotan + Xe 
y"2 = — (1 + 3¢) sin? doy. + (1 — 3a) cos Dox’, + 


where X_ and Y2 indicate a set of higher order terms. The 
characteristic equation has the form 


x? + 30 (1 — 3) cotan 


—(1 — 3¢) cosdy:x x? + (1+ 30) sin? tied 


or 
pr? +q = 0 
where the quantities p and q denote 
p = 380 + (1 + 3a) sin?d + (1 — 30)? cos? = pi + pot ps 
q = 3o0(1 + 3a) sin? Jo = pipe 


For the elongated body, i.e., in the case a> 0, it is obvious 
that p > 0 and qg> 0. Forming the discriminant of the 
quadratic equation, we obtain 


44 = + po + ps)? — 4pipo = (pi — pe + ps)? + 4pops> 0 


Hence both values of x* are real and negative, and it follows 
that for ¢ > 0 the regular motion with vertical vehicle atti- 
tude is stable in a first approximation, and librations will 
occur. 

For the flattened body, i.e., in the case ¢ < 0, the question 
cannot be so easily settled and additional examination is re- 
quired. For o <Oand 1+ 0 it follows that < 0 and 
4A = p? — 4q> 0. This indicates that both values of x? 
will be real so that one is positive and the other negative. 
Hence the characteristic equation will have one positive 
root. It follows that in this case the regular motion with 
the vertical attitude will be unstable. 

Now let ¢ < 0 and 1 + 30 < 0. Then qg> 0 and the 
nature of the motion depends on the sign of the discriminant 
which we write in the form 


4A = b?s? + 2(ab + 2c)s + a® 


where 
s = sin? Bo 
a= 1- 30+ 9o0?> 0 
b = + 0 
c = 380(1 + 30)> 0 


Since for s = 0 we have A > 0, then for 
0 < sin? d < = (1/b%)[ab + 2c — V(ab + 2c)? — a%®?] 


both values of x will be real and negative indicating that for 
this condition the regular motion of a satellite with vertical 
attitude will be stable in a first approximation. If, however 


sin? > 5 


then A < 0 and both values of x? will be complex conjugates; 
hence the characteristic equation will, in this case, have roots 
with positive real parts. The instability of regular motion 
of the vertical attitude vehicle for this condition is here evi- 
dent. With this, the question of stability of regular motion 
of the vertical attitude vehicle is completely settled. 

The angular velocity of the proper rotation for the per- 
turbed motion is defined by the formula 


= go — Vo + nyo sin + 


Again it follows that in the case of instability in the vertical 
attitude motion, ¢ increases without restriction, and in the 
case of stability it remains (at least to the first approximation 
indefinitely constant. As discussed previously, we shall 
determine the initial terms of the series presenting the per- 
turbed motion in the two limiting cases of ¢ = +1 and for 
the conditions that cos 3% = 0. We have: For = 1 


Zo = COS V3 + (2x'o sin V3 
Yo = Yo cos 27 + (1/2) yo sin 27+... 
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and for ¢ = —1 


= cosh V3 + sinh V3 7+ ... 


Y2 = cosh V2r + 2) sinh V2 7+... 

5. In conclusion, we will consider recent particular solu- 
tions corresponding to a regular motion which we call ‘the 
buoy.” In this particular solution sin % = 0, and hence 
qs. 8 are unsuitable. Let us take Eqs. 4, which have the 
trivial solution z = 0, y = 0 corresponding to regular motion 
n which the axis of symmetry of the body remains perpendic- 
iar to the plane of the circular orbit of the point G. Also, 
it is obvious that the angular velocity ¢ of the natural rotation 
remains constant. The right-hand members of Eqs. 4 
are arranged by powers of x, y, 2’, y’. We obtain the equa- 
tions of perturbed motion in the form 


= 4x 
n n 


where X and Y designate terms of higher order. The charac- 
teristic equation is written 


(1 + 30 - £) 


where 


p=1- 380+ (1 — €/n)? 
q = (1 — &/n)? + 3a0(1 — 


Hence if the parameters o and ¢ are such that they fulfil! 
simultaneously the inequalities 


p> 90 q> 9 


both values x? will be real and negative, and the characteristic 
equation has two pairs of imaginary roots. Hence the regu- 
lar motion of “the buoy”’ will be stable in the first approxima- 
tion. If, however, only one of the conditions of Eq. 8 fails to 
be fulfilled, the regular motion of “the buoy” will be unstable. 
For example, in the limiting cases in which o = +1 and e = 
2n, there will be an instability for the prolate body and a 
stability for the oblate body. Perturbed motions are de- 
fined in this case by the series: For o = 1, € = 2n 


p? — 0 [8] 


+... 


= Yo cosh V2r + (y'o/V 2) sinh V2 7 +... 


+. 
Y = yo cos 27 + (1/2)y’o sin 27 + . 


It should be noted also that Eqs. 4 are useful for the in- 
vestigation of regular motion of vehicles with vertical 
and horizontal attitudes providing cos }) # 0. Actually, it is 
easy to see that Eqs. 4 are satisfied for 

c= + sind, y= 0 cos = €/n 
which corresponds to the motion with horizontal attitude and 
y = F cos Jy = €/n(30 + 1) 


which corresponds to the motion of “the buoy’ type. 
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Hence with the aid of Eqs. 4 it is also possible to consider the 
case of the stability of these regular motions relative to the 
quantities z, y, x’ and y’. 


6. In conclusion we will consider briefly the question of 
the influence of a foreshortening of the central body M, on the 
rotational motion of the body M. As previously mentioned, 
we shall display the solution of Eqs. 2 in the form of the series 


= 


and we shall define the functions ¥, 3, Yi, d1.. .by the initial 
conditions 


YO) = vo ¥(0) = % 50) =% 80) = 
=0 ¥,(0) =0 =0 = 0 
(k =1,...) 


Then ¥ and @ are defined by means of Eqs. 3; for the 
definition of the functions y, and J; it is necessary to form 
the corresponding equations. In order to obtain a precise 
idea of the effect of the size of the foreshortened body My 
on the rotational motion of an artificial body, we will consider 
only that case in which the unperturbed equations have the 
simplest possible form. 

For this purpose we shall limit ourselves to the case for 
which @ = 0 and the point G@ describes a circular Keplerian 
orbit. The rotational motion of the body M for a = 0 
corresponds to a regular motion of the vehicle with a vertical 
or horizontal attitude. Then the equations defining the first 
approximation (with which we will be satisfied) are: For 
the horizontal attitude 

for the vertical attitude 
vy", = — 

From here it can immediately be seen that in both cases 
vd; = 0; ie., the angle of nutation does not experience per- 
turbations as a result of the foreshortening of the body Mo. 
In order to define y, it is necessary to take from the classical 
theory v;, which, as is known (5), is the sum of secular terms 
and an infinite series of periodic terms. 

For our purpose it will be sufficient to consider, in the 
expression for v;, only the infinite terms, i.e., to put u = m7, 
where n; is a certain constant. Then the equation defining 
the function y, can be written in the form 


vi = — mr) 


and its solution satisfying the initial conditions can be ob- 
tained by means of the formula 


= [rt — (1/k) sinh kr] = +30 


when +3¢ is positive, and by means of the formula 
vy = mir — (1/k) sin kr] k? = + 30 


when +3¢ is negative. In both cases the foreshortening of 
the body causes infinite perturbations to the angle of pre- 
cession y. 


—Received June 27, 1959 
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Reviewer’s Comments 


The results obtained in this paper are developed from 
fundamentals; thus many well-known results are rederived. 
In particular, the approximation for the gravitational inter- 
action of two asymmetrical bodies (which the author actually 
refers to in one of his previous publications on the subject) 
is part of the classical treatment of the nutation and libration 
of celestial bodies. The transformation of coordinates used 
to establish Eq. 4 is a standard procedure used to eliminate 
similarities arising in the treatment of a motion of heavy 
gyroscopes for which the axis becomes vertical. The stability 
of the rotation of space vehicles is obtained by Routh’s 


method. 


The reviewer feels that the results established by the 
author could have been obtained more rapidly by employing 
a Lagrangian rather than a Newtonian formulation. How- 
ever, the paper is very well written and systematically 
organized. Indeed, the inclusion of fundamental details 
which would have been eliminated for reasons of space by a 
Western editor, except in a review paper, does have the 
desirable effect of making a paper of this kind readable with- 
out an undue investment of time. 


—EuioT T. BENEDIKT 
Astro Sciences 
Northrop Corp 


Density of Liguid Ozone 


HE VARIOUS values given (1 to 3)! for the density of 

ozone differ by 8%; they are, in gm per cm’, 1.752 (1), 
1.614 + 0.004 (2), and 1.613 + 0.0004 (3). We have meas- 
ured the density of ozone at — 195.6 C and find it to be 1.619 
+ 0.004. The gaseous ozone (100%) was admitted into a 
known volume and its pressure was measured; then it was 
condensed at — 195.6 C into a calibrated tube 2 mm in diame- 
ter (one division was equivalent to 0.01 em’). The error in the 
volume measurement was not more than +0.002 %. 


Translated from Vestnik Moskovskogo Universiteta (Bulletin of 
Moscow University), 1960, no. 4, p. 11. Translated by Research 
Information Service, New York. 

! Numbers in parentheses indicate References at end of paper. 


Technical Comment 
D. T. SWIFT-HOOK! 


Central Electricity Generating Board, London, England 


T WAS suggested in a paper recently published in this 
Journal (1)? “that a sphere passing the observer with the 
relativistic speed appears noticeably flattened in the direction 

of motion.” 
I should like to draw attention to a paper by Penrose (2) in 


Received June 29, 1961. 

‘Group Leader, Magnetohydrodynamic Generation Group, 
Research and Development Dept. 

? Numbers in parentheses indicate References at end of paper. 


T. I. EMEL’IANOVA, 
B. V. STRAKHOV and 
V. P. LEBEDEV 


Department of Physical Chemistry 
Moscow University 


The method was tested with liquid oxygen, whose density at 
— 195.6 deg was found to be 1.204 gm per em’, which agrees 
well with the currently accepted (4) value of 1.2037 gm per 
cm’. Measurements on ozone yielded values of 1.620, 1.621, 
1.619 and 1.617 (mean 1.619 + 0.004), which agree well with 
values published recently (3). 


—Received March 15, 1960 
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which it is shown that the appearance of a sphere, no matter 
how it is moving, is always such as to present a circular out- 
line to any observer. A discussion at a more elementary 
level has been given by Weisskopf (3), who points out that, 
although an object is described as contracted in a coordinate 
system with respect to which it is moving, it is not seen as con- 
tracted by an observer at rest in the coordinate system. 
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USSR, GEOPHYSICS SERIES (Izvestiia Akademii 
Vauk SSSR, Seriia Geofizicheskaia). Published 
by American Geophysical Union, Washington, 


D.C. 
No. 4, April 1959. 


Stability of Planetary Scale Atmospheric Motion, M. B. (alin, 
pp. 393-400. 

Abstract: A two-layer spherical model is used to investigate the 
stability of atmospheric motion. It is considered how a difference 
between the circulation indexes for the different levels affects the 
stability conditions and characteristics of the motion. It is shown 
how the theory can be used to predict a pressure field. 


No. 7, July 1959. 


P-Waves Diffracted at Earth’s Core and Rigidity of the Core, 
IL. N. Rykunov, pp. 678-684. 

Abstract: We investigate the possibility of evaluating the shear 
modulus of Earth’s core on the basis of studying the character 
of decrease of amplitudes of P-waves in the shadow zone. In 
view of the considerable difficulties which arise in a theoretical 
study of the diffraction of seismic waves at Earth’s core, the 
investigation will be carried but by the method of models. 


Determination of Types of Reflected Multiple Waves by Their 
Kinematic and Dynamic Characteristics, A. M. Epinat’eva and N. 
G. Mikhailova, pp. 685-695. 

Abstract: This paper presents the experimental data on the 
kinematic characteristics of multiple reflected waves recorded 
with a medium frequency (f = 37 eps) and high frequency (f = 
105 eps) instruments. These instruments were used to determine 
types of multiple waves. By using the wave amplitude, it is 
possible to define more accurately the wave propagation with 
multiple reflections. 

Conclusions: Based on experimental data obtained in the region 
of intensive multiple reflected waves, the latter were classified as 
multiple and simple waves. 

The presence of various types of multiple waves, fully and 
partially multiple, was established. 

A predominant number of recorded multiple waves are reflected 
at least once from the same interface. This boundary occurs at a 
depth of H ~ 800 m and is characterized by a high reflection co- 
efficient (q ~ 0.3). 

It was established that for the fully multiple waves, which are 
reflected from the lower boundary at a depth of H ~ 800m, the 
upper reflecting interface occurs near the base of the low velocity 
layer. 

Using vertical travel time curves, it is possible to determine 
possible paths of multiple waves which undergo reflections at 
different interfaces. It has been shown that kinematic charac- 
teristics of the waves in a series of cases are not sufficient for the 
determination of the propagation path of a multiple wave. 

A method is suggested which allows a refined determination of 
the position of the reflecting boundaries for multiple waves based 
on the use of amplitude ratios of multiple and simple waves. 

The types of multiple waves which are recorded with seismic 
instruments tuned for recording different frequencies (MF—37 
cps, HF—105 eps) are different. With the experimental data, it 
has been clarified that the number of simple reflected waves which 
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are recorded with HF instruments is considerably higher than 
with MF instruments in the same time interval. 

The difference in wave types recorded with HF and MF 
instruments probably is explained by unequal absorption of 
oscillations of different frequencies and with different absorbing 
properties of the medium at different depths. 


Feasibility of Using Shielded Detectors in Airborne Radiometric 
Prospecting, V. A. Vorob’ev, pp. 705-709. 

Abstract: The author analyzes the problem of raising the 
effectiveness of recording local aero-gamma anomalies with deeply 
shielded detectors. The separation factors for ray approximation 
are computed and the relations for selecting the optimum shield 
aperture (@pt) and the radiometer time constant (7op:) are ob- 
tained. It is shown here, that the use of round, conical shields is 
feasible if In, < 0.3 Ip & (uh), where Jy ® (uh) is the intensity of 
Karth’s gamma field at flight aititudes, and J, is the residual 
background. The influence of the shape of the shield on the 
quality of the gamma survey is discussed. 


Determination of Turbulent Flow of the Heat, A. V. Perepelkina, 
pp. 727-733. 

Abstract: Results of the test of the instrumentation con- 
structed at the Institute of Physics of the Atmosphere of the 
Academy of Sciences of the USSR, for direct measurements of 
turbulent flow of the heat by pulsation methods, are given. The 
instrumentation permits a direct determination of the turbulent 
flow of the heat due to the lineality of instruments (acoustic micro- 
anemometer and a pulsation resistance thermometer) and the 
application of a correlometer. Pulsation measurements were 
accompanied by the gradient measurements. With the results of 
measurements, dimensionless empirical functions are derived, 
which enter into the formulas for the determination of turbulent 
flow of the heat with the data of the gradient observations. The 
case of instability of stratification is investigated and the value of 
the universal constant characterizing the conditions of free 
convection close to the value determined by Priestley is obtained. 

Conclusions: Analysis of the data obtained in testing the new 
instrumentation for direct determination of the turbulent heat 
flow, permits an empirical development of the first approximation 
of the dimensionless function of Richardson’s number, included 
in the formula for the determination of q with the data of the 
gradient measurements. 

Especially, the case of unstable atmospheric stratification was 
investigated and the value of constant Fre for the free convection, 
close to the value proposed by Priestley, was obtained. 

In relation to the transition conditions there was found some 
deviation from the results of observations in Australia. It can be 
hoped that further observations will clarify this problem. Due 
to the lack of information it was impossible to analyze the case of 
a stable stratification. 


Some Results of Experiments Conducted With a Nuclear 
Magnetometer, A. Ya. Rotstein and V. 8. Zierel, pp. 763-765. 


No. 8, Aug. 1959. 


1. Propagation of Elastic Waves in a Layered Medium, («. ». 
Pod’’yapol’ski, pp. 788-793. 

Abstract: The paper analyzes the application to a multilayered 
medium of functional-invariant solutions of the dynamic prob- 
lems of the theory of elasticity as suggested by Smirnov and So- 
bolev, as well as of the asymptotic evaluations of the wave fields 
in the neighborhood of the fronts of reflected and head waves, ob- 
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tained by Zvolinsky on the basis of this method. The analysis is 
limited only to those phenomena which are related to the charac- 
teristics of the wave field of the unit elementary wave. The phe- 
nomena which are basically related to the interference of the ele- 
mentary waves are not analyzed in this paper. 


Double Refraction in Earth’s Mantle, A. V. Vvedenskaya and L. 
M. Balakina, pp. 814-819. 

Abstract: Observations of the field of displacement of longi- 
tudinal (P) and transverse (SV, SH) waves propagating in the 
mantle of Earth, disclose an increase in the amplitudes of P and 
of SV waves in relation to SH waves during the penetration of 
seismic rays to the depths of: 250-500, 900-1000, 1200-1300, 
1800 and about 2200 km. These peculiarities of the displace- 
ment field may be connected with a polarization of the transverse 
waves during double refraction in the anisotropic layers of Earth’s 
mantle which correspond to the depths enumerated. 


Penetration of Fast Charged Particles From the Outer 
Atmosphere Into the Ionosphere, V. D. Pletnev, pp. 833-834. 

Abstract: The oscillatory motion of charged particles trapped 
in Earth’s magnetic field is considered. It is shown that in such 
streams of oscillating particles the directions of motion will 
exhibit pronounced spatial anisotropy because of absorption in 
the lower atmosphere. 


Ultraviolet Solar Spectrum in the 2636-2937 A Region, V. P. 
Kachalov, N. A. Pavlenko and A. V. Yakovleva, p. 844. 

Abstract: Absorption lines of the solar ultraviolet spectrum in 
the 2636-2937 A sector were photographed with a concave diffrac- 
tion grating, measured and identified. 


Acoustic and Gravitational Oscillations in the Atmosphere, |.. A 
Dikii, pp. 849-853. 

Abstract: Wave motion is investigated in the atmosphere with 
vertical temperature stiatification (a polytropic atmosphere). It 
is shown that this type of motion can be divided naturally into 
“acoustic-type”’ waves and “gravitational-ty pe”’ waves, similar to 
those observed in an isothermal atmosphere. The frequency 
spectra of these waves are also investigated 


Observations of Short-Periodic Variations in Earth’s Electro- 
magnetic Field, B. E. Brunelli, M. N. Berdichevsky, A. M. 
Alekseev and O. A. Burdo, pp. 864-871. 


Fig. 1. General Circuit of Measuring Device 


I - oscillograph; II — magnetometer 


Introduction: The present article is devoted to a description of 
observations conducted from September 4-14, 1957, in one of the 
southern regions of the Tumen’ region (latitude 56°40’, longitude 
67°). The observations were conducted in one point (which we 
shall henceforth term the magneto-telluric station) upon the 
initiative of Zagarmistr for verifying the correctness of the 
hypotheses of the magneto-telluric profiling proposed by Berdi- 
chevsky. 


JOURNAL OF APPLIED MATHEMATICS AND 
MECHANICS (Prikladnaia Matematika i Mek- 
hanika). Published by Pergamon Institute, New 
York, in conjunction with American Society of 
Mechanical Engineers. 


Vol. 24, no. 2, 1960. 


Three-Dimensional Hypersonic Gas Flow Past Slender Bodies 
at High Angles of Attack, V. V. Sychev, pp. 296-306. 


Introduction: We consider three-dimensional hypersonic flow 
past bodies whose transverse dimensions are substantially smaller 
than their length. Making use of the smallness of a parameter 
characterizing the relative thickness of the body, it is possible to 
put the problem of flow past such a body approximately in a form 
that generalizes the similitude of hypersonic small-disturbance 
flow to the case of arbitrary angle of attack. 

Approximate formulas are obtained for the calculation of the 
aerodynamic characteristics of slender bodies at high angles of 
attack, containing as unknowns only certain constants depending 
on the cross-sectional form of the body. 


Approximate Computation of Optimal Control by the Direct 
Method, N. N. Krasovskii, pp. 390-397. 

Abstract: A description is given of an approximate method for 
computing the optimal (in regard to speed) trajectory in a linear 
control system. The optimal problem is replaced by an auxiliary 
“smoothed out’”’ problem which is treated by the usual methods of 
variational calculus. 


A Limiting Process in the Solution of an Optimal Control Prob- 
lem, F. M. Kirillova, pp. 398-405. 

Abstract: An optimal problem is considered for a linear system 
of differential equations in which the integrals of the pth power 
(p> 1) of the modulus of the control function are bounded. 

Text: We shall assume that the control function u(t) satisfies 
the condition 


|pdrS (p> 1) [1.2] 

The optimal problem consists of the following: Among all con- 
trol functions satisfying condition 1.2, it is required to find such a 
u(r, p) that a point moving along the trajectory of Eq. 1.1 will 
move from the initial point z(t) = x to the origin of the co- 
ordinate system in the shortest time 7'(p). 


Two Classes of Plane Extremal Motions of a Rocket in Vacuum, 
Iu. A. Gorelov, pp. 434-442. 

Introduction: A study of the conditions is given, the fulfilment 
of which will secure a motion of rockets along curvilinear trajec- 
tories which is extremal in time and expenditure of mass. 

The motion is considered for an active segment of a trajectory, 
i.e., the mass of the rocket is assumed to vary in time. 

Within the framework of the indicated hypotheses, the varia- 
tional problem is formulated, regarding the determination of the 
characteristics of extreme maneuvers [turns, rotations] in time 
of rockets by given angles, for given initial and final velocities of 
motion, initial and final weight of rockets. The solution of the 
problem is sought for plane motions in horizontal and vertical 
planes separately, under the assumption of absence of the 
influence of aerodynamic forces. 


JOURNAL OF PHYSICAL CHEMISTRY (Zhurnal 
Fizicheskoi Khimii). Published by The Chemical 
Society, London. 


Vol. 34, no. 3, March 1960. 


Thermodynamic Theory of Critical Phenomena in Ternary 
Systems. I. Thermodynamic Inequalities Characteristic of the 
Critical State, A. V. Storonkin and A. I. Rusanov, pp. 251-254. 

Summary: 

1 The principal thermodynamic inequalities expressing the 
conditions for the stability of the critical phase and defining the 
type of critical point have been considered. 

2 The concepts “type of critical point with respect to tem- 
perature and pressure”’ and “type of critical point with respect to 
the concentration of one of the components’’ have been discussed. 


The Kinetics and Mechanism of the Oxidation of Methane. 
I. Fundamental Macrokinetic Laws, L. V. Karmilova, N. S. 
Enikolopyan and A. B. Nalbandyan, pp. 261-265. 

Introduction: This paper describes a detailed study of the 
kinetics and composition of the products of oxidation of methane 
throughout its course in a reaction vessel whose surface is of a 
strictly specified nature. The laws for the accumulation of 
intermediate products (CH.O and H:02), some new phenomena 
observed during the oxidation, and the probable mechanism of 
the reaction will be considered in subsequent papers in this series. 

Summary: 

1 The kinetics of the consumption of the original reactants 
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(CH, and O,) and the accumulation of the reaction products (CO, 
CO2, He, HO) when methane is oxidized to obtain a large measure 
of conversion, were studied in a quartz vessel with a stable surface 
under a wide range of conditions (7' = 423°-513°, Py) = 117-375 
mm Hg, CH,:02 = 0.5, 1, 2). 

2 It was shown that the maximum rate at which the original 
reactants are converted and the final products accumulated 
corresponds with the stoichiometric proportions WcH, = Wco. 
Wo, */2WecH, Wu.0, = 2Won,, and WcH, 2AP dt. 

3 The effective energy of activation for the oxidation of 
methane was found to be Kefg = 42.5 + 1 keal per mole, inde- 
pendent of the way in which it is derived [from the temperature 
variation of the pressure, AP, or (WCH,)max or (Wo:)max], and it 
agrees with the effective energy of formation of CO and H.0. 

4 The overall order of the reaction (n = 2.7) and the order 
with respect to methane (a = 1.62) and oxygen (8 = 0.96), and 
the total pressure y & 0.1 were determined. 

5 The temperature coefficient (#7; = 36 + 1 keal per mole) of 
the induction period, as well as its variation with the pressure of 
the initial mixture (m = 1), were determined. 


Distorting Effect of Layer Adjacent to Electrode on Measure- 
ments of Dielectric Constants at Sonic Frequencies, A. N. 
Efremov, pp. 268-271. 

Summary: 

1 The method described enables correct values of dielectric 
parameters, free from the distorting effect of the layer adjacent 
to the electrode, to be obtained at low frequencies. 

2 Both exact and approximate formulas are given for calculat- 
ing the dielectric parameters of a substance and of its electrode 
laver. 

3 The practical applicability of the method has been shown 
with water, alcohol and acetone as examples, and it has been 
found that the dispersion, observable by the usual method of 
measurement with these substances, does not in fact exist in the 
sonic range right down to a measured frequency of 40 eps. 

4 The anomalous dispersion observable in the sonic range in 
the dielectrics studied results from the occurrence of dispersion in 
the electrode layer. 

5 The proposed method yields a comparatively accurate 
value of the dielectric constant without the necessity for a high 
degree of purification of the dielectric. 

6 It is possible to investigate the properties of the layer ad- 
jacent to the electrode by the dielectric method. 


Effect of Initial Temperature on Critical Diameter of Nitro- 
glycerin and Trinitrotoluene, A. F. Belyaev and R. Kh. Kur- 
bangalina, pp. 285-289. 

Introduction: The critical (limiting) charge diameter of an 
explosive is taken to be the minimum diameter at which steady 
propagation of detonation is still possible under given conditions: 
if the diameter of the charge is less than the critical value, detona- 
tion is extinguished. The physical significance of the critical 
diameter was demonstrated by Khariton in his theory of detona- 
tion capacity. According to him, the critical diameter is propor- 
tional to the time of reaction in the detonation wave front. A 
short time of chemical reaction facilitates the occurrence of 
detonation and makes an existing detonation more persistent. 

We have carried out experiments on the effect of the initial 
temperature on the critical diameter of nitroglycerin and TNT; 
change in critical diameter should indicate change in the time of 
chemical reaction in the detonation wave front. There is a sug- 
gestion in the literature of a considerable change in the critical 
diameter of liquid nitromethane with change in temperature. 
Nitromethane, however, is not a typical explosive with well 
established kinetic constants, so these data cannot be regarded as 
in any way adequate. 

Nitroglycerin and TNT were chosen as typical of the most 
thoroughly investigated explosives, and their behavior in this 
respect is of considerable interest. The work was carried out 
mainly with liquid nitroglycerin and TNT; some exploratory 
experiments were carried out with TNT in the form of a low 
density powder. 

Summary: 

1 The critical diameter falls rapidly with increase in tempera- 
ture for both nitroglycerin and liquid trinitrotoluene. 

2 Near its melting point TNT has a critical diameter which is 
10-100 times as great as that of nitroglycerin at room tempera- 
ture. In conformity with this, its explosive properties (detona- 
tion capacity, sensitivity) differ very greatly from those of nitro- 
glycerin. Near the flash point the critical diameter of TNT is 
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only about three times as great as that of nitroglycerin, and conse- 
quently its explosive properties should approach those of nitro- 
glycerin. 

3 If it is assumed that liquid TNT and nitroglycerin detonate 
because of the heating effect due to compression with subsequent 
homogeneous reaction, the temperature produced by the shock 
wave (initial temperature of the chemical reaction) should be of 
the order of 1100-1200 K. 

4 The mechanism of the chemical change occurring with TNT 
powder is apparently different from that for liquid TNT. 


Relaxation Mechanism of Propagation of Combustion in Hetero- 
geneous Exothermic System, Z. I. Fur, pp. 289-293. 

Introduction: Kinetic theories of chemical reactions which occur 
in the propagation of combustion are inadequate when the com- 
bustion of certain solid systems, such as smoky gun powders, 
thermites, etc. is considered, since in these cases, theories that sug- 
gest that the propagation of combustion depends on the tem- 
perature and dispersion of the system cannot be used. 

In this paper an attempt is made to link the rate at which com- 
bustion is propagated with the rate at which heat is transferred 
into adjacent layers of the solid phase. This type of theory was 
first applied, without success, to the burning of gas mixtures by 
Mallard and Le Chatelier. 


Magnetochemistry of Active Centers. Stabilization of Free 
Radicals on a Surface and Electron Paramagnetic Resonance in 
Quinhydrone Salts, V. B. Golubev, Yu. M. Boyarchuk and V. B. 
Evdokimov, pp. 331-332. 

Conclusion: Our work has shown that: 1) The surface stabiliza- 
tion of semiquinones takes place by a chemical reaction analogous 
to the homogeneous formation of quinhydrone salts. 2) Under 
normal conditions a quinhydrone salt is a diamagnetic semi- 
quinone dimer; the degree and heat of dissociation of the dimeric 
form have been determined for the free quinhydrone salt, its 
solution in alcohol, and the system quinhydrone-Ba(OH)--8H,0. 
3) The kinetics of the polymerization of semiquinone in alcoholic 
solution have been studied, and the order of the reaction and the 
energy of activation determined. 


Vol. 34, no. 4, April 1960. 


Structure of Crystalline Polymers, E. I. Kolobov, pp. 339-344. 

Abstract: Polymers are not completely crystalline, but contain 
amorphous regions in addition to the ordered regions. 

To clarify the nature of the distribution of amorphous and 
ordered regions in a crystalline polymer, a study has been made of 
the thermomechanical properties of such polymers in both 
isotropic and oriented states, and the thermal expansion of 
polymer films having different degrees of orientation has been 
determined. The results have been interpreted on the basis ef 
the x-ray data of Bunn and Alcock. 

Summary: On the basis of the present investigations it may be 
assumed that the crystalline phase is separated into individual 
localized microscopic regions at the temperatures of the second 
thermomechanical transition (numbered in sequence with rise in 
temperature). 


Effect of Sodium Molybdate on the Catalytic Decomposition of 
Hydrogen Peroxide in the Presence of Dichromate, E. Z. Nano- 
bashvili, pp. 351-354. 

Abstract: The present paper is devoted to an investigation of 
the combined catalytic effect of sodium molybdate and potassium 
dichromate on the decomposition of hydrogen peroxide in a 
homogeneous medium. 

Summary: 

1 Sodium molybdate causes practically no decomposition of 
hydrogen peroxide in the absence of dichromate, either in neutral 
or in acid medium. Increasing the molybdate concentration 
while keeping the dichromate content constant leads to an 
acceleration of the reaction in neutral medium, but to a retarda- 
tion in acid medium. 

2 Under the joint action of dichromate and molybdate, in- 
creasing the initial concentration of hydrogen peroxide accelerates 
its decomposition in neutral medium, but retards it in acid 
medium. 

3 The activity of dichromate and molybdate catalysts falls 
with time, and after reaching a limiting value remains practically 
unchanged. 

4 Addition of the MoO{~ ion to the Cr.0;- ion increases the 
activity of the latter by a factor of 1.54 in neutral medium, but 
diminishes it by a factor of 1.11 in acid medium (0.233 N H2SQO,) 


1353 


low 
ler 
ater 
> to 
orm 
nce 
the 
3 of 
ing 
rect 
for 
ear 
ary 
s of 
ob- 
em 
wer 
fies 
on- 
ha 
vill 
co- 
im, 
ent 
ec- 
ry, q 
‘ja- 
the : 
me 
of 
the 
cal 
the 
al 
al 
iry 
he 
he 
he 
m- 
to 
ad. 
1e. 
S. 
he 
ne 
fa 
of 
na 
of 
es. 
its 
UT 


and 1.22 in a more strongly acid medium (2.23 N H2SO,). 

5 The reaction represents a case of irreversible catalysis, as is 
demonstrated by the dependence of the rate of decomposition on 
the initial concentration of the hydrogen peroxide and by the 
change in the kinetic course of the reaction on repeating the re- 
action (second reaction). The activities of the mixed catalysts 
are clearly nonadditive. 


Thermodynamic Theory of Critical Phenomena in Three- 
Component Systems. II. Shape of the Constant Pressure 
Isotherms for the Coexistence of Two Phases in the Vicinity 
of the Critical Point, A. V. Storonkin and A. I. Rusanov, pp. 
354-356. 

Introduction: Analysis of the shape of coexistence curves in the 
vicinity of the critical point occupies an important place in the 
theory of critical phenomena in binary systems. On passing to 
ternary systems it is of the greatest interest to examine the curve 
for coexistence at constant temperature and pressure. 

Summary: 

1 It has been demonstrated that the curves for the coexistence 
of two phases in a ternary system at constant temperature and 
pressure are symmetrical in the vicinity of the critical point, and 
the arrangement of the tie-lines in this region has been discussed. 

2 It has been shown that in general the curvature of these 
curves and the changes in second differentials of the thermo- 
dynamic potential at the critical point of a ternary system have 
values which are finite and differ from zero. 


Heat of Formation and Physical Properties of Semiconducting 
Intermetallic Compounds, P. P. Otopkov and A. M. Evseev, pp. 
388-390. 

Introduction: A number of papers consider relationships be- 
tween crystal-chemical, physical and energetic properties of 
semiconducting compounds. Although several relationships 
have in fact been found, others would be most useful. Thus 
there is only a single rule for the energy gap AE: Namely that in a 
series of compounds, for example GaP-GaAs-GaSb or AlSb- 
GaSb-InSb, AZ) decreases with increasing atomic weight of the 
variable constituent. Ormont and Goryunova, have also noted 
certain correlations in AZ) values for the group of compounds 
with zine blende structure. In order to obtain a better relation- 
ship between AZ) and the properties of a solid we have made use of 
a hypothesis of Ioffe: That the structure of energy levels and the 
energy gap is determined by short-range order, and further that 
the properties of semiconductors depend on the nature of the 
forces acting between the constituent atoms. 

Summary: An empirical relationship was found between the 
energy gap and the energy of atomization for intermetallic semi- 
conducting compounds. It is of optimum type. [Translator’s 
note: This is a literal translation—the precise meaning is not 
clear. 


Similarity Method in Mass Transfer Studies, A. G. Usmanov and 
A. N. Berezhnoi, pp. 430-435. 

Abstract: Mass transfer by diffusion or thermal diffusion is 
often the principal mechanism of various processes, such as 
evaporation, and mixing and separation of gases. The kinetics 
of diffusion mass transfer govern partly or wholly the rate of 
chemical processes occurring in the intermediate or diffusion 
ranges. Study of thermal diffusion in gases provides information 
on intermolecular forces. Thermal diffusion is used in industry 
to separate gas mixtures. For these reasons diffusion and 
thermal diffusion have been the subjects of numerous experi- 
mental and theoretical investigations. These problems are con- 
sidered here from the viewpoint of similarity of molecular proc- 
esses. 

Summary: 

1 Generalized relationships for diffusion have been derived 
which permit the determination of coefficients of diffusion of 
various vapors in temperature and concentration ranges for which 
no experimental data are available. 

2 Generalized relationships for thermal diffusion of gases 
have been derived which may be used to calculate the separation 
of mixtures by thermal diffusion over a wide temperature and 
concentration range. 

3 Separation of binary mixtures by thermal diffusion is 
described by a linear relationship. 


Effect of Temperature on the Flashback Velocity of Laminar 
Flames, F. B. Moin and V. U. Shevchuk, pp. 441-443. 
Introduction: According to Lewis and Elbe, laminar-flame 
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flashback is described by the equation 


uL 


"Dz 


where g is the velocity gradient of the gas at the tube-wall, w the 
average velocity of the gas at flashback (flashback velocity), D 
the burner diameter, uz the normal velocity of propagation of the 
flame, and z the “penetration distance,”’ i.e., the distance from the 
wall at which the extinction of the flame becomes appreciable. 
Since uy depends on the temperature of the mixture and x may, in 
addition, depend on the wall temperature, these factors influence 
the flashback velocity of the flame. 

Few papers have been published on the influence of tempera- 
ture conditions on flashback. The turbulent flashback of 
hydrogen-oxygen flames has been studied, and the authors state 
that the temperature of the burner tip is an important factor. 

We have studied the influence of the temperature of the mix- 
ture and of the tube wall on the flashback velocity of laminar 
flames in a methane-oxygen mixture in 8-mm diameter quartz 
tubes. The apparatus consisted of a vertical burner placed in an 
electric heater, and measuring instruments. The tip of the 
burner was fitted with a special electric heater to maintain the 
burner wall at the required temperature, which was kept con- 
stant with a photoelectric regulator with an accuracy of +10 deg. 


OPTICS AND SPECTROSCOPY (Optika i Spek- 
troskopiia). Published by Optical Society of 
America, Washington, D. C. 


Vol. 8, no. 6, June 1960. 


Spectroscopic Investigation of the State of the Gas Behind a 
Shock Wave. III, F. S. Faizullov, N. N. Sobolev and E. M. 
Kudryavtsev, pp. 400-404. 

Abstract: The generalized spectrum line reversal method was 
used to measure the temperature of nitrogen and air behind « 
shock wave in the range of shock wave velocities from 1.9 to 4.3 
km per sec, and for gas pressures in the plug from 0.1 to 4 atm. 
It was established that there was satisfactory agreement between 
the experimental and the theoretical variation of the temperature 
of the plug with the shock wave velocity. The distribution of 
temperature along the plug was studied for various shock wave 
velocities. 


Theory of the Luminescence of Activated Crystal Phosphors, \. 
D. Potekhina, pp. 437-440. 

Abstract: The Williams method for calculating the spectra of 
erystal phosphors is analyzed quantum-mechanically, and it is 
shown that the “repulsion energy’’ determined by the Williams 
method for the excited state corresponds to the sum of the 5- 
energy and the exchange energy, and in contradistinction to the 
ground state, does not include the energy of the Coulomb overlap. 
Thus it is shown that the Williams method correctly accounts for 
the shift of the excited level of the impurity ion in the crystal field. 
The calculations carried out also prove that the 4d"-4d%5p transi- 
tion in the silver ion is responsible for the luminescence bands in 
the phosphor NaCl:Ag. 


Temperature Effect on the Photoluminescence of the Evaporated 
ZnS: Mn Phosphor, N. A. Vlasenko, pp. 445-448. 

Abstract: A study was made of the effect of temperature on the 
absorption and emission spectra of an evaporated ZnS:Mn phos- 
phor between 100 and 550 K. The temperature dependence of 
the relative quantum yield of luminescence for samples with 
different activator concentrations was also studied. On the 
hasis of these data, certain conclusions were drawn regarding the 
excitation mechanism of the luminescence center, the character of 
luminescent and the nonradiative transitions, and the nature of 
concentration extinction in the phosphor studied. 


Photo-Ionization Process in Silicon, K. I. Britsyn and V. %. 
Vavilov, pp. 452-455. 

Abstract: Data are presented on the spectral and temperature 
dependence of the quantum yield (generation of electron hole 
pairs) in silicon during absorption of photons with energies up to 
1.9 ev. A method was used which excluded the effects of capture 
centers. The temperature dependence of the quantum yield was 
observed in the region of collision ionization effected by carriers 
liberated in the absorption of photons. As the temperature in 
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creased, the boundary of this region shifted toward lower photon 
energies, partly because of a narrowing of the energy gap in silicon. 


Calculation of Intensities in Infrared Spectra of Molecules, B. N. 
Snegirev and M. A. Kovner, pp. 462-463. 
Vol. 9, no. 3, Sept. 1960. 


Concentration of Atoms in the Arc Plasma as Determined by the 
Pulse Absorption Method, L.. D. Kondrasheva and I. V. Pod- 
moshenskii, pp. 149-151. 


Fig. 1. Schematic diagram of the equipment. 


Abstract: This paper describes a pulse method which gives the 
instantaneous absorption spectrum of a plasma at temperatures 
up to 6000 K. It was found that in the Fe-Cu direct current are, 
the volume occupied by atoms excited to about 3 ev comprises 
only about 30% of the volume occupied by unexcited atoms. The 
instantaneous concentrations of Fe, Cr and Mn were determined 


Radiation Transfer in a Homogeneous Medium for the Non- 
conservative Case, V. I. Barkov, pp. 195-200. 

Abstract: A solution is given for the integro-differential equa- 
tion of radiation transfer. The case of a point source located 
within a homogeneous medium filling all space is examined. A 
general expression for the average intensity and flux of radiation 
is obtained in the nth approximation. 

Introduction: The method of solving the integro-differential 
equation of transfer used in this paper is based on the replacement 
of integrals by finite sums. The particular merit of this method, 
which is essentially equivalent to the spherical harmonie method, 
is its rapid convergence in approximations of low order. 


Principles of the Theory of the Synthesis of Coated Absorbing 
Layers, P. G. Kard, pp. 200-204. 

Abstract: Optical layers are investigated which consist of a 
number of absorbing layers separated by nonabsorbing inter- 
layers. Necessary and sufficient conditions are found for the 
composition of the system in order that the transmitted light be 
maximized, The number, location, thickness and optical param- 
eters of the absorbing layers are assumed to be prescribed; 
whereas the conditions for the coating are related to the non- 
absorbing interlayers. It is shown that the effect of the coating 
depends only on the properties of the absorbing layers, and a 
simple relation is found for this dependence. 


Optical Systems With Enhanced Resolving Power, A. I. Karts- 
shev, pp. 204-206. 

Abstract: Two new methods of optical image formation are 
described which under some conditions enable resolving powers to 
be attained that are higher than those obtained by the usual 
methods of observation of optical images. 


The Franck-Condon Factors for the Schumann-Runge Band 
System in the Oxygen Molecule, I. T. Yakubov, pp. 212-213. 
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-1 o-2 o-4 o-5 «-6 
Comparison of Franck-Condon factors cal- 
culated by different methods. v’-value: 1— 
3:v’ = 3; 4—6:v’=0. Method of calculation: 
1, 4—ref. 4; 2, 5—ref. 3; 3, 6—this paper. 
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PHYSICS OF METALS AND METALLOGRAPHY 
(Fizika Metallov i Metallovedenia). Published by 
Pergamon Institute, New York. 


Vol. 8, no. 5, 1959. 

Effect of Fluctuations in the Parameters of Correlation Upon 
Scattering of X-Rays and Thermal Neutrons by Solid Solutions, 
M. A. Krivoglaz, pp. 8-26. 

Abstract: Scattering of x-rays and thermal neutrons by sub- 
stitutional solid solutions associated with a variety of atomic 
scattering factors and geometric distortions are discussed. 
Heterogeneity associated with fluctuating changes not only in the 
composition but also in the parameters of correlation is con- 
sidered. It is shown that accounting for fluctuations in the 
parameters of correlation order may lead to substantial changes 
in the shape of isodiffusion curves. Ideal solutions and solutions 
near the point of the phase transition of the second order are dis- 
cussed in detail. Formulas are given for the quantities which 
determine the weakening due to static distortions in intensity of 
normal reflections. 

Conclusions: Simultaneous calculations of scattering from 
fluctuations in the composition and in the parameters of correla- 
tion (p.c.) lead, firstly, to a certain change of the coefficient correl- 
ating the lattice distortions with the fluctuations in the composi- 
tion, and, secondly, to the occurrence of an additional scattering. 
Isodiffusive curves of scattering from fluctuations in the p.c. are 
not of lemniscate but close to oval shape. The effects associated 
with fluctuations in the p.c. play a relatively important part when 
the lattice spacing greatly varies with the p.c. and deviations from 
Vegard’s principle are large. In ideal solutions the scattering 
power of fluctuations in the p.c. is proportional to c? (1—c)? and 
becomes insignificant for small concentrations. In non-ideal solu- 
tions the part played by these fluctuations may strengthen. Near 
the phase transition point, in particular, they may lead to the oc- 
currence of strong diffusive scattering in the vicinity of structural 
reflections (in addition to the usual scattering in the vicinity of 
superlattice reflections) the intensity of which tends to infinity 
with 

Adiabatic Galvanomagnetic and Thermomagnetic Phenomena in 
Semiconductors. II. Mixed Conductivity, V. Ye. Vzdornov and 
I. M. Tsidil’kovskii, pp. 31-37. 

Abstract: A theoretical examination is made of adiabatic 
galvanomagnetic and thermomagnetic phenomena in semicon- 
ductors with mixed conductivity, based on the hypothesis of the 
isotropy of effective masses and the relaxation times of con- 
ductors. 

Conclusions: 

1 The nature of the dependence of all transverse galvano- 
magnetic and thermomagnetic effects on the magnetic field for 
semiconductors with unequal concentrations of electrons and holes 
is the same as in the case of impurity conductivity: in the region 
of weak fields these effects are proportional to H, while in the 
region of strong fields they are inversely proportional to H (in the 
case of the Hall effect this is true not of the field Fy, but of the 
current jy). 

2 Where there are equal concentrations of conductors the 
transverse effects may be divided into two groups: Hall and 
Righi-Ledue effects, for which the nature of the magnetic field 
dependence is the same as in the case of unequal concentrations, 
and Ettingshausen and Nernst-Ettingshausen effects, which are 
directly proportional to H whatever the magnetic field. 

3. The nature of the dependence of all linear effects on the 
magnetic field is the same as in the case of impurity conductivity. 
Characteristic properties of the Maggi-Righi-Leduc effect are to 
be observed in the presence of strong fields and equal concentra- 
tions of electrons and holes. While the dependence of electronic 
thermal conductivity on the magnetic fields, obtained when Ey = 
0, is identical for different correlations of concentrations, when 
Jy = Oand 


Ni = N- = Oand [(u + H)/c]? > 1 


the electronic thermal conductivity tends not toward zero, as in 
the case of N; # N_, but toward some constant value. When the 
values of € are sufficiently high (as for example in Ge and Si) the 
symbol for the Maggi-Ledue effect may even vary; i.e., there 
must be an increase in thermal conductivity when the transition 
from impurity to natural conductivity occurs. 

4 Attention must be paid to the essential difference between 
adiabatic effects and corresponding isothermal effects, based on 
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their relation to the width of the forbidden energy band. Thus, 
the Hall field and resistance variation in a magnetic field will 
contain terms proportional to ¢?, while in isothermal conditions 
they are not dependent on «. The Nernst-Ettingshausen trans- 
verse effect in weak magnetic fields and in strong fields when N+ + 
N—- contains terms ~ e*, but in isothermal conditions it contains 
terms ~ e. Electronic thermal conductivity in both weak and 
strong magnetic fields when N+ contains components ~ 
while in isothermal conditions Aen ~e?. 

Thus, in equations for adiabatic effects ¢ is at a power twice as 
high as in equations for corresponding isothermal effects. 


Dependence of the Field of a Defect of the Crack Type on the 
Thickness of the Covering Metal Layer, V. V. Vlasov and R. Ye. 
Yershov, pp. 47-51. 


Determination of Surface Tension on the Boundaries of Phases, 
I. V. Salli, pp. 76-78. 


High Temperature Study of Alloys the Phases of Which Have 
Identical Lattices With Nearly the Same Parameters, I. I. 
Titarenko and B. M. Rovinskii, pp. 84-87. 

Abstract: The results of a study of refractory nickel alloys at 
temperatures of up to 1250 deg are given. 

Conclusions: 

1 The maximum difference in the parameters of the lattices of 
the 7 and 7’ phases in the alloy ZhS6 are about 0.002 kX both at 
high temperatures (1125 deg and after 1 hr holding at 1200 deg) 
and in their electrolytic separation. 

2 Diffusion processes which occur in the alloy ZhS6 at a tem- 
perature of 1125 deg lead to the segregation of the lattices of the 
and phases. 

3 In the process of holding the alloy at a temperature of 1200 
deg diffusion displacement of the components from the y’ phase 
into the solid solution occurs. 

4 On the cooling the solid solution from a temperature of 1250 
deg the segregated y’ phase is first precipitated, but at a lower 
temperature the “‘interlinking”’ of the lattices of the y and y’ 
phases occurs. 


RADIO ENGINEERING AND ELECTRONICS 
(Radiotekhnika i Elektronika). Published by Per- 
gamon Institute. New York, in conjunction with 
Massachusetts Institute of Technology. 


Vol. 4, no. 11, 1959, 


Structure of the Matrix and Properties of Reactive Four-Terminal 
Networks, L. R. Iavich, pp. 57-64. 

Abstract: The general case of irreversible reactive four- 
terminal networks is considered. It is proved that if employing 
an ideal power transformer it is possible to extract from it a 
reversible four-terminal network having known properties. 


Theory of Double-Curvature Reflectors, L. B. Tartakovskii, pp. 
95-110. 

Abstract: The properties of the envelope of a single-parametric 
family of confocal paraboloids are studied and simple formulas 
are proposed for the design of reflector antennas with fan-shaped 
beams. The differences between the vertical patterns of a double- 
curvature reflector and a cylindrical reflector with the same curve 
in the central section are analyzed. The position of the equiv- 
alent linear antenna is found for which it is possible to apply the 
aperture method to correct the curve of the central section. A 
specific method for this correction based on the application of the 
phase diagrams of linear antenna theory is illustrated. 


Plasma Electron Resonance, M. D. Gabovich and L. L. Pasech- 
nik, pp. 142-148. 

Abstract: Plasma electron resonance observed in interaction of 
a modulated electron beam with an independently established 
plasma is described. 


~oe i With an Electronic Commutator, Iu. E. Aptek, pp. 
184-193. 

Abstract: An RC commutator (synchronous filter) with a 
commutator tube is described. The precision of detecting 
periodic signals is analyzed. The requirements for the commu- 
tator tube are derived. 

Conclusions: The analysis and the experimental verification 
have shown that an RC commutator with a commutator tube is, 
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in principle, similar to the mechanical RC filters and to filters with 
individual switching tubes. The circuit described was used for 
studying the quality of modulation in a radiometer. In spite of 
the small modulation index (the intensity of the hot standard was 
only a small percentage of the noise intensity of the receiver) a 
good reproduction of the signal form was obtained. However, 
the technical characteristics of the commutator tube are such that 
it is difficult to obtain 7; min << 2kQ, ri max > 500 MQ; while in the 
case of mechanical commutators 7; min < 2Q, ri max> 105 MQ. In 
an electronic commutator it is also impossible to prevent com- 
pletely an overlap between the chambers. Therefore the preci- 
sion of reproducing a signal in an RC commutator with commu- 
tator tube is considerably lower than with a mechanical commu- 
tator. The advantages of the system are: Simplicity, facilities 
for synchronization, a wide range of frequencies fo. It can be 
assumed that with industrial production of commutator tubes 
(necessary also for other purposes) the RC filter circuit with such a 
tube can find applications. 


Some Problems of Producing Stable Secondary Electron Multi- 
pliers Resistant to Exposure to Gaseous Media, A. M. Tiutikov, 
pp. 195-203. 

Abstract: The results are given of an investigation into the 
stability of the amplification factor of secondary electron multi- 
pliers based on copper-beryllium and copper-magnesium all sys as 
a function of long-term exposure to air or in some other gaseous 
medium and as a function of the activation conditions; rec- 
ommendations are given for the selection of the alloy and for the 
technology of manufacture insuring the highest stability. 


Initial Energy Distribution of Photo-Electrons Emitted From 
Semitransparent Silver-Cesium-Oxide Photocathodes, N. A. 
Soboleva, pp. 223-236. 

Abstract; As a consequence of distortions of the form of the 
retarding-potential curves, caused by the presence on the surface 
of spots with unequal work function, the usual method of the 
spherical collector is unsuitable for the analysis of the energy dis- 
tribution of photo-electrons emitted by cesium-oxide photo- 
cathodes. 

To obtain the initial energy distribution of photo-electrons from 
these cathodes, a method was employed using a three-electrode 
spherical system with an additional accelerating grid in the 
vicinity of the cathode. 

Distribution curves of photo-electrons within a wide spectral 
range were obtained by illuminating semitransparent cathodes 
both from the front and the rear. 


Contribution to the Theory of the Barkhausen Effect in a Periodi- 
cally Varying Field, F. V. Bunkin, pp. 237-247. 

Abstract: A necessary criterion is formulated for the physical 
realizability of a model, describing the statistical nature of a 
cylical magnetic polarity reversal of ferromagnets. Two possible 
models for a single-domain specimen are discussed. Expressions 
are obtained for the noise spectrum of the induction emf. 


Oscillography of Relativistic Electron Beams, L. E. Tsopp, pp. 
281-284. 


Stabilization of the Frequency of a Quartz Oscillator by a Maser, 
I. D. Murin, pp. 293-297. 
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FIG. 1. Block diagram of the installation for the automatic 
control of the frequency of a quartz oscillator by a 
molecular oscillator 


ARS JourNaL SUPPLEMENT 


: els 4U=-xaf 
= le Rk, 
=f 
nis bh 
le 
! 1 M 


with 
d for 
ite of 
1 was 
er) a 
ever, 
that 
1 the 

In 
com- 
reci- 
amu- 
amu- 
lities 
n be 
ubes 
ich a 


[ulti- 
ikov, 


» the 
\ulti- 
as 
cous 

rec- 
r the 


‘rom 
A. 


the 
rface 
the 
dis- 
10to- 


from 
rode 

the 
etral 
odes 
iodi- 
sical 
of a 


sible 
ions 


pp. 


ENT 


Introduction: In the present article the description of a system 
for stabilizing the frequency of a quartz oscillator by a maser is 
given; a system developed for the molecular time (frequency) 
standard established at the Lebediev Institute of Physics of the 
Academy of Science of the USSR in 1958. The basis for this 
standard is the maser, which has an absolute frequency stability 
of the order of 10-*. Since the power of the maser is small (10~%- 
10-“w) and because it is difficult to realize the amplification and 
direct scaling of its frequency (23, 870, 129, 420 eps), a quartz 
oscillator with its frequency regulated by the frequency of the 
maser is used as the frequency standard. 


SOVIET PHYSICS—ACOUSTICS (Akusticheskii 
Zhurnal). Published by American Institute of 
Physics, New York. 


Vol. 6, no. 3, Jan.-March 1961. 


Flow Regulation of Neural Impulses in the Auditory System, G. V. 
Gershuni, pp. 298-305. 

Abstract: On the basis of an electrophysiological investigation of 
animals to study the responses of various branches of the auditory 
system (cochlea, inner geniculate body, auditory cortex), the 
author treats the problem of an organism’s paths of regulation of 
the flow of neural impulses carrying information on an acoustic 
signal. Data are presented showing the significance of the various 
paths of this regulation namely: 1) Limiting of the flow of im- 
pulses as the result of the very process of traversal by the neural 
impulse; 2) limiting of the flow of impulses as the result of action 
of a special system of (reverse) centrifugal connections; 3) varia- 
tion of the flow of impulses in the higher branches of the auditory 
system as the result of influences coming from nonauditory 
branches of the central nerve system. Two mechanisms are 
considered for the regulation of information transmitted by the 
flow of nerve impulses arising in the acoustically stimulated 
auditory system: 1) Variation in the number of elements partici- 
pating; 2) variation in the intrinsic noise levels of the system. 


Optimum Operating Conditions for a High Power Concentrator, 
K. A. Naugol’nykh and L. D. Rozenberg, pp. 352-355. 

Abstract: The authors discuss the operation of a spherical 
focusing system when the intensity of the focused sound is high, in 
which case nonlinear effects distorting the wave form become im- 
portant. They determine the working conditions of a concentra- 
tor that correspond to maximum velocity amplitude of the wave 
at the focus of the system. 


Propagation of Finite-Amplitude Waves in Relaxing Media, A. L. 
Polyakova, pp. 356-359. 

Abstract: The propagation of sound having a finite amplitude in 
a relaxing medium is considered in the second approximation, and 
it is shown that if p:/po<2Ac/c, i.e., if the sound amplitude is less 
than a value which is characteristic of dispersion and relaxation 
absorption, then the behavior of the second harmonic will be 
determined primarily by absorption. In this case those effects 
due to phase shift as the result of dispersion appear only at fre- 
quencies near wr = 1/ »/2, when the influence of dispersion comes 
to a maximum, the influence of absorption to a minimum. 


Ultrasonic Waves Acting on Biomacromolecules, [. FE. El’piner. 
pp. 399-406. 

Abstract: Under the action of ultrasonic waves, biomacro- 
molecules are not only subjected to depolymerization processes, 
but under certain definite conditions an increase in molecular 
weight of these polymers takes place. The nature of these 
changes depends on the nature of the gas present in the sound- 
irradiated solution of the substances investigated. When oxygen 
is present, the biomacromolecules are predominantly observed to 
suffer processes of degradation, which in a number of cases cause 
oppression of the functional characteristics (biocatalytic activity ) 
of the macromolecules. When hydrogen is present, on the other 
hand, an increase in molecular weight of the irradiated biopolymers 
is observed. In this case the intrinsic (natural) properties of the 
molecules are preserved. The investigations are of interest from 
the point of view of studying the essential connections between 
structure and form of biomacromolecules and their biological 
function, which is one of the foremost problems of modern bio- 
physics. 
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Vol. 6, no. 4, April-June 1961. 


Propagation of Tone Signals in a Plane Layer of Water With 
Waves on the Surface, L. N. Zakharov, pp. 455-461. 

Abstract: The effect of a wavy surface on the amplitude and 
phase difference of acoustic tone signals propagating in a plane 
layer of water is considered. On the basis of a statistical analysis 
of the experimental data obtained, quantitative relations are 
established between fluctuations in the acoustic signal and ampli- 
tude of the surface wave. 

Introduction: The principal aim of the present experimental 
work was to establish a connection between amplitude and phase 
fluctuations in the acoustic signal and state of the free surface 
The investigations were made in a plane (in the horizontal), 
almost unbounded layer of water under so-called shallow water 
conditions, when the depth of the reservoir is comparable with 
the acoustic wave length and is much less than its own length. 

Summary: 

With the stipulation of a plane shallow water layer (fresh water 
reservoir), unbounded in the horizontal, a fixed receiver and 
radiator, and total isothermicity of the layer, the main cause of 
fluctuation in the stationary acoustic signal is roughness of the 
water surface. The following results were established in the 
work: 

1 The mean periods of the amplitude and phase difference 
fluctuations agree with the mean periods of the surface waviness. 
At low frequencies (f = 0.5-2 ke) the curves for surface roughness 
and amplitude and phase fluctuations are in complete agreement 
with one another. 

2 When the roughness is slight the amplitude and phase 
fluctuations in the acoustic signal increase in proportion to the 
height of the wave. 

3 With an increase in separation (B) between pickups (the 
experiment was carried out at frequencies of f = 0.9-1.6 ke with a 


roughness AH? = 3 cm), an increase in the phase difference 
fluctuations is observed at first, then, beginning with B = 80-100 
em, the rise in fluctuations eases off and saturation sets in. 

4 The largest fluctuations in the acoustic signal were observed 


for the smallest separation between source and pickups (L. = 15 
m). With an increase in distance the fluctuations fall off, and at 
L = 150-200 m their difference in values becomes comparable 


with the spread in experimental points. 

5 The normalized autocorrelation functions for the amplitude, 
phase difference and surface roughness have an_ oscillating 
character. 

6 With increasing frequency and constant roughness or with 
increasing wave height and constant signal frequency (the dis- 
tance between bases in these cases remaining fixed), the phase 
difference and amplitude autocorrelation time intervals for the 
acoustic signal decrease. 


Analysis of Acoustic Radiation From a Turbulent Aerodynamic 
Flow, L. M. Lyamshev, pp. 472-476. 

Abstract: Starting from Lighthill’s equation, a method is 
presented for calculating the acoustic radiation of a turbulent 
aerodynamic flow in free space and in the presence of elastic sur- 
faces outside the flow. It is shown that if one knows the solution 
to the nonstatistical problem of diffraction of a spherical (or 
plane) sound wave by the said elastic surfaces, and the correlation 
function for pressure fluctuations in the flow, calculation of the 
spectral intensity of the radiation from the turbulent flow reduces 
to quadraties. 


Radiation of Sound in a Liquid With Cavitation Present, |.. ID). 
Rozenberg and M., G. Sirotyuk, pp. 477-479. 

Abstract: It is experimentally demonstrated at a frequency of 
21 ke that the radiation resistance of a sinusoidally oscillating 
adiator, as measured on the electrical side, begins to fall off 
sharply at the instant that cavitation sets in. With a further 
increase in particle velocity and corresponding growth in cavita- 
tion, the specific radiation resistance attains a value of 0.3 poc and 
then becomes stabilized. The radiated acoustic intensity at first 
increases in proportion to the square of the radiator particle 
velocity, then its growth slackens and ceases altogether (in our 
experiments) at a level of 1.5 w per cm?. With an increase in 
particle velocity above 25 em per sec radiation intensity once 
again begins to increase in proportion to the square of the particle 
velocity, but with one third the slope as before. 


Dependence of Cavitation Erosion on the Solubility of a Gas 
Under a Liquid, A. S. Bebchuk and L. D. Rozenberg, pp. 496-498. 
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SOVIET PHYSICS—DOKLADY (Doklady Aka- 
demii Nauk SSSR, Otdel. Fiziki). Published by 
American Institute of Physics, New York. 


Vol. 5, no. 2, Sept.-Oct. 1960. 


Periodic Modes in Systems With Nonlinear Pulse Elements, M. 
M. Simkin, pp. 227-230. 


Calculation of Oblique Shock Waves in Magnetogasdynamics, 
M. I. Kiselev and N. I. Kolosnitsyn, pp. 246-249. 


Isothermal Scattering of a Gas Cloud, V. S. Imshennik pp. 253- 
256. 

Introduction: It is of interest in many astrophysical problems 
to study the diffusion of a cloud of gas, and where this diffusion 
occurs under the influence of internal irradiation. The law of 
conservation of energy does not hold in the cloud under these 
circumstances, and the entropy of the gas particles increases with 
increasing diffusion of the gas. Diffusion of this type has, in 
general, a complex character, but there is one simple limiting case 
of nonisentropic diffusion—that of isothermal diffusion. In 
investigating this, it will be assumed that the temperature of the 
gas does not depend on the space coordinates, but depends in an 
arbitrary known way on the time, i.e., 9 = 0(t). The dependence 
of the gas temperature on the time is determined by the character 
of the external radiation. An isothermal approximation has been 
used before to obtain certain accurate solutions of the one- 
dimensional gasdynamic equations. In the present paper, we will 
show that, in the isothermal case, there is an asymptotic solution 
of the gas cloud diffusion problem. Self-gravitational forces and 
magnetic fields will be assumed to be absent. 


Numerical Solution and Numerical Analysis of the Solutions to 
Homogeneous Systems of Linear Algebraic Equations of General 
Type, L. S. Mayants, pp. 257-260. 

Introduction: The numerical solution of homogeneous systems 
of linear algebraic equations having the form (in matrix notation) 


(A — =0 


(where A and B are nth order matrixes, X is a column vector with 
n components, is a number) is especially difficult when both 
matrixes A and B are singular. The method presented is useful 
not only in the case of nonsingular matrixes A and B, but also in 
cases when one or both simultaneously are singular. 


New Kinetic Equations in the Theory of Monatomic Gases, 8. V. 
Vallander, pp. 269-271. 

Introduction: In the present paper, we state some new kinetic 
equations in the theory of monatomic gases. Unlike the Boltz- 
mann kinetic equation, our equations are pure integral equations 
and contain boundary conditions at the surfaces of the bodies 
around which the gas flows. Owing to the inclusion of the 
boundary conditions in the kinetic equations, these equations 
can be regarded as the equations of the kinetics (aerodynamics) of 
rarefied gases. The Boltzmann kinetic equation is a decidedly 
impoverished consequence of our equations. To avoid giving 
attention to secondary details, we shall assume that there is no 
external field of mass forces, and that all the atoms of the gas are 
identical, and shall take into account only collisions of pairs of 
atoms. 


Laws of Thermal Emission of a Flame, B. I. Plyukhin, pp. 279- 
282. 

Introduction: According to experimental evidence, the laws 
often used for temperature measurements—Wien’s law, the 
Stefan-Boltzmann law, etc.—are far from being always valid in 
the conditions of a flame. The theoretical examinations of the 
problem are essentially limited to general assumptions regarding 
the blackness of a flame. 


Acceleration of the Plasma of a Discharge and Production of 
Strong Shock Waves in a Camera With Coaxial Electrodes, S. R. 
Kholev and D. 8. Poltavchenko, pp. 356-360. 

Introduction: In a series of experimental investigations 
different methods for accelerating a plasma interacting with an 
electromagnetic field were studied. The method described here 
for accelerating a plasma consisted basically in a comparatively 
long interaction (lasting, e.g., one or a few periods of the current) 
of a radial discharge current with a concentric magnetic field. 
The purpose of the investigation was not only the study of a 
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Fig. 1. E) Accelerating elec- 
trodes; P) and P,) dischargers; 
initiating capacitor; T) 
pulse transformer; 4-4 
connections of the transformer 
with the discharger P. 


method of plasma acceleration, but the production of strong 
shock waves. In this case, a particular energy redistribution can 
take place: The gas, condensed by the plasma, can itself become a 
plasma with a still higher temperature. 

The accelerating device was built in the form of two coaxial 
electrodes having a cylindrical shape (Fig. 1A): a central rod abc 
and a tube def. The discharge of the condenser battery C caused 
a current, whose direction is indicated in Fig. 1A by arrows, to 
flow in the circuit. The radial discharge current J interacted 
with the concentric magnetic field H of the central electrode. 
The ponderomotive force F = [IH] thus arising was always 
directed outward from the accelerating electrodes, independent of 
the direction of the current. Therefore, the magnetic field must 
act on the discharge plasma as a piston would 


Study of the Interplanetary Ionized Gas, High Energy Electrons, 
and Corpuscular Radiation From the Sun by Means of the Three- 
Electrode Trap for Charged Particles on the Second Soviet 
Cosmic Rocket, K.I. Gringauz, V. V. Bezrukikh, V.D. Ozerov and 
R. E. Rybchinskii, pp. 361-364. 


Diffraction Modulator of Light With Several Ultrasonic Counter 
radiators, Yu. V. Popov and I. I. Adrianova, pp. 390-392. 


a b c 
Absorber 


ROOY 


Reflector Radiator Radiator 


| 


Radiator 


Absorber 


Fig. 1. Varieties of high-frequency diffrac- 
tion light modulators: a) modulator with re- 
flector; b) modulator with traveling waves; 
c) modulator with ultrasonic counterradi{ators. 
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Introduction: There are two well-known types of high frequency 
diffraction modulators of light: A modulator with a reflector in 
which standing waves are formed as a result of superposition of 
waves incident on, and reflected from, the reflector (Fig. la); anda 
modulator in which standing waves are simulated by two traveling 
waves propagated in opposite directions (Fig. 1b). In the re- 
Hector-type modulator, there is a large phase inhomogeneity in 
the modulated light beam, and the depth of modulation depends 
strongly on the modulation frequency and on external conditions; 
in the traveling wave modulator, the inhomogeneity is less, but 
the modulation depth, although not depending on frequency and 
external conditions, amounts to about 30% of the maximum 
depth for the reflector-ty pe modulator. 

The diffraction modulator with several ultrasonic counter- 
radiators has, in comparison with the reflector-type modulator, a 
higher efficiency and improved phase characteristics and, in 
‘omparison with the traveling wave modulator, a greater modu- 
lation depth, the magnitude of which, even for work at fixed 
frequencies, can be held constant. 


Stability of a Cantilever Cylindrical Shell With Reinforced Edge 
Subjected to External Pressure, V. M. Darevskii and R. I. 
{X\shnyakin, pp. 423-426. 

Introduction: The question of the stability of a cantilever 
circular cylindrical shell whose movable edge is reinforeed by an 
elastic ring, subjected to an external pressure g > 0 is investi- 
gated. For simplicity, we consider the cross section of the ring 
to be a rectangle with base a and height H. The interaction of 
the ring and the shell is characterized by forces 7), S, N, and 
moments M,, A, acting in the cross section separating the shell 
irom the ring. In the subcritical state, of these forces only N, and 
\f, are different from zero. Upon loss of stability, all of them, 
in general, do not vanish, and may be represented in the form of 
sums of fundamental and supplementary values. Separating the 
mentioned forces and moments into directions corresponding to 
the shell and the ring, one may consider then isolated from each 
other. 


TELECOMMUNICATIONS (Elektrosviaz’). Pub- 
lished by Pergamon Institute, New York, in 
conjunction with Massachusetts Institute of 
Technology. 


No. 12, 1959. 


Bilateral Amplitude Modulation and Its Application to the 
Simultaneous Transmission of Two Signals, B. M. Pevzner, pp. 
1329-1341. 

Abstract: Fundamentals are given of bilateral amplitude modu- 
lation, as a method of transmitting two signals with a single 
carrier and with the level selection at the receiving end. The 
following in particular are analyzed: cross-talk distortion be- 
tween the transmitted signal in the quasi-stationary state and its 
dependence on the width of the frequency spectrum, on the 
depth of modulation and on the type of the sub-carrier of the 
system used. Possibilities for compensating distortion are 
examined and methods of obtaining bilateral modulation discussed. 

Conclusions: The advantages of bilateral modulation are: 
Simple signal separation at the receiver; noncritical phase 
characteristics of the channel. The disadvantages of bilateral 
modulation are: Considerable cross-modulation distortion; 
ciitical frequency characteristics of the channel (wrong charac- 
teristics will lead to changes in the magnitude of a); impossi- 
bility of mixing-in of low frequency (third) signals. It is rec- 
ommended to use bilateral modulation in cases when the two 
transmitted signals are being obtained from the same source and 
are to be used simultaneously (color television, stereophonic 
sound, ete.) and consequently when relatively large cross-modula- 
tion distortions are permissible and the simplification of the 
receiver is required. The efficiency and the interference-killing 
properties of bilateral modulation are intended to be discussed, 
together with other methods of multiplex transmission, later on. 


Determining the Stability of Transistor Triode and Valve Circuits 
from the Real Component of their Terminal Imittance, A. A. 
Kulikovskii, pp. 1368-1380. 

Abstract: It is shown that a consideration of the real com- 
ponent of the terminal immittance of an active linear circuit enables 
us to determine the stabilitity of the latter, when the terminals 
across which the immittance is measured are correctly chosen. 
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The conditions for which potential instability is converted into 
actual instability, analogous to the conditions for amplitude 
und phase balance, are obtained; a generalized stability margin 
coeflicient is introduced; it is shown that the choice of the 
parameter matrix and equivalent circuit of an actual four- 
terminal network depends on the type of instability inherent in 
the network. 


Input and Output Reactive Circuits of Amplifiers, Taking Into 
Consideration the Finite Value of the Transformer Inductance, 
A. E. Znamenskii, pp. 1381-1387. 

Abstract: Input and output reactive amplifier circuits con- 
taining a finite number of circuit elements are considered from the 
point of view of obtaining the maximum current and voltage 
amplification. The inductance of the transformer windings is 
assumed to be finite. 


No. 4, 1960. 


Certain Questions of the Theory of Oral Information, Yu. F. 
Pelegov, pp. 366-376. 

Abstract: On the basis of the proposed equivalent circuit for 
speech communication, individual elements of the theory of 
speech have been developed. The question of redundencey and its 
filtration has been rendered more exact. It has been shown that 
the problem of the objective recognition of a phoneme can be 
analyzed by the methods of the theory of the “‘ideal receiver.” 

Discussion of Results: 

The elements of the theory of speech communication set out 
here allow us to explain a number of experimental facts which, in 
their turn, confirm the proof of the fundamental assumptions. 

1 The clipping of speech. It is well known that in clipping 
there is an inconsiderable distortion of the correlation function 
determining, as we have shown, the orientation of the vector of 
the controlling signals. Consequently the amount of meaningful 
information changes very little, and in the final analvsis this leads 
to a high degree of distinctness. Note that a possible objective 
method of estimating the distortion consists of the measurement 
of the rms errors of the instantaneous correlation function. 

2 Preservation of distinctness on Gabor transformation. The 
Gabor transformation of speech provides for the preservation of 
portions of a duration in one or several periods of the fundamental 
note. On the other hand the duration of the short-time correla- 
tion function is either considerably less (unvoiced consonants, 
whispered speech), or this function is periodic with a period equal 
to that of the frequency of the fundamental tone. In both cases 
the Gabor transformation does not distort the correlation func- 
tion which ensures the conservation of distinctness. 

3 Compression of the spectrum by the “forming vocoder.” 
We have seen that although the complete spectrum of the con- 
trolling signals can be fairly wide, information as to the ‘‘form”’ 
of the phonemes is transmitted by only part of the neurons. 
Their number includes a group which characterizes the funda- 
mental resonance frequency and the goodness of the resonator of 
the mouth cavity, and also the degree of tension of the vocal 
chords. The construction of the ‘forming vocoder’’ is based on 
the transmission of exactly this group of signals, which also causes 
the compression of the spectrum. 

4 High distinguishability of the oscillograms of ‘‘visible’’ 
speech. Visible speech is the description of the short-term energy 
spectrum, but the latter is analogous to the transmission function 
of the linear circuit for unvoiced consonants, or contains supple- 
mentary information as to the frequency of the fundamental tone. 
In both cases we have a complete characteristic of the vector of 
the controlling signals which carries all the information. 


Application of Nonminimal Phase Circuits in Certain Systems 
With Phase Transformations, A. D. Artyn, pp. 377-388. 

Abstract: Practically useful nonminimal phase circuits have 
been described, convenient, in particular, for wide-band phase 
transformation. The possibility has been shown of using such 
circuits as effective phase modulators. 


Calculation of a Phase Inversion Cascade With Semi-Conducting 
Triode, IX. V. Gribov and Yu. G. Kryukov, pp. 388-401. 

Abstract: On the basis of the analysis of the equivalent circuit 
of a plane semiconducting triode, the fundamental relations are 
deduced for the calculation of a phase inversion cascade. Ques- 
tions of matching with the output twofold cascade of class B using 
nonidentical semi-conducting triodes are considered. It is shown 
that, in this case, it is possible to work with small nonlinear distor- 
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tion. The succiessve steps in engineering calculation are given, 
and the results of its experimental verification over a wide range 
of temperature. 

Conclusions: Investigation and experimental verification have 
shown that the use of a phase inverter with semi-conducting 
triodes makes it possible to use, in the ouput, double-contact 
cascade semi-conducting triodes without special selection, giving 
small nonlinear distortion and considerable improvement in the 
frequency characteristic of the amplifier of low frequencies. The 
circuit described can be used within a wide temperature range 
from — 50 to 70 C and could be applied both to radio broadcasting 
and also to special apparatus. 


No. 5, 1960. 


Transient Processes in a Trigger Circuit With Semi-Conducting 
Triodes, P. N. Faizulaev, pp. 402-415. 

Abstract: Transient processes in a symmetrical trigger cir- 
cuit with plane semiconducting triodes on a “‘strong’’ starting sig- 
nal are considered. The optimal value of the accelerating capacity 
from the point of view of the greatest speed of action has been de- 
termined. An estimate of the limiting speed of action of the trig- 
ger circuit is given. 


Increasing the Effectiveness of Multichannel Systems With Time 
Separation of Channels, G. A. Levin, B. R. Levin, K. I. Aizen- 
berg and V. 8. Rozanov, pp. 487-496. 

Abstract: Possibilities are indicated of increasing the carrying 
power of impulse systems of communication by using the statis- 
tical instantaneous values of the signal and the activity of the 
telephonic channels. A theoretical investigation of interval time 
impulse modulation is carried out making it possible to realize 
this increase. 


Low Frequency Spectrum of Noise at the Output of a Clipper, 
B. N. Mityashev, pp. 497-504. 

Abstract: The width of the low frequency portion of the spec- 
trum of the noise passing through a high frequency amplifier and a 
clipper has been calculated. The calculations have been carried 
out for one-sided and two-sided clipping. 


Loss of Amplification of High Directional Antennas in the Path of 
Tropospheric Scattering, M. F. Bakhareva and B. E. Kinber, pp. 
573-577. 

Introduction: When high directional antennas are used in the 
path of tropospheric scattering their amplification is not com- 
pletely realized, and the so-called losses of amplification appear. 
These losses have been shown to be due to the fact that the 
volume of scattering diminishes with increasing amplification 
of the antennas. 

On the other hand, it has been shown that the losses of amplifi- 
cation are due to insufficient correlation of the scattered field in 
the aperture of the receiving antenna. 

We shall show that both these causes of the loss of amplification 
of antennas are identical; i.e., they are different ways of describing 
this phenomenon. 


No. 6, 1960. 


Stability of the Modulus of Amplification of a Feedback Amplifier, 
B. Ia. Lur’e, pp. 615-620. 

Abstract: The noncorrespondence between the formula and the 
versal definition of the magnitude of the sensitivity of a feedback 
system in the well-known book by Bode has been considered. 
The conception has been introduced of the modulus of sensitivity 
of the system, which determines the stability of the modulus of 
amplification of a feedback amplifier. 


Limiting Values of the Resistances of Iterated RC Generators, 
E. N. Garmash, pp. 630-636. 

Abstract: The conditions of the self-excitation of iterated RC 
generators taking account of the ouptut resistance of the amplifier 
have been considered. Formulas have been obtained to deter- 
mine the minimal limiting values of the resistances of the phasing 
circuits. 

Conclusions: 

1 From a comparison of the curves of Fig. 3 for R-parallel and 

'-paralle] networks it follows that, other conditions being equal, an 
RC generator with a C-parallel network is excited at lesser values 
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of R/ Ra than with the network #-parallel. 

2 Beginning with an amplification factor equal to 38 and 
more, an RC generator with a four-link network F-parallel be- 
comes excited at greater values of R/R, than an RC generator 
with a four-link network C-parallel. 

3 For a coefficient of amplification greater than the defined 
value (different for different networks), the resistances of the 
phasing networks can be taken less than the resistance Ry. 

4 Incalculating an RC generator it is necessary to choose such 
a value for the resistance R that, with the given coefficient of 
amplification, the relation R/Ra is greater than the minimum 
value Rnin/ Re. 


Specialized Calculating Machine for Harmonic Synthesis and 
Analysis, A. M. Zaezdnyi and L. M. Rakhovich, pp. 695-700. 

Introduction: The importance of the device of the Fourier 
series in radio technology and electrical communications is well 
known; however, the place occupied by Fourier series in the 
practice of engineering calculation is far from corresponding to 
the possibilities it offers. It is sufficient to say that in the study 
and calculation of linear systems there is no problem which 
cannot be solved by means of Fourier’s series. The reason why 
the use of this device does not correspond to its possibilities is to 
be sought primarily in the difficulties which are met with in 
carrying out the operations of harmonic synthesis. The present 
state of the technique of calculation renders it possible to remove 
these difficulties completely and to place in the hands of investi- 
gators a completely universal method for the calculation and 
study of linear systems. 

The present article gives a short description of a specialized 
calculating machine called ‘‘Synthesis’’ designed in the Depart- 
ment of Theoretical Radio Technology of the Bonch-Bruevich 
Leningrad Electro-technical Communications Institute, in use 
since November 1958. The machine is intended for harmonic 
synthesis and analysis and its application in scientific-engineering 
practice opens up wide new possibilities. 


No. 7, 1960. 


A Receiving Apparatus Permitting of Statistically Noise-Stable 
Coding, Yu. N. Mel’nikov, pp. 737-742. 

Abstract: The principles of a receiving apparatus are described 
which permits of the correction of different numbers of errors for 
different code combinations. A description of a practical model 
of a receiver working on this principle is added. 

Introduction: The noise stability of communication systems 
can be increased by using codes which reveal and correct errors. 
Hence it it easy to understand the interest which is being taken at 
the present time in the investigation of systems which render 
possible the increase of noise stability by the use of corrective 
codes. Here the immediate interest centers in the pursuit of 
rational methods of constructing apparatus for noise-stable 
coding. The present article is devoted to one of the questions 
which arise in the practical designing of such systems; a new re- 
ceiving apparatus is described for both noise-stable and statisti- 
cally noise-stable coding. 


Analysis of the Stability of Valve Circuits and Transistor Circuits, 
E. V. Zelyakh, pp. 787-804. 

Abstract: A method is described for finding the characteristic 
equation of a linear electric circuit, required for the analysis of its 
stability. The method is based on the representation of the cir- 
cuit in the form of a four-pole network, in which the poles of the 
left- and right-hand ends are either short circuited, or discon- 
nected, or connected straight through. The theoretical founda- 
tion of the method is given, together with examples of its applica- 
tion to valve circuits and transistor circuits. 

Conclusion: In conclusion we would point out an important 
result arising from the theory given here. If a four-pole network 
which enters into the composition of a linear circuit is replaced by 
another four-pole network, having in the general case all the four 
characteristic parameters respectively equal to the four param- 
eters of the first four-pole network in a frequency range of zero 
to infinity, then the characteristic equation of the whole circuit 
remains unchanged. This means that the parameters of the four- 
pole circuit given as functions of the frequency characterize it 
equally both under conditions of forced oscillation and also condi- 
tions of free oscillation. 
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by ABLESTAR and AGENA'B 


with DIMAZINE 


(unsymmetrical dimethyihydrazine, UDMH) 


the storable high-performance fuel! 


There are no launching crews or support equipment miles 
above the earth’s surface where AbleStar and Agena-B begin 
their crucial missions. Yet out of a combined total of 52 oppor- 
tunities to perform, from initial flight tests through 12 July 
1961, these advanced upper stages and their earlier configu- 
rations have racked up at least 44 unqualified successes with 
only 4 known failures involving their propulsion systems. 


These powerful, versatile liquid-fuel vehicles, mated with 
rugged Thor and Atlas boosters, have led the race for space 
with an impressive array of firsts. An AbleStar was the first 
rocket engine to be successfully restarted in space and a prede- 
cessor Able helped send the first re-entry vehicle over full inter- 
continental distance. The Agena vehicles were the first to put 
U.S. orbital payloads into the thousand-pound class and the 
first to send capsules back to earth from orbit. Together, they 
have helped orbit no less than 39 of the 48 U.S. satellites! 
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Both are reliable, restartable and applicable to a variety of 
space missions. Both use Dimazine in their bi-propellant main 
propulsion systems. 


Dimazine provides a near-ideal combination of physical and 
chemical properties for rocket applications. It is remarkably 
stable to heat, contamination (catalytic decomposition) and 
shock. Its non-corrosiveness, low freezing point, storability and 
ease of handling minimize the complexities of pre-firing opera- 
tions. It has excellent characteristics of ignition, combustion; 
fluid-flow and cooling. Density, liquid range and vapor pres- 
sure are also favorable. 


Dimazine can be used with a variety of storable and higher- 
energy cryogenic oxidants. Fast expandability of production, 
to meet new and larger requirements with increasingly favor- 
able supply economics, has been repeatedly demonstrated. 
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‘meet these demands, Aerojet 
Qeaility a variety of advanced liquid ro 
ABLESTAR, the first rocket engine 
space...simple pressure-fed engine 
n-board propulsion ...pulse rockets 


nto, Californ 


tomorrow's spacecraft seek the master 
be severe, imposing stringent require 
control, variable t t, and multiple re 
 Aerojet's lea in liquid rocket poy 
Engineers, scientists: investigate outstanding opportunities at Aerojet. 


